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PART Il 
DAMAGE CRITERIA Qi 


Chapter 9 


INTRODUCTION TO DAMAGE CRITERIA 


|| Part J of this manual describes the basic 
phenomena associated with a nuclear explosion 
for various burst conditions. These phenomena 
include: blast and shock. thermal radiation, 
X-ray radiation. nuclear radiation. transient radi- 
ation effects on electronics (TREE). electromag- 
netic pulse (EMP) phenomena, and phenomena 
affecting electromagnetic wave propagation. Part 
1] treats the mechanisms of casualty production 
and damage to military targets, and describes the 
response of these targets by correlating the basic 


physical phenomena. with various defined 


degrees of damage. 


SECTION I 


& CONTENT AND LIMITATIONS 
OF PART I 


9-1 Introduction to Chapter 9 nd 


The information in Part I is divided ac- 
cording to the phenomenology, with one chap- 
ter being devoted to each of the seven phenom- 
ena listed above. This chapter provides a general 
discussion of the physical damage mechanisms 
associated with each of the first six phenomena 
listed above (the last phenomenon does not pro- 
duce physical damage: the degradation of signals 
from radio and radar systems is discussed in 
Chapter 17). Separate sections describe the dam- 
age mechanisms associated with each of the six 
phenomena. In addition, Section IV of this 
chapter discusses the degradation of equipment 


by thermal radiation that might affect its re- 
sponse to the blast wave. 


9-2 Organization of Part il = 


Part II of this manual! is divided into 
chapters according to types of targets that ex- 
hibit similar response characteristics. Further 
subdivisions within the chapters separate the 
types of targets into subtypes, and for each sub- 
type there are frequently further subdivisions 
according to the phenomena that cause the 


damage. 
| The data presented here are interpreta- 


tidhs of complex results of the nuclear weapons 
effects research and test programs of the Depart- 
ment of Defense. A constant effort is made to 
deduce theoretical models and scaling laws for 
the various weapons effects that permit a quanti- 
tative prediction of the extent of a given effect 
from a weapon of one yield related to weapons 
of other yields. Since the initiation of the 
limited nuclear test ban treaty, a large amount 
of effort has been devoted fo the development 
of complex computer codes to predict the en- 
vironments created by the various phenomena 
resulting from nuclear explosions and the inter- 
actions of these environments with personnel 
and military systems. A large number of scaling 
laws presented in Part I that are useful in pre- 
dicting the environment from a given explosion 
were derived from the calculations performed 
with these codes. Many additional scaling laws, 
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derived from field and laboratory experiments, 
and theoretical studies, as well as the codes, are 
presented in the chapters of Part Il to aid in the 
predictions of the response of personnel and 
systems to the environment. Graphical presenta- 
tion is used in preference to tabular or purely 
computational presentations wherever possible. 
The damage curves are drawn for a probability 
of 50 percent of producing the indicated dam- 
. age. Curves of 90 and 10 percent probability are 
included when the quantity and quality of the 
data permit. In addition, Appendix C provides 
the methodology for estimating any probability 
of damage, provided that the median values of 
the response of a given system are known. Esti- 
mates of such median values for the equipment 
discussed in Sections I through II] of Chapter 11 
are provided in Appendix C. Effort has been 
made to provide a comprehensive set of data in a 
readi!, usable form; however, certain categories 
of damage are not amenable to generalization. 
Some limitations in the content of Part II are 
discussed in paragraph 9-3. 
An estimate of the degree of reliability 
_ accompanies most of the data presented herein. 
Statements of the reliability of damage data 
only pertain to the basic response data, which, 
for the target analyst represent the “radius of 
effect.” They should not be confused with the 
terms variability and probability of damage, 
-which pertain to target response; nor do these 
estimates include operational considerations 
such as linear, circular, or spherical aiming and 


fuzing errors, yield variations, and target intelli-’ 


gence. 
93 Limitations in Part Il eS 


As mentioned in paragraph 9-2, the re- 
maining chapters of Part IJ each contain descrip- 
tions of damage to a category of targets that 
exhibit similar response characteristics. Unfor- 
‘tunately, the response to certain of the phenom- 
ena depends so strongly on a specific system 
design that it is impossible to present response 
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data for a “category” of targets (e.g., communi- 
cations equipment, aircraft, missiles). These phe- 
nomena include X-ray radiation and EMP. 


Therefore, general discussions of the types of — 


systems that are likely to be affected by these 
phenomena and the general nature of the 
Tesponses are given in Sections V and VIII of 
this chapter. No response criteria are provided in 
subsequent chapters; however, the potential 


importance of these phenomena are mentioned 


where appropriate. TREE response is-also very 
dependent on specific system design. The re- 
sponse of electronic components and circuits is 


described in Section VII of this chapter. Section | 


IV of Chapter 14 contains a brief review of com- 
ponent and circuit response supplemented with 
discussions of ‘general electrical. responses of 
classes of systems. This latter discussion is 
intended to include a cross section of systems 
that should provide some basis for estimating 
the radiation damage threshold of other similar 
equipment. 


SECTION II 


BLAST AND SHOCK 
DAMAGE @ 


1 | When a blast or shock wave strikes a tar- 


- get, the target may be damaged (distorted suffi- 


ciently to impair usefulness) by the blast or 


_ shock wave itself, by being translated by the 


blast wave and striking another object or the 
ground, or by being struck by another object 
translated by the blast wave. For example, the 
air blast wave can shatter windows, dish in walls, 
collapse roofs, deflect structural frames of build- 
ings, and bend or rupture aircraft panels and 
frames. Vehicles, tanks, artillery pieces, and per- 


‘sonnel can strike other objects on the ground 


while being hurled through the air or tumbled 
along the ground by the blast wave. Ship hulls 
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may be split or crushed by the water shock 2ioatio 
wave. Buried structures or structural founda, Bator ieee 


. py. 
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tions can be displaced. collapsed, or ruptured by . 


the ground shock wave. Usually, the degree of 
damage sustained by a particular target cannot 
be correlated specifically to a single blast or 
shock parameter. The total damage received by 
the target may depend on a combination of air 
blast and ground or water shock parameters. The 
orientation of the target with respect to the 
blast wave. and the type of surface (the topog- 
raphy or the type of soil) associated with the 
target also determine damage. 


Loapine 


The blast loading on an object is a fune- 
tion not only of the blast characteristics of the 
incident wave (rise time, peak overpressure. peak 
dynamic pressure. decay. and duration). but also 
of the size, shape. orientation. and response of 
the object. The influence of the target character- 
istics on the loading is discussed below. with 
emphasis on air blast loading. 


9-4 Air Blast Loading in_the Mach 
Reflection Region 


' The loading on an object exposed to air 
blast is a combination of the forces exerted by 
the overpressure and the dynamic pressure of 
the incident blast wave. The loading at any point 
on a surface of an object can be described as the 
sum of the dynamic pressure. multiplied by a 
local drag coefficient. and the overpressure after 
any initial reflections have cleared the structure. 
Since the loading changes rapidly while the blast 
wave is reflecting from the front surfaces and 
diffracting around the object. loading generally 
comprises two distinct phases: the initial diffrac- 
tion phase: and the phase following diffraction 
when the object is completely engulfed by the 
blast wave. This latter phase approaches a steady 
State and usually is referred to as the drag phase, 
because during this phase the drag forces (i.e., 
the forces resulting from the dynamic pressures) 
are the predominant factors in the production of 


a net translational force on the object. The fol- 
lowing discussion of the loading process is based 
on an ideal blast wave as described in Section I. 
(Figure 2-1), Chapter 2.. Where nonideal blast 
waves, with slow rise time, irregular shapes, and 
high dynamic’ pressures (paragraphs 2-21 . 
through 2-31, Chapter 2) introduce complica- 
tions into the loading process. further explana- 
tion is provided. The loading on an object can be 
described conveniently in three. parts: diffrac- 
tion loading. drag loading. and net loading. 
These are discussed separately below. 

Diffraction Loading. The side of an 
object facing the shock front of an air blast wave 
bears overpressures several times that of the inci- 
dent overpressure because it both receives and 
reflects the shock. In the Mach reflection region 
the overpressure incident on the object is actual- 
ly that of the original free air blast wave which 
has been reflected from the ground surface to a 
higher value. The reflection off the object there- 
fore constitutes a second reflection process. In 
the regular reflection region, the incident over- 
pressure is that of the free air blast wave (see 
paragraph 9-6). The magnitude of the reflected 
overpressure depends on the angle between the 
shock front and the face of the object. the rise 
time of the incident blast wave. and the initial 
incident shock strength. The greatest reflected 
overpressures occur when the direction of propa- 
gation of the shock front is normal to the face 
of the object. when the rise to the peak over- 
pressure is essentially instantaneous. and when 
the incident shock strength is high. As the blast 
wave progresses it bends or diffracts around the 
object, eventually exerting overpressures on all 
sides. Before the object is entirely engulfed in 
the pressure region, however, overpressure is 
exerted on the front side of the object, whereas 
only ambient air pressure exists on the back 
side. During the diffraction phase this pressure 
differential produces a translational force on the 
object in the direction of blast wave propaga- 
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tion. When the blast wave has completely sur- 
rounded a small object, the translational force 
that results from diffraction loading is reduced 
essentially to zero, because the pressures on the 
front and on the back are almost equal. In the 
cases of long objects or short duration blast 
waves, the net force may actually reverse, be- 
cause the overpressure on the front face may 
decay to a value lower than that on the rear 


face. The importance of this translational load- 


ing in the production of damage to the target 
depends on the duration of the loading or on the 
time required for the shock front to traverse the 
target and. therefore, on the size of the target. 
The effects of.the translational load decrease as 
the duration of the load is decreased until, in 
some instances, translational load effects can be 
ignored. The overpressures continue on all sides 
of the object until the positive phase of the blast 
wave has passed. These pressures may be suffi- 
cient to crush an object (a 55-gallon drum may 
be damaged in this manner in addition to dam- 
age that might be incurred by translation). Thus, 
the diffraction phase translational loading de- 
pends primarily on the object size, pressure 


pulse duration, and on increases in differential - 


overpressures resulting from reflection on the 
front face. 

Ss Drag Loading. During the time of dif- 
fraction and until the blast wave has passed, the 
wind behind the shock front causes dynamic 
pressures to be exerted on the object as drag 
loading. Except in the case where shock 
strengths are high, these pressures are much low- 
er than the reflected overpressures; however, 
they produce a translational force that the target 
component receives for the entire positive phase 
duration of the blast wave. For a given blast 
wave, the loading that results from dynamic 
pressure depends principally on the shape and 
orientation of the object, ranging from less than 
four-tenths of the dynamic pressures in the case 
of a cylinder (when normal to the cylindrical 
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axis), to over twice the dynamic pressures for an 


irregular, sharp-edged object. 
- Net Loading. Net loading is the com- 
ine 


b load on the element that tends to translate 
it in the direction of propagation of the blast 
wave. Thus, it is the difference between the load 
on the front face and the load on the back face; 
the loads on the sides are of no effect in produc- 
ing translation. 


Net Loading - 


4 The net load on a target can be develop- 
ed most simply by considering the idealized case 


shown in Figure 9-1, in which a classical, sharp 
fronted blast wave moving along the surface of 
the ground encounters a simple, rigid, fixed 
cube. When the blast wave arrives at the front 
face of the target, this face experiences a sudden 
rise in overpressure to a value p, that is greater 
than the peak overpressure Ap of the incident 
blast wave (also frequently called “‘side on” 
overpressure). As the shock wave moves over the 
cube, subjecting its top and side faces to side on 
overpressure, the pressure on the front face of 
the cube begins to drop as the result of rarefac- 
tion waves that are generated at the edges of the 
front face, and which move across that face. 
When the wave encounters the back face of the 
cube, that face experiences a gradual pressure 
increase. The air behind the front of the incident 
shock wave is in motion, and as the wave en- 
velopes the cube this air motion also affects the 
pressures experienced by the various faces of the 
cube. — 
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After a time that is related to the cube 


a’ and to the velocity of the incident | 


shock wave, the pressure on the front face be- 
comes 
Py = Ap(t) + C(t) 
where 
Ap(t) = the overpressure in the incident 
shock wave as a function of time, f, 


AiR SHOCK 


Figure 9-1, 


Ca, = the drag coefficient of the front 
face of the cube, 

the dynamic pressure in the incident 
shock wave as a function of time, f. 


gtt) 


After the shock wave has passed over the cube 
and has reached the back face, the pressure on 
that face begins to rise and, at a later time also 
related to the cube dimensions and the velocity 
of the incident shock wave, the pressure 
becomes , 


P, = Aplt) ~ Cy,qft) 


where 


Cy, = the drag coefficient of the back face 


of the cube. 


Note that both expressions contain overpressure 
and dynamic pressure components. As far as the 
net horizontal (translational) loading on the 


TARGET 


Initial Conditions for Loading 
of a Rigid, Fixed Cube 


cube is concerned (i.e., the pressures tending to 
move the cube to the right in Figure 9-1), the 
overpressure contributions must be subtracted 
from one another (pressure on the back face 
tends to move the cube to the left) while the 
dynamic pressure contributions must be added - 
to_each other. 

It is convenient to consider the two con- 
tributions separately. Figure 9-2 illustrates the 
overpressure loadings only (dynamic pressure 


_ contributions are not included) experienced by 


the front and back faces of the cube both during 
the early (diffraction) phases — up to time ft, on 
the front face and time ft, on the back face of 
the cube — and after the times that the equa- 
tions given above apply. 


Laan The incident shock wave is of the clas- 
sical (p 


eaked) form with peak overpressure and 
overpressure duration of Ap and rf respectively; 
t, is the time after which the total front face 
loading is represented by the equation for Py 
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Figtire 9-2. | Target Loading by Overpressure ake 


above: 7, is the time for the shock wave to reach 
the back face, and 7, is the time after which the 
total back face loading is represented by the 


e 


uation for p, above. 
« From Figure 9-2b. the total impulse on 


the ront face. 7, (i.e.. the area under the curve 
multiplied by the area of the cube face. A) is 
equal to the sum of the areas a, b, and c. 


I= latb +o 


r 
The total impulse on the back face. J,. is (Fig- 
ure 9-2¢) 

i, =(d+ eA 
The net overpressure impulse is /,- /,. If the 


‘cube is not too large compared to the duration 
of the blast wave. the following approximations 
hold among the areas: 


b d +f, 


Thus. the net overpressure impulse is 


I, = (a + f)A. 


If the initial drop of pressure on the front face 
and the build up of pressure on the rear face are 
assumed to be linear. the net impulse is 


be = A[l/2 (2, - Sp}t, + 1/2 Ap (3 -t3}] 


in which the term containing r, is area “a,’’ and 
the term containing t, and 1, is area ‘f” of 
Figure 9-2. Note in particular that the only 
times are f,, f,, and 7, Le.. those times related 
to the initial envelopment of the cube by the 
blast wave. /,, does not, in this SuipIE analysis, 
depend on Sverpressure duration i 


. ty 


Substitution of commonly accepted 
.t,. and 1, for a curve leads to the 
following expression for the net impulse from 
overpressure diffraction around the target: 


I =A E (37, + ar)| 


where 


¥ = one-half the width or the height. which- 
ever is smaller, 
U = velocity of the shock wave. 


These values of 7,. t, and ¢, originally were 
derived from early shock tube experiments and 
full scale nuclear tests on structures at low pres- 
sure levels. Consequently. they are used for illus- 
trative purposes only in this simple analysis. 
More recent work with two and three dimen- 
sional models in the shock tube has shown that 
the exact values for these times depend strongly 
on model shape and pressure level as well as 
sound velocity in the reflected pressure region 
on the face of the model rather than shock front 
velocity. 

The impulse on the cube due to dynamic 
pressure fe is -_ 


t+ 
IA i gttiCydt 
Qa 


q 
where 
f = drag impulse. 
q = dynamic pressure, 
_C, = drag coefficient of the entire object. a 


combination of Car and Cap 


L 4 By expressing p, in terms of Ap and 
shock strength £, and U in terms of the sound 
velocity. c and the shock strength &. the drag 
impulse equation may be solved in terms of Ap, 
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t, C,, and shock wave duration* ¢*. The com- 


bined impulse is 
-1/2 
Pt (SE ( + 8) 


2c ee 


0.54ApC, ft € 
a Ee 


The first expression in brackets, which is inde- 
pendent of duration, is the overpressure con- 
tribution. The second expression is the contribu- 
tion of dynamic pressure. 

If height of burst and ground distance 
a aled as the cube root of the yield W for 
weapons of different yields, Ap remains con- 
stant, but the shock wave duration f* varies as 
the -ube root of the yield. Thus, the yield de- 
penaznice is 


I, =A E + c(wir)| 


where B, the overpressure contribution dy fora. 


1 kt yield), and C, the dynamic pressure contri- 
bution (. for a 1 kt yield), are both constants. 
Thus, the contribution to total impulse from 
overpressure remains constant, while that from 
dynamic pressure increases as the cube root of 
the yield. For very low fractional kiloton yields, 
the loading is highly impulsive with most of the 
load coming from the overpressure contribution. 
As the yield increases, at a constant scaled 
height of burst and ground distance, the total 
impulse also increases, with an increasing frac- 
tion resulting from dynamic pressure. 

Figures 9-3a and b show representative 
n adings for two classes of weapon yield and 
two classes of structures, one small and one 
large. A small element would be about the size 
of a telephone pole or a jeep; a large element 
would be the size of a house or larger. Since the 
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reflected overpressure is more than twice the 
incident pressure on the front face of the ele- 
ment, the loading displays an initial peak value. 
The reflected pressure decays or clears the front 
face at a time that depends on the size of the 
element. The rapid decay for the small element 
may make the reflected pressure spike of no sig- 
nificance, whereas the slow decay for the large 
element creates a load that may govern the re- 
sponse of the target entirely. For the representa- 
tive cases shown, the diffraction phase 
terminates at time fj i, the time at which the 
reflected pressure has decayed to the incident 
pressure. At this time the drag phase begins. It 
continues until the end of the positive phase of 


the incident blast wave. The load during the drag © 
. phase is shown to be equal to the dynamic pres- 


sure, i.e., the drag coefficients of the elements 
are equal to 1.0. The characteristics of the target 
element determine whether the response of the 
element is governed primarily by the diffraction 
phase or the drag phase. Figures 9-3a and b show 
that for medium and high yield weapons and 
small elements, a much greater impulse (the area 
under the loading curve) occurs during the drag 
phase than during the diffraction phase. As the 
yield increases the drag phase impulse increases 
in importance. For large elements and large yield 
weapons, the diffraction phase and drag phase 
impulses are about equal. In this latter case. the 


drag phase impulse may still be of no impor-— 


tance, because the significant target response 
may occur during the diffraction phase. The dif- 
fraction. phase impulses are not changed by the 
yield of the weapon (this is true for alJ but very 
large structures exposed to low yield weapons), 
whereas the drag phase impulses are directly re- 
lated to the weapon yield (for the same peak 
dynamic pressures). 


5 In this solution it is assumed that the difference between 24 
an % may be negiected, and » = % = ft, 
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Figure 9-3. 1 4 Net Blast Loading on Representative Structures eS 
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9-6 Regular and Mach Reflection ww 


Specific aspects of the blast wave propa- 
gation must be considered when computing the 


load on a target. The loading of a surface target 


in the regular reflection region (pargraph 2-18, 
Chapter 2) is complicated by the vertical compo- 
nent of the incident blast wave, which causes 
multiple reflections between the ground and the 
target and additional reflected pressures on hori- 
zontal surfaces. In the Mach reflection region 
(paragraph 2-18), the loading is simplified be- 
cause the blast wave propagation is horizontal. 
Near the surface of the ground, the vertical com- 
ponent of. the drag forces in the regular reflec- 
tion region is cancelled by the reflected wave 
rapidly; therefore. the brief vertical drag loading 
is ignored. except when the target is near the 
ground zero of an air burst. For aircraft in flight, 
the l-ding may be a single horizontal shock 
from u Mach stem or two separate shocks; the 
first from the free air wave and the second from 
the ground reflected wave. In establishing the 
damage curves for surface targets, the loadings 
on targets in the regular reflection region during 
the diffraction phase are considered separately 
from the loadings on similar targets in the Mach 
reflection region. The surface conditions are as- 
sumed to be average unless otherwise indicated 
in the figures that provide the damage curves. 


Objects that are primarily susceptible to hori- ~ 


. zontal drag loading in the Mach reflection region 
may become primarily susceptible to crushing 


action if they are in the early regular reflection 


region. 


9-7 Target Motion a 


When air blast loading is considered, the 
movement of the target component during load- 
ing is assumed to have negligible effect on the 
- loading itself. Aircraft and missiles in flight are 
exceptions. Their speed, orientation, and move- 
ment during loading assume increased impor- 
tance (see Chapters 13 and 16). 
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9-8 Nonideal Waveforms 


As discussed in paragraphs 2-31 and 
2-32, ideal waveforms are seldom found along 
the surface for overpressure levels above 6 psi. 
The description of the diffraction and drag 
phases given in preceding paragraphs is not true 
in regions of nonideal waveforms. If the over- 
pressure wave has a long rise. time (30 to 40 
msec) to a peak value, full reflection of the wave 
from the surface of a. structure will not occur. © 
At the same time, the relationship between 
dynamic pressure and overpressure is different 
from that described for the ideal blast wave. 


‘During the diffraction phase, the drag forces 


caused by high dynamic pressures may be. the 
predominant damage-producing criteria. Since 
many conventional surface structures sustain 
severe damage at low. peak overpressure levels, 
and nonideal waveforms occur only in the higher 
overpressure regions, such waveforms were not 
considered in determining damage criteria for 
such structures. Careful consideration must be 
given to nonideal waveforms when assessing the 
vulnerability of protective structures that were 
designed to withstand high overpressures. The 
dynamic pressures for such waveforms generally 
will be higher than would be expected if the 
blast wave were ideal. There are few data on 
blast loading in the high pressure regions. 


9-9 Underwater Shock Wave 
Loading 


The air blast wave operates in a com- 
pressible fluid, while the water shock wave oper- 
ates in a noncompressible fluid. This difference 
in medium accounts for the differences in detail 


_betweeen air blast and water shock loading. The 


peak values in water are higher than they would 
be at the same distance from an explosion of the 
same yield in air. However, the duration of the 
shock wave in water is shorter than in air. When 
the shock wave strikes a rigid, submerged sur- 
face, such as the sea bottom, it is reflected. 


When the shock wave reaches the air-water sur- 
fuee (a less rigid medium). however, a refraction 
(or negative pressure) wave occurs. The com- 
bination of the surface reflected shock wave and 
the direct shock wave produces a sharp decrease 
or “cutoff” in the water shock overpressure. as 
shown in Figure 9-4. Incidence of a shock wave 
on a Ship or structure produces damage through 
the direct mechanism of overpressure (or excess 
impulse) and by imparting a transitional velocity 
to the target structure. , 


9-10 Ground Shock Loading Qi 


The Joading of buried structures by 
ground shock is connected intimately to the re- 
sponse of the structures. For certain under- 
ground structures. serious damage will only 
occur if the ground shock is so intense that the 
damage area for those structures is confined 
closely to the crater area of a surface or under- 
ground explosion. For structures near the sur- 
face. air blast induced ground shock may cause 
significant damage. The damage to these struc- 
tures is more closely related to air blast pressures 
than it is to crater dimensions. Loading pressures 
are numerically equal to the ground stress nor- 
mal to the structure. Such pressures do not pro- 
duce detectable reflected pressures. Internal 
equipment of a structure may be subjected to 
ground shock accelerations that will severely 
damage the equipment without damaging the 
structure. For a discussion of accelerations re- 
sulting from ground shock see Section I]. Chap- 
ter 2. and Sections II and II]. Chapter 11. 
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Le Damage to a target is closely related to. 
and is a direct derivative of. its response. For 
targets anchored to the ground. damage is usual- 
ly the result of displacement of one part of the 
target with respect to another part, resulting in 
permanent distortion, collapse, or toppling. For 
movable targets, however, the target may be 


moved by the loading with or without a result- 
ing damage. In the latter case, the damage to the 
target is governed primarily by the manner in 
which the moving target comes to rest. Whether 
drag phase loading or diffraction phase loading 
causes the greater damage will depend upon the 
weapon yield, target characteristics. and the 


' damage level considered. 


9-11 Surface Structures ee 


The predominant cause’ of failure to 
large targets. such as buildings that have small 
window areas compared to wall areas that either 
support the structure or are as. strong as the 
structural frames. is the pressure differential 
between the front and rear faces that exists for a 
relatively long period of time. If the window 
area is large, the pressure on each wall is equal- 


-ized quickly by the entry of the blast. wave 


through the windows. The pressures exerted on 
the inside of the wall reduce the translational 
force on the wall. This translational force also is 
reduced because of a smaller wall area on which 
the pressures can act: however, the force exerted 
on interior partitions and rear. walls tends to off- 
set the reduction in front face loading in the 
production of total damage. When the overpres- 
sures causing translational force on the struc- 
tural component are equalized quickly as a 
result of the geometry or construction of the 
building. the primary damaging forces are those 
that are significant in damaging structural com- 
ponents that have fairly small cross sections. 
such as columns and beams. Structures that are 
normally damaged by drag forces include smoke 
Stacks, telephone poles, truss bridges, and steel 
or reinforced concrete frame buildings with light 
walls. These buildings are drag sensitive, because 
the light walls of corrugated steel. asbestos, or 
cinder block fail at low reflected pressures. and 
they do not transmit a significant load to the 
structure frame. Only the frame is exposed to 
the blast, and, being composed of small cross 
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9-12 


section structural elements, it is distorted pri- 
marily by drag forces. These buildings are not 
considered severely damaged unless the struc- 
tural frame has collapsed or is near the point of 
collapse. A tree is a good example of a drag 
sensitive target, because the duration of the dif- 
fraction phase is extremely short and there is 
considerable force applied by the high wind 
velocity drag loading. Most military field equip- 
ment is drag sensitive. because damage generally 
results from the tumbling or byesuLoine caused 
by the drag forces. 


9-12 Shielded Structures , | 


If.the target is shielded from the drag 
forces. or if it lies within the early regular reflec- 
tion region. high overpressures may become the 
damage producing criteria. For blast resistant 


aboveground structures designed to resist more — 


than 5 to 10 psi overpressure. the distinction 
between diffraction and drag sensitivity cannot 
be defined well. because full reflection from the 
surface of the structure does not occur. and 
dynamic pressures exceed those expected in the 
case of the incident waveform. As a result, drag 
forces may be predominant in producing dam- 
age, even during the diffraction piave? 


9-13 Aircraft 1 f 


Aircraft may be damaged by the forces 
developed in the diffraction phase, in the drag 
loading phase, or in both. Parked aircraft can 
receive light. crushing forces corresponding to 
low overpressures. For example, light skins and 
frames are easily dished and buckled at relatively 
low overpressures. At higher overpressure levels, 
drag loading (referred to as “gust loading” with 
respect to aircraft) adds to the damage. At these 
levels, much of the damage may result from 
translation and overturning of the-aircraft. For 
aircraft in flight, the diffraction and drag forces 
combine with the existing aerodynamic forces to 
develop destructive loads on airfoils at low over- 


pressure levels. The diffraction or crushing over- 
pressure effects on the fuselage and other thin 
skinned components, however, are usually of 
secondary importance for most in-flight aircraft. 


Responses of aircraft are discussed in more de- 


tail in Chapter 13: 


SECTION III 


| | THERMAL RADIATION 
DAMAGE 


INTRODUCTION. GR 


Two important effects of thermal radia- 
tion are injuries to personne! (burns) and fires 
that might be ignited in the target area. Burns to 
personnel are treated in Section I] of Chapter 
10. Section VII of Chapter 11 discusses fires in 
urban areas. Section II] of Chapter 15 discusses 
fires in forest stands. The effects of thermal radi- . 
ation on various classes of equipment are de- 
scribed in other chapters of Part JI. This section 
contains a description of the properties of ma- 
terials that might result in ignition or degrada- 
tion of their physical properties. A brief discus- 
sion of survival in fire areas is also provided. 
Section IV describes the degradation of struc- 
tural resistance to air blast of materials, primari- 
ly metals, as a result of thermal radiation. 

w Chapter 3 provides the data necessary to 
estimate the thermal environment. In general, 
thermal damage is likely to be more important 
than blast damage for surface targets for high 
yield weapons that are air bursts rather than sur- 
face bursts. As discussed in Chapter 3, the influ- 


_ence of weather makes thermal effects much less 


predictable than blast effects. Clouds or haze 


can provide a protective screen, and moisture 


from an earlier rain can reduce the number of 
fires in outdoor targets. In some cases, cloud 
cover can enhance thermal effects on the 
ground. The criteria for thermal damage given in 
succeeding paragraphs are representative values 
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for dry materials directly exposed to the thermal — 


source. These values should be applied with the 
understanding that they can be modified by a 
variety of factors discussed briefly in the follow- 
ing paragraphs. 


914 Radiant Exposure | 


The thermal damage inflicted on a target 
depends upon the incident energy per unit area. 
This quantity is called radiant exposure and, as 
described in Chapter 3, it usually is expressed in 
calories per square centimeter (cal/cm?). Radi- 
ant exposures below 2 cal/cm? will produce 
little damage other than possible eye injury (Sec- 
tion Il, Chapter 10). Radiant exposures above 
10 cal/em? usually produce significant damage 
in unprotected target areas. 


9-15 Thermal Pulse Duration Ga 


The radiant exposure required to dam- 
age a particular target varies with the duration of 
the thermal pulse. The reason is readily seen 
from the following example: about 3 to 5 
_cal/cm? from a short thermal pulse, e.g., from a 
1 kt detonation, will produce a second degree 
burn on bare skin, but direct sunlight delivers 
the same radiant exposure in a little over 2 min- 
u ith no serious effects. 
“ The time scale of the thermal pulse from 
low altitude nuclear weapon bursts may be char- 
acterized by the parameter t Saas the time dur- 
ing the final pulse at which the fireball is radiat- 
ing maximum power. This time increases with 
increasing yield and decreases with increasing al- 
titude. The relation of t,,, to pulse duration 
and the method for calculating ¢,, ,, are given in 
Chapter 3. As described in Chapter 3, the time 
scale of the thermal pulse may be specified indi- 
rectly in terms of yield for low altitude bursts. 
The thermal pulse for high altitude bursts (up to 
about 100,000 feet) may be described in terms 
of an equivalent sea level burst. Specifying the 
thermal pulse duration in terms of yield is con- 
venient since it eliminates the necessity to calcu- 
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late toes therefore, many damage criteria are 
given in terms of radiant exposure and yield 
rather than radiant exposure and time to final 
maximum. _- 


9-16 Target Response (ae 

For most materials, the heat absorbed 
from the thermal pulse initially is confined to a 
thin surface layer. Damage usually results from 
high surface temperatures. The nature and the. 


degree of the damage depend not only on the 


intensity and duration of the thermal pulse, but 
also on several properties of the material that are 
described below. : 

Thickness. Thick organic materials such 
as wood, plastics, and heavy fabrics char and 
may burst into flame while exposed to the 
thermal pulse; however, this flaming is only a 
transient effect. As the radiant pulse decays, the 
absorbed thermal energy continues to penetrate 
the material. This flow of heat allows the surface 
to cool to a temperature too low to support 
combustion. es 

Thin materials, such as light fabrics, 
newspaper, dry leaves, and dry grass tend to be- 
come hotter than the surfaces of thicker ma- 
terials since the absorbed thermal energy is con- 
fined to a relatively small volume of material. 
Also, since these materials are heated through- 
out, they usually continue to burn once they are 
ignited. Moreover, subsequent arrival of the blast 
wave frequently will fail to extinguish the fire. 
Thermal Conductivity. Most. metals 
allow rapid penetration of thermal energy to 
depths below the surface as a result of their high 
values of thermal conductivity. Consequently, 
the peak surface temperatures induced in thick 
pieces of metal are considerably less than the 
peak surface temperatures induced in most non- 
metallic materials. Thin pieces of metal may fail 
as a result of the combined effects of the ther- 
mal pulse, which reduces the strength of the 
metal by heating it, and the blast wave, which 
causes it to bend or collapse. This effect is dis- 


Se age es 


cussed in more detail in Section IV. Dry rotted 
wood (or punk) has .a lower thermal conductiv- 
ity than sound wood as a result of its porous 
structure. When the surface of punk is ignited by 
a thermal pulse. conduction of heat to the in- 
terior is too slow to allow significant cooling of 
t irface. and burning continues. 

Color. Light colored objects of a given 
thickness are more resistant to thermal radiation 
than dark colored objects. because they reflect 
more of the incident energy. Color has little ef- 
fect on the response of materials that blacken 
(char) early in the thermal pulse. because the 
energy delivered during the remainder of the 
pulse is absorbed efficiently by the charred sur- 
face. 


Transparency, Transparent materials are 
relatively resistant to. thermal damage (even 
though transpurency in the visible region does 
not assure that the material is transparent to 
infrared) because the radiant energy that passes 
through them usually does not contribute to 
heating.. Partially transparent materials are ther- 
mally resistant -because the incident thermal 
energy is deposited over a range of depths rather 
than being confined to a thin surface layer. Peak 
temperatures induced in partially transparent 
materials therefore are likely to be lower than 
the peak temperatures produced at the surfaces 


of opaque materials. 
Moisture Conteur. Thermal damage to 


materials that absorb moisture depends on the 
percentage of water in such materials. Usually, 
the moisture content varies with the prevailing 
relative humidity. Exposure to recent rain, how- 
ever. may alter the moisture content sighificant- 
ly. Scorching or charring of an organic surface 
by radiant energy is preceded by vaporization of 
the water. Consequently. more energy is requir- 
ed to produce a given damage effect on wet sur- 
faces or on targets in highly humid atmospheres. 
Materials located indoors and exposed to ther- 
mal radiation through windows are damaged 
more readily during the latter part of the heating 


season (late winter and early spring), largely be- 


cause of decreased interior humidities. 


9-17 Target Orientation , | 


In clear atmospheres and in the absence 
of reflecting surfaces, the amount of energy inci- 
dent on a unit area of a target surface is greatest 
when the surface is directly facing. the burst. If 
the target surface is not facing.the burst. the 


‘energy per unit area depends on the angle be- 


tween the perpendicular to the surface and the 
direction of the incoming radiation. — . 

When the atmosphere is hazy, much of 
the energy received at the target is scattered by 
atmospheric particles, arrives from all directions. 
and reaches portions of the target that are not 
exposed to direct radiation. Scattered energy is 


- likely to be a large fraction of the total energy 


received when the slant range from the burst to 
the target exceeds about half the visual range. 
Reflection from clouds and from the ground can 
produce similar effects. 


9.18 Shielding QB 


Any object that casts a shadow is capa- 
ble of shielding objects behind it from the direct 
component of thermal radiation. Trees, build- 
ings. foxholes. hills, etc. offer effective protec- 
tion except when the scattered component of 
thermal radiation is large. After leaves have been 
stripped by the dynamic pressure (wind) of a 
nuclear burst. the trees offer little shielding from 
the thermal radiation produced by. a second 
burst. Reflection of thermal radiation from ex- 
posed walls of foxholes is about 5 percent. 


THERMAL RESPONSE OF 


MATERIALS (ay 


The amount of thermal radiation that a 
material will absorb depends on the properties 
of the material discussed in paragraph 9-16. The 
amount. of energy absorbed usually is deter- 
mined by multiplying the incident energy by an 
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absorption coefficient. The absorption coeffi- 
cient cannot be measured directly, but a good 
estimate can be made by performing spectral 
measurements of reflectance (the fraction of 
incident energy reflected from the surface) over 
the range of wave lengths of the nuclear spectral 
distribution and by finding an average value of 
reflectance. The coefficient of absorptance is 
equal to one minus the average value of the re- 
flectance. The absorptance determined in this 
manner is valid only as long as the surface does 
not change as a result of the absorbed thermal 
energy. Another approach is to divide the quan- 
tity of energy absorbed by the incident energy. 
The absorbed energy is determined from experi- 
mentally determined temperature vs time rela- 
tionships and appropriate theoretical. energy 
balance relationships. Absorptance is a function 
of time and temperature although an average 
value is usually employed. This function when 
known may be found as an input specification 
for certain computer programs. In the absence 
of more definitive data for a specific systems, an 
absorption coefficient of 0.5 is a reasonable 
value to assume for many materials, particularly 
metals. For a conservative defensive assumption, 
the absorption coefficient may be taken to be 
1.0 (i.e., total absorption), particularly for dark 
porous materials. 

In many cases, it is appropriate to con- 
sider only one or two fuel types, assuming that 
these fuels are the critical items that will deter- 
mine whether fires are started in a particular 
area. In other cases, more detail is required. Fot 
example, an estimate of the thermal energy that 
will ignite a particular type of material used in 
military uniforms may be desired. The data .con- 
tained in the succeeding paragraphs is intended 


as a guide in the solution of more specific — 


problems. 
919 Thickness Effects @ 


The time required for thermal energy to 
penetrate very thin materials is short compared 
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to the duration of the thermal pulse. At the end 
of the pulse, the absorbed thermal energy is dis- 
tributed more or less uniformly throughout the 
material. In very thick materials, the end of the 
pulse finds most of the absorbed thermal energy 
in a surface layer, with the bulk of the material 
a t unaffected. 

If the terms thermally “thin” and ther- 
mally “thick” are used to distinguish between 
materials that are heated throughout from those 
that are initially heated only at the surface, it is 
apparent that physical thickness is not the only 
parameter involved. For example, as a result of 


the rapid penetration of heat into metals, a sheet 
of metal is more likely to be thermally thin than 
is a sheet of insulating material of the same 


p sical thickness. 
The ability of a short pulse of energy to 


penetrate a target material is most readily. mea- 


‘sured in the laboratory and most readily treated 


analytically if the thermal pulse is given a rec- 
tangular waveform rather than the more com- 
plex waveform produced by a nuclear weapon. 


_A parameter that is useful for calculating ther- 


mal response of materials is the characteristic 
thermal response time 7,, given by the equation 


° 


7 = pC, L? /k sec, 


where k is thermal conductivity (cal-sec! 
cm!°C?!), pC. is heat capacity per unit volume 
a rp 4 i ; 

(op = density in g-cm’” and C. = specific heat at 
constant pressure in cal-g"! C Vy, and Z is the 
thickness, in centimeters, of the layer of 
material. 

The quantity 


is called thermal diffusivity (cm/sec). Use of 
this quantity simplies the previous equation to 


L- 
T. =— sec.” 
< a 


Thermal diffusivity and other properties of a 
n er of materials are shown in Table 9-1. 

Characteristic time may be related in 
several ways to the thermal response of a target. 
Most important is that the ignition of thin fuels 
by a rectangular thermal pulse requires the least 
radiant exposure when pulse duration is. about 


equal to 7,. 
= Two other relations apply to a thick slab 


of material. For any particular material exposed 
to a rectangular pulse of length 7. the previous 
equation can be transformed to-give a character- 
istic thickness 


& = \/oar cm, 


for which the characteristic time is equal to the 
pulse duration. If a thick slab of this material is 
exposed to a pulse of length 7. the temperature 
rise at the surface is the same as would be pro- 
duced by uniformly distributing the absorbed 
thermal energy in a slab of thickness 5, and the 
peak temperature rise at depth 6 in the thick 
slab is about half as great as the peak tempera- 
t ise at the surface. 

For example, consider a block of red 
pine that is exposed to 15 cal/em? from a rec- 
tangular pulse of 3 seconds duration. From 
Table 9-1. the properties of red pine are 


p = 0.51 gicm?, 
Cie 0.4 cal/g + °C, 
a = 24 x 10% cm?/sec. 


The characteristic depth is 


§ = Var = V(24 x 107)(3) = 0.085 cm. 


The mass of wood per unit area in a slab of this 
thickness is 


pd = (0.51)(0.085) = 0.043 g/cm?. 
The heat absorbed by the wood before it begins 


to scorch is equal to the product of the incident 
radiant energy. Q, and the absorption coeffi- 


cient, A. If the absorption coefficient is assumed 


to be 0.5. 


QA = (1590.5) = 7.5 cal/em?. 


Absorption of this amount of energy in-a layer 
of thickness 6 would result in an energy density 
of . 


= ———— = 174 cal/g. 


If the energy were evenly distributed through 
this layer, the resulting temperature rise would 
be 


and the peak temperature rise at the surface of 
the wood would be about the same. 

es The result obtained above may be gen- 
eralized as follows: 


ar =«@4. 224 __ = —_@4__ 
s pac, pC, 9/ OF pC, Vrkip, 


* This equation is useful, but it is.by no means exact. The 
‘2 heat-flow analysis from which this equation is derived 


“neglects the effects of radiation and convection heat losses from 


the surfaces of the exposed sample. It also assumes an isotropic 
medium. i.e., a medium whose structure and properties in the 
neighborhood of any point are the same relative to all directions 
through the point. It also neglects the changes in thermal prop- 
erties that occur as the exposed material heats, volatilizes, chars, 
and bursts into flame. 
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Table 9-1. @ Thermal Properties of Materials Ai 


Specific 
Density, p Heat, C Conductivity, k ' Diffusivity, a 
’ Materials (gm/em?) (cal/gm + °C)  .(cal/sec * cm * °C) (cm?/sec) 
Insulating Materials 
Air ' 946 x 107 0.24 0.55 x 104 0.22 
Asbestos 0.58 0.20 4.6 x 104 40. x 104 
Balsa 0.12 04 | 1.2 x.107 - 25. x 107 
Brick (common red) 1.8 0.2 16. x 10% 18. x 104 
Celluloid 1.4 0.35 $0 x 104 10. x 104 
Cotton, sateen, green 0.70 0.35 1.5 x 104 25 x 104 
Fir, Douglas- 
spring growth 0.29 0.4 2. x 104 17. x 104 
summer growth - 1.00 0.4 5. x 104 12. x 107 
Fir, white 0.45 0.4 26 x 104 14. x 104 
Glass, window ray. 0.2 i9. x 10% 43. x10“ 
Granite 25 0.19 66. -x 104 140. x 10 
Leather sole 1.0 0.36 3.8 x 104 11. x 107 
Mahogany 0.53 0.36 3.1 x 104 16. x 104 
Maple 0.72 0.4 , Shes RAO 16. x 1074 
Oak 0.82 04. — - 50 x 104 15. x 107 
Pine, white 0.54 0.33 3.6 x 104 18. x 104 
Pine, red 0.51 0.4 5. eie* 24. x 1074 
Rubber, hard 1.2 0.5 3.6 x 104 60. x 107 
Teak 0.64 0.4 4.1 x 104 16. x 10% 
Metals (100°C) 
Aluminum 35 0.22 0.49 1.0 
Cadmium 8.65 0,057 0.20 0.45 
Copper 8.92 0.094 0.92 1.1 
Gold 19.3 0.031 0.75 | 1.2 
Lead 11.34 0.031 0.081 0.23 
Magnesium 1.74 0.25. . 0.38 0.87 — 
Platinum 21.45 . 0.027, 0.17 0.29 
Silver 10.5 0.056 0.96 1.6 
Steel, mild 7.8 0.11 0.107 1.2 
Tin 6.55 0.056 0.14 0.38 
Miscellaneous Materials 
Ice (0°C) 0.92 0.492 54. x 104 120. x 104 
Water 1.00 1.00 14. x 104 14. x 104 
Skin (porcine, dermis, dead) 1.06 0.77 9. x 104 ll. x 104 
Skin (human, living, averaged  —s-:1.06 0.75 8. x.104 30. x 104 
for upper 0.1 cm) 
Polyethylene (black) 0.92 0.55 8. x 104 17. x 104 
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where A7_ is the peak temperature rise at the 
surface. The parameters that define the thermal 
pulse may be separated from those that define 
the material properties. and. 


ar, = (Fa) yae): 


For a fixed rectangular pulse. Q/\/7 is a con- 
stant. and the equation may be written 


In practice, the surface temperature rise 
produced in a thick material is approximately 
proportional to 4 Vie pC over a wide range of 
pulse shapes and pulse times. Thus. this param- 
eter may be used us a measure of the relative 
susceptibility of various materials to surface 
hegting. 

Direct measurements of the ignition 
properties of a representative cellulosic fuel have 
been reported. The material is a-cellulose. black- 
ened by the addition of carbon. and made into 
sheets of various thicknesses. This material has 
the uniformity required to obtain repeatable test 
results. and its thermal properties are similar to 
those of common cellulosic fuels. such as paper 
and drv leaves, 

Figure 9-5 shows the results of exposing 
this material to the idealized rectangular thermal 
pulse. Figure 9-6 shows the results obtained in 
the more practical case of a simulated-weapon 
pulse. These curves show data for black a- 
cellulose sheets of a variety of thicknesses by 
using normalized coordinates. The ordinate is 
normalized radiant exposure. One unit on this 
scale corresponds to the energy per unit area 
that will raise the temperature of the entire 
sheet 1°C. Thus. a given value of the ordinate 
represents a higher radiant exposure for a thick- 
er sheet of cellulose. The abscissa of Figure 9-5 


is normalized pulse duration, the duration of the 
rectangular pulse divided by 7,, the character- 
istic thermal response’ time of the sheet of ma- 
terial. The abscissa of Figure 9-6 is a similar nor- 
malized time, /,.,,/T,- Normalization fails in 
the case of long pulses, and the curves break into 
families of curves for different thicknesses of 
material. 


If the pulse has such a short duration 
that the exposed surface reaches the ignition 
temperature while the back surface remains 
cool. two distinct ignition thresholds appear. 
Low values of radiant exposure produce high 
surface temperatures and flames, but the flames 
do not persist after the end of the thermal pulse. 
This region is labeled transient ignition in Fig- 
ures 9-5 and 9-6. Sustained ignition. only occurs 
when higher radiant exposures raise the tempera- 
ture throughout the thickness of the cellulose to 
a level that is sufficiently high to sustain the 
flow of combustible gases from breakdown of 
the fuel. Jt is difficult to supply sufficient. 
energy with short pulses, since a large amount of 
the energy that is deposited is carried away by: 


‘the rapid ablation of the thin surface layer. This 


transient flaming phenomenon is typical of the 
response of sound wooden boards to a thermal 
pulse. 


If the pulse is of long duration. the igni- 
tion threshold rises because the exposed material 
can dissipate an appreciable fraction of the 
energy while it is being received. For very long 
rectangular pulses an irradiance of about 0.5 cal > 
em*sec! is required to ignite the cellulose. 
Heat supplied to the material at a slow rate is 
just sufficient to offset radiative and convective 
heat losses, while maintaining the cellulose at 
the ignition temperature of about 300°C. 

Materials respond to the thermal pulse 
fr a nuclear weapon in much the same way 
that they respond to a rectangular thermal pulse; 
however. a number of practical effects make a 
thorough analysis of ordinary fuels ignited by 
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weapon pulses difficult. While the fuel absorbs 
heat and its surface chars, its absorption coeffi- 
cient changes. As a result, the shape of the 
energy pulse actually absorbed by the material 
differs from that of the incident energy pulse 
(this effect is avoided in the a-cellulose by mak- 
ing the material black initially). A weapon pulse 
has no single rectangular wave equivalent, since 
the nature of material response to a thermal 
pulse varies with pulse shape as well as pulse 
length. However, for many calculations a rough 
equivalence may be. obtained by assigning be- 
tween 2.5 and 3 times the irradiance of the rec- 
tangular pulse to the weapon pulse. | 

Both rectangular and weapon pulse data 
s an optimum pulse length for igniting ma- 
terials. For the rectangular pulse, this optimum 
occurs with a pulse length slightly shorter than 


the characteristic thermal response time of the . 


material. A pulse of this optimum duration lasts 
long enough to allow an appreciable amount of 
heat to penetrate to the back face, but it is short 
enough that only a small heat loss occurs during 


the pulse. 
a™ example, consider newsprint that is 


0. cm thick, which has a thermal diffusivity, 


a, of 10° cm?/sec. From the equation given pre- 


viously 
L 2 


70 ee 
From Figure 9-6, the threshold for sustained 
ignition of black a-cellulose is a minimum when 


Imax _ max 


T) L 
The properties of newspaper are similar to those 


of a-cellulose, so the same relation may be 
assumed to hold for the newsprint. Thus, 


L? _ (0.045)(0.008)? 
~ 0.001 


= 0.045. 


0.045 


max a 


2.9 x 10° sec. 


~ 
il 


™ 
" 


From the equation given in Chapter 3. 
tax 7 9.043 yo.43 (lp, 042 sec, 
At sea level p = p,, and 


W = (t,,,,/0.043)?> 


max 


4\23 
ji fe 0.000 KE: 
4.3 x 102 


This yield is so smali that little confidence may 
be placed in the result; however, this example 
demonstrates the general rule that short thermal 
pulses are more likely to ignite newsprint than 
long pulses over the ‘range of yields that 
probably will be of practical importance. 


920 Color [ij | 
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Although thermal damage depends on 
the radiant exposure to which a target is subject- 
ed, it depends in a much more direct manner on 
the amount of thermal energy absorbed. Thus 
color, which indicates the amount and spectral 
distribution of the visible energy reflected, is 
one indication of the resistance of a material to 
thermal damage. 

White materials reflect most of the 
energy in the visible portion of the spectrum 
that is incident on them, but all cellulosic fuels 
(which constitute the commonest potential igni- 
tion sites) absorb energy in the near infrared 
region. Cellulosic materials also char in response 
to high radiant exposures, and once the surface 
begins to blacken they absorb strongly. The 
total amount of energy absorbed depends on the 
absorption coefficient averaged not only over 
the spectrum of energies contained in the ther- 
mal pulse, but also over the interval of time dur- 
ing which the pulse is received. For a thin cellu- 
losic material exposed to a moderately high level 
of radiant energy (just enough to ignite the ma- 


terial) blackening will not occur immediately, 
and the following rules give an effective average 
for the absorption coefficient that is roughly 
correct. 

@ if the material is white or nearly white, the 
effective absorption coefficient for the 
thermal pulse ts about the square root of 
the absorption coefficient for visible light. 
Example: Jf the visible absorption coeffi- 
cient is 0.1. the effective absorption coeffi- 
cient for the thermal pulse is about 0.3. 

@ Jf the material has a color that is as dark as 
or darker than dove gray -(a fairly light 
shade of gray). the effective absorption co- 
efficient for thermal radiation. is about 
equal to the visible absorption coefficient. 

@ The response of most of the nonwhite 
materials. that provide potential ignition 
sites may be approximated by assigning an 
effective absorption coefficient of 0.5. 


9-21 Transparency 


The energy that passes through a partial- 
ly transparent material cannot heat it. Thus. the 


energy deposited in the material is the incident — 


energy minus the reflected energy: minus the 
transmitted energy. 


The depth of penetration of thermal 
energy into a partially transparent. material is 
determined by two mechanisms: heat absorbed 
at or near the surface penetrates by diffusion: 
and heat penetrates as radiant energy before it is 
absorbed. As shown in paragraph 9-19, the diffu- 
sion of heat may be associated with a transient 
thickness 6. Similarly, the penetration of ther- 
mal energy due to diathermancy (partial trans- 
parency) may be associated with a characteristic 
depth of 1/8, where 6 is the extinction coeffi- 
cient. In a practical situation. both of these 
mechanisms work simultaneously, and the heat 
penetrates deeper than it would 7 only one 
mechanism were effective. 


Usually. one or the other of the two 
mechanisms is dominant. If the product £6 
exceeds 3 or 4. penetration will occur princi- 
pally by diffusion. If the product is less than 
about 1/4, penetration will be principally by 


.transmission of radiant energy. Temperature 
_ profiles in a semi-infinite solid are shown in Fig- 


ure 9-7 in terms of dimensionless parameters, 
where AT is temperature rise, k is thermal con- 
ductivity, p is density, C_ is specific heat. H is 
energy absorbed (not the incident energy) per 
unit area, 7 is rectangular pulse duration, and.6 
is characteristic depth defined in paragraph 9-19. 


9-22 Effect of Humidity 


ae The water content of thin fuels responds 
pi to changes in humidity. The radiant 
exposure that is required to ignite a thin cellu- 


- Josic fuel is approximately 


C= 0). Ue h/2), 


where /t is relative humidity and Q, is the radi- 

ant exposure that will ignite the fuel when it is 
etely dry. 

In general, ignition data give radiant 


exposure thresholds for average rather than ex- 


tremely dry conditions. If the tabulated value of — 
radiant exposure Q, applies to a relative hu- 
midity of /,, the equation given above becomes 


1+ h/? 


as ae Tek 


For example, if the radiant exposure 
from a 1 Mt explosion that is required to ignite 
new white bond paper is 30 cal/cm? when the 
relative humidity is 65 percent, the ignition 
threshold for the same explosion when the rela- 
tive humidity is 15 percent is 


1 + 0.15/2 


= s 2 
Q = (30) (} - ps) 24 cal/em-. 
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DIMENSIONLESS TEMPERATURE RISE, AT VipC, /JHV/ 7: 4 


DIMENSIONLESS DEPTH, x/5 


Figure 9-7. we Dimensionless Temperature Profiles in Opaque and 
Diathermanous Semi-Infinite Solids 
Exposed to a Rectangular Thermal Pulse wa 
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9-23 Response ‘of Metal Sheets | 


Both the tensile strength and the elastic 
modulus of metals tend to decrease with increas- 
ing temperature. Thus, thermal energy has the 
potential to weaken metallic structural mem- 
bers. If exposed metal parts are essential to the 
strength of a structure or a piece of equipment. 
weakening the metal by heating it may make the 
unit more susceptible to damage by the blast 
wave , 

Most structural metal is relatively im- 
mune to this form of damage in the loading 
levels of interest. If a heavy part such as a steel 
I-beam is subjected to a thermal pulse that is 
sufficiently intense to weaken it appreciably. it 
is probably so close to the burst that the blast 
wave could demolish the structure without assis- 


tance from the thermal pulse. Walls or roofs 


usually shield heavy structural members from 


the thermal pulse. . 
- Light weight units are more likely to be 


vulnerable to this synergistic damage. The outer 
covering and the materials that give these units 
mechanical strength are frequently one and the 
same. The material is likely to be an aluminum 
alloy. which loses strength at a much lower tem- 
perature than steel. 

(U) The degradation of the response of 
metals to the blast wave as a result of thermal 
exposure is: discussed in more detail in Section 
IV of this chapter. 


Materials 


eS The genetal.trends discussed in the pre- 
ceding paragraphs are illustrated by a number of 
specific examples in Tables 9-2 and 9-3. These 
tables show the critical radiant exposures for 
specified damage to fabrics and to various other 
materials, respectively. The radiant exposures 
are shown for three weapon pulse durations. The 
corresponding yields are shown in the footnote 
for low air bursts (note that altitude scaling is 
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required for a 24 Mt burst to be considered an 
air burst (see Chapter 3)). Various yield/altitude 
combinations that will give the same values of 
tnax may be obtained by the methods described 
in Chapter 3. 

The values presented in Table 9-2 for 
fabrics apply for an ambient relative humidity of 
65 percent and an ambient temperature of 20°C. 
For extremely dry conditions, the values shown 
for fabrics should be reduced by 25 percent. For 
extremely high relative humidities, near. 100 per- 
cent (at 20°C), the values for fabrics should be 
increased by 25 percent. If the fabrics are water- 
soaked. the critical radiant exposures should be . 
increased by 300 percent. The values for uni- 
forms in Table 9-2 refer to damage to the ma- 
terial itself. They are not applicable for predict- 
ing skin burns under uniforms. However, when 
the outer garment ignites. the probability of skin 
burns should be considered. 


The effect, “tears on flexing,” indicates 
the radiant exposure level at which the outer 
garment ceases to sérve as-a reliable barrier to © 
the environment. Under certain environmental 
conditions the garment would require immediate 
replacement. . 


Cotton or rayon fabrics generally can be 
ignited by exposure. to thermal radiation. 
Nylons, dacrons, and similar synthetics generally 
are not ignited but melt at relatively low radiant 
exposures. Woolens do not ignite, but the fibers 
scorch and coalesce. When a sufficient thickness 
of the woolen fabric is scorched it becomes 
brittle and can be torn or crumbled readily. 


Most thick, dense materials that ordinar- 
ily are considered inflammable do not ignite to 
persistent flaming ignition when exposed to 
transient thermal radiation pulses. Wood, in the 
form of siding or beams, may. flame during the 
exposure but the flame is extinguished when the 
exposure ceases. 
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Table 9-2. zs Approximate Radiant Exposures for Ignition of Fabrics | | 


Radiant Exposure* 


(cal/cm?) 
Weight Effect tax max ae 
Material (oz'yd*) Color on, Material G.2 sec 10 sec 3.2 sec 
Clothing Fabrics 
Cotton 8 White Ignites 32 48 85 
Khaki Tears on flexing 17 27 34 
Khaki Ignites 20 30 39 
Olive Tears on flexing 9 14 21 
Olive Ignites 14 19 2 
Dark blue Tears on flexing li 14 17 
; Dark blue ‘Ignites 14 19 21 
Cotton corduray 8 Brown Ignites i 16 ae 
Cotton denim. new 10 Blue Ignites 2 2? 44 
Cotton shirting 3 Khak: Ignites 14 21 28 
Cotton-nylon mixture 5 Olive Tears. on flexing & 15 17 
5 Olive ignites 12 28 53 
Wool & White Tears on flexing 14 25 38 
Khaki Tears on flexing 14 24 34 
Olive Tears on flexing 9 13 19 
Dark blue Tears on flexing 8 12 18 
20 Dark blue Tears on flexing 14 20 26 
Rainwear (double-neoprene- 9 Olive Begins to melt 5 9 13 
coated nylon twill) 9 Olive Tears on flexing 8 14 22 
Drapery Fabrics ; 
Rayon gabardine 6 Black Ignites 9 20 26 
Rayon-acetale drapery 5 Wine ignites 9 22 28 
Rayon gabardine 7 Gold ignites =e 247 28+ 
Rayon twil) lining 3 Black lgnites 7 17 23 
Rayon twill fining z Biege Ignites 13 20 28 
Acetate-shantung 3 Black - Ignites 107 ait 357 
Cotton heavy draperies 13 Dark colors — Ignites 1S 18 34 
Tent’ Fabrics 
Canvas (cotton) 12 White Ignites 3 28 51 
Canvas 12 Olive drab lgnites 2 18 28 
Other Fabrics 
Cotton chenille bedspread Light blue Ignites abd iit 15+ 
Cotton venetian blind tape. White Ignites 10 18 22 
dirty ; 
Cotton venetian blind tape White Ignites 13% 27+ 31t 
Cotton muslin window shade 8 Green Ignites 7: 


13 19 


: i exposures for the indicated responses (except where marked +) are estimated to be valid to 225% under 
dard laboratory conditions. Under typical field conditions the values are estimated to be valid within 250% with a 
greater likelihood of higher rather than lower values. For materials marked +, ignition tevels are estimated to be valid 
within 450% under laboratory conditions and within 100% under field conditions. For low air bursts, values of ‘max of 


0.2. 1.0. and 3.2 sec correspond roughly to yields of 40 kt, 2 Mi, and 24 ML, respectively. 


: are not available ot appropriate scaling not known, 
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Table 9-3. , Approximate Radiant Exposures for Ignition of Various Materials 


Radiant Exposure’ 


" 
(cal:em-} 
Weich: ; fect 
; eight Eftec: ; nay tinay bnas 
Muteric! oes faz vd-) Color of Materia! O.2 sec 1.0 seo 3.2 see 
Howsenoid Tings: 
Newspaper. shredded = Ignites ' 4 o 1] 
Newspaper, dark picture ares 2 lentes $s pe ee 
Newspupe:. printed (ext ares = Ignites 6 & 1s 
i Greer: fenitcs B 6 9 th 
x Tan ignites 10. 1s 20 
iva Dar Ignites 6 art dae 40 
lo Brower, Ignites 16 20 40 
2 Jenites 24° 307 50% 
lentes a IS 20 
dgnites io] ey i 
Ignites OF 18+ “212 
Ignites 107 Jur Qe 
Ignites I= }4e 2 
bv ot perdetesa pine Pte ie Light vellaw Ignites 23° 23 
Tider Maternal”? 
mated woes pum. Git : I ar oF at 
jucuis deaves (beech: gaites 4 6 & 
Brass tchewt) Ignites 5 & 1G 
< Giass (sedge) lenites Oo 2 y “yy 
Pine needies. brown (ponderosy? fenites 1G 16 ol 
tenais 
. mineral surface lgnites ime >is Pile 
2 smooth surface Ignites =e 30 ta 
Piywoe gigs Fr : Fiaming during a Ib 20 
exposure 
Rubber. pal2 Jatex lenites Su & 110 
Rubber. black [ letiites 1G tu 25 


‘Other Materials 


Alumitain aircraft. skin (0.020 in, ‘ Blisters 15 30 40 
thiskhy cogied with 0.002 in. ot 
standard whpte aircraft puss: 


Cotton canvas sendbugs. dey filled Failure 10 Ie 32 
. Coral sand Explodes 18 a 47 
. (popcorning) 

Sihceous sand Explodes lt 9 35 

{popcorning} 


§:. exposures for the indicated responses (except where marked 7) are estimated to be valid to 4257 under 

ndard laboratory conditions. Under typical field conditions the values are estimated.to be valid within S407 with 3 
greater likehhood of higher rather than lower values. For materiats marked 7, ignition levels ase estimated to be valid 
within +50% under laboratory conditions and within 2100¢ under held conditions. For low air burst, values of tnan 
of 0.2, 1.0. and 3.2 sec correspond roughly to yields of 40 k1, 2 Mi. and 24 M1. respectivels 

= Data ate not available or appropriate scaling not known, 

or Kadsant exposures for ignition of these substances are highly dependent on the moisture content. 
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_ WB survive iw rire areas | | 


The best documented fire storm in his- 
tory (but not the one causing the greatest loss of 
life) occurred in Hamburg, Germany during the 
night of July 27-28, 1943, as a result of an 
_incendiary raid by Allied forces. Factors that 
contributed to the fire included the high fuel 
loading of the area and the large number of 
bujldings ignited within a short period of time. 
The main raid lasted about 30 minutes. 
Since the air raid warning and the first high ex- 
plosive bombs caused most people to seek shel- 
ter, few fires were extinguished during the at- 
tack. By the time the raid ended, roughly half 
the buildings in the 5 square-mile fire storm area 
were burning. many of .them intensely. The fire 
storm developed rapidly and reached its peak in 


two or three hours. 
a Many people were driven from their 
h 


shelters and then found that nearly everything 

was burning. Some people escaped through the 

streets; others died in the attempt; others return- 

ed to their shelters and succumbed to carbon 
oxide poisoning. 

Estimates of the number that were killed 
range trom about 40,000 to 55,000. Most of the 
deaths resulted from the fire storm. Two equally 
heavy raids on the same city (one occurred two 
nights earlier; the other, one night later) did not 
produce fire storms, and they resulted in death 
rates that have been estimated to be nearly an 
order of magnitude lower. 

More surprising than the number killed 
is the number of ‘survivors. The population of 
_the fire storm area was roughly 280,000. Esti- 
mates have been made that about 45,000 were 
rescued, 53,000 survived in non-basement shel- 
ters; and 140,000 either survived in basement 
shelters or escaped by their own initiative. 


9-25 Causes of Death || 


The evidence that can be reconstructed 
from such catastrophes as the Hamburg fire 
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storm indicates that carbon monoxide and ex- 
cessive heat are the most frequent causes of 
death in mass fires. Since the conditions that 
offer protection from these two hazards gener- 
ally provide protection from other hazards as 
well, the following discussion is limited to these 
two causes of death. 
Carbon Monoxide. Burning consists of a 
series of physical and chemical reactions. For 
most common fuels, one of the last of the reac- 
tions is the burning of carbon monoxide to form 
carbon dioxide near the tips of the flames. If the 
supply of air is limited, as it is likely to be if the 
fire is in a closed room or at the bottom of a pile 
of debris from a collapsed building, the carbon 
monoxide will not burn. completely. Fumes from 
the fire will contain a large amount of this taste- 
less, odorless, toxic gas. _ 

During the Hamburg fire, many base- 


ment shelters were exposed to fumes. Imperfect-  - 


ly fitting doors and cracks produced by explod- 
ing bombs allowed carbon monoxide to pene- 
trate these shelters. The natural positions of 
many of the bodies recovered after the raid indi- 
cated that death had often come without. warn- 
ing, as is frequently the case for carbon. mon- 
oxide poisoning. - 

Carbon monoxide kills by forming a 
more stable compound with hemoglobin than 
either oxygen or carbon dioxide will form. 
These latter are the two substances that hemo- 
globin ordinarily carries through the blood 
stream. Carbon monoxide that is absorbed by 
the blood reduces the oxygen carrying capacity 
of the blood, and the victim dies from oxygen 
deficiency. 

As a result of the manner that carbon 
monoxide acts, it can contribute to the death of 
a person who leaves a contaminated shelter to 
attempt escape through the streets of a burning 
city. A person recovering from a moderate case 
of carbon monoxide poisoning may feel well 
while he is resting, but his blood may be unable 


to supply the oxygen his body needs when he 
exerts himself. After the air raid at Hamburg. 
victims of carbon monoxide’ poisoning. appar- 
ently in good health. collapsed and died from 
the strain of walking away from a shelter. It is 
suspected that many of the people who died in 
the streets of Hamburg were suffering from 
incipient carbon monoxide poisoning. 

Heat, The body cools itself by perspira- 
tion. When the environment is so hot that this 
_ method fails. body temperature rises. Shortly 
thereafter, the rate of perspiration decreases 
rapidly. and. unless the victim finds immediate 
relief from the heat, he dies of heat exhaustion. 
Death from excessive heat may occur in an in- 
adequately insulated shelter: it also may occur in 
the streets if a safe area cannot be located in a 
short time. 


9.26 Shelters QB 


The results of the Hamburg fire storm. 


illustrate the value of shelters during an intense 
mass fire. The public air raid shelters in Ham- 
burg had very heavy walls to resist large bombs. 
Reinforced concrete three feet thick represented 
typical walls. Some of these shelters were. fitted 
with gas proof doors to provide protection from 
poison gas. These two features offered good pro- 
tection from the heat and toxic gases generated 
by the fire storm. 

The public shelters were of three tpes: 
Bunkers. These were large buildings of 
several shapes and sizes, designed to with- 
stand direct hits by large bombs. The fire 
storm area included 19 bunkers designed to 
hold a total of about 15,000 people. Prob- 
ably twice this number occupied the 
bunkers during the fire storm, and all of 
these people survived. 


® Splinterproof Shelters. These were long 


single story shelters standing free of other - 


buildings and protected by walls of rein- 
forced concrete at least 2-1/2 feet thick. 


No deaths. resulting from the fire storm 
were reported among occupants of these 
shelters. These structures were not gas- 
_ proof. Distance from burning structures 
and low height of the shelters probably 
provided protection: from carbon mon- 
oxide. 

@ Basement Shelters. The public shelters that 
were constructed in large basements had 
ceilings of reinforced concrete 2 to 5 feet 
thick. Although reports indicate that some. 
of the occupants of these shelters survived 
and some did not, statistics to indicate the 
‘chance of survival in such structures are not 
available. 

@ Private Basement Shelters, Private — base- 
ments were constructed solidly, but most 
of them lacked the insulating value of verv 
thick walls and the protection of gas-tight 
construction. Emergency exits. (usually 
leading to another shelter in an adjacent — 
building) could be broken if collapse of the 
building caused the normal exit to be 
blocked. As a result of the total destruction 
in the fire storm area, this precaution was 
of limited value. Many deaths occurred in 
these shelters as a result of carbon mon- 
oxide poisoning, and the condition of the 
bodies indicated that intolerable heat fol- 
lowed the carbon monoxide frequently. In 
some cases, the heat preceded the poison- 
ous gas and was the cause of death. Gen- 
erally, these shelters offered such a small 
amount of protection that the occupants 
were forced out within 10 to 30 minutes. 
Most of these people were able to move 
through the streets and escape. Others were 
forced out later when the fire storm was 
nearer its peak intensity, and few of these 
escaped. A few people survived in private | 
basement shelters. 

Experience in Hamburg and other mass 


fire areas suggests the following requirements for 
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shelters for fire protection: 


® Location. Shelters should be located as far 
away from combustible structures as pos- 
sible. The bottom of a mass of burning 
debris is very hot and can remain so for 
days. An open area also may be exposed to 
considerable heat, but not so much as the 
basement of a burned-out building. Three 
feet of earth will provide sufficient insula- 

' tion for a shelter if no heavy fuels are near- 
by. However, material scattered by the 
blast wave could fall on a shelter and 
render it less safe. The choice of a suitable 
location for a shelter is most difficult in 
heavily built-up areas, where mass fires are 
most likely to occur. Shelters in such areas 
would have to be designed to withstand the 
severe environment to which they would be 
subjected in case of a mass fire. When the 
shelters in Hamburg were built, the prob- 
lems of mass fires were not anticipated. 

@ Ventilation. Any shelter that is large 
enough to house its occupants in reason- 
able comfort contains sufficient air to sus- 
tain life for the duration of a mass fire. The 
shelter should be constructed so that it can 
be sealed from the entry of gases from the 
outside during the period of active burning. 
After the fire subsides, air will be available 
inside the shelter if burning debris has not 
fallen on or near the air intake. 

@ Provisions. \f escape from a shelter is nec- 
essary, wet coats or blankets would im: 
prove the chance of survival in the open. 
Since any shelter built for protection 
against a mass fire would logically serve 
also as a fallout shelter, such items as 
blankets and water normally would be 
available. 


9-27 Escape from the Fire Areas @ 


A large number of people, chiefly the 
occupants of basement shelters, escaped from 
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the Hamburg fire simply by leaving while the 
streets were still passable. In areas where damage 
is sufficiently light that people can attempt 
escape, the rate at which the fire builds up is 
expected to be slow enough to allow 30 minutes 
or more before movement through the streets 


will become dangerous. 
Escape from the Hamburg fire storm 


area was. simplified by the limited size of the 
fire. In the event of a mass fire caused by a 


Nuclear attack, the fire area probably will. be 


much larger. Although escape to the edge of the 
fire may be impossible, parks, bodies of water, 
or even areas that are not heavily built up may 
offer relative safety from the fire. 

As the fire grows in intensity, selection 
of'a suitable escape route becomes important. 
Moving through a narrow street, with tall build- 
ings burning intensely on both sides, will subject 
people to an excessive amount of heat. 

Preplanned escape routes are desirable. 
The safest routes are those that minimize ex- 


- posure to thermal radiation from burning struc- 


tures. A safe street would be wide compared to 
the heights of the buildings facing it. Masonry 
buildings with few windows would be a lesser 
threat than buildings with many windows. The 
solid wall blocks radiant exposure from fires 
burning inside the building, and-it also protects 
the contents of the building from direct ignition 
by the thermal pulse from the nuclear weapon. 

Simple rules for distinguishing safe from 
unsafe streets are not available; however, calcula- 
tions of the street width for which the thermal 
radiation would be intense enough to ignite 
clothing in a short time have been made, assum- 
ing that buildings on both sides of the street 
were burning. The results of these calculations 
are shown in Table 9-4. 

Whether or not such streets are safe will 
depend on many factors, such as fire intensity, 
percentage of the building fronts occupied by 
windows, wind speed and direction, protective 


eee 


clothing. and exposure time. In general, the 
Streets of Hamburg were 45 to 60 feet wide. and 
the buildings were 3 to 5 stories high. Table 
9-4 predicts that such streets would be danger- 
ous. and experience shows that they were. How- 
ever. even these streets provided escape for 


some. 


Table 9-4, oe Minimum Structure Separation 
for Escape Route 


Building Height 


(feet) (feet) 
20 40 
30 57 
40 67 
SQ &3 
60 — 96 
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Experience indicates that large open 
areas within a fire storm area probably are safe. 
A park 1.000 feet in diameter provided adequate 
protection from the Hamburg fire storm: a simi- 
lar area of 400 x 400 feet did not. 

@ In Tokyo. during the Kanto earthquake 
and fire of 1923. fire whirlwinds sweeping across 
some large open areas killed many who would 
otherwise have survived. Thus. safety in a partic- 
ular location depends to some extent on unpre- 
dictable aspects of the fire. 


SECTION IV 


hace THERMAL RADIATION 
DEGRADATION OF STRUCTURAL 


RESISTANCE TO AIR BLAST 


Section II] of this chapter contains a 


description of the properties of materials that 
might result in ignition or degradation of their 
physical properties. The emphasis of the discus- 


Distance Between Buildings” 


sion in Section III concerns the ignition of com- 
bustible materials. This section provides a some- 
what expanded treatment of the degradation of 
structural resistance to air blast that is brought 
about by exposure to thermal radiation. | 

The problem of integrated thermal/blast - 
effects can be divided into several elements. One 
element is the free field thermal and air blast 
environments to which a system element is 
exposed. Methods of predicting these environ- 
ments are presented. in Chapters 2 and 3. As 
pointed out in Chapter 3, prediction of the ther- 
mal environment is difficult because of the 
importance of climatic conditions. A. statistical 


or probabilistic description of weather condi-. ~ 


tions frequently will have to be used to form 
some estimate of the validity of the results. 

The first element of the overall problem 
treated in this section is the coupling of thermal 
energy into the structure or system element of 


~ interest. The condition of the target surface. its 


orientation to the detonation. and coatings 
employed all can affect the amount of thermal 
energy that is absorbed by the target. © 

The. second eleiment is the mechanisms 
for the loss of energy from the structure through 
convective heat losses or reradiation.: The third 
element is the effect of absorbed energy on the 
state of the target material. The primary effect 
of this change of state is the fourth element. 
which is the change in material properties. 

The final and fifth element of the prob- 
Jem is the effect of changes in material 
properties on the structural resistance to air 
blast loading. 

a Insufficient data are available concerning 
the thermophysical properties of nonmetallic 
materials to provide a realistic discussion of the 
combined blast/thermal effects on such ma- 
terials. Any current prediction would be highly - 
uncertain, both with regard to the occurrence of 
any specific physical phenomena and to the 
quantitative evaluation of the phenomena. Con- 
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sequently the discussion in this section will pe 
limited to metallic materials. 


GBR thermar enercy assorsen §& 


9-29 Absorption a 


The amount of thermal radiation that a 
metallic element will absorb depends on the 
properties of the metal, color, surface condition, 
orientation to the burst, the absence or presence 
of some type of coating (either accidental such 
as dirt or grease or intentional such as a paint 
system).* The amount of thermal energy ab- 
sorbed by a target usually is determined by 
multiplying the incident energy by an absorp- 
tion coefficient. The absorption coefficient can- 
not be directly measured; however, a good esti- 
mate can be made by performing spectral mea- 
surements of reflectance over the range of wave- 
lengths representative of the nuclear radiation 
spectral distribution, and finding an average 
value of reflectance. The coefficient of absorp- 
tance then equals one minus this average value 
of reflectance. The absorptance determined in 


this manner is valid only as long as the surface - 


does not change as a result of the absorbed ther- 
mal energy. Another approach is to divide the 
quantity of energy absorbed by the incident 
energy. The absorbed energy is determined from 
experimentally determined temperature vs time 
relationships and appropriate theoretical energy 
balance relationships.’ Some typical values for 
absorption coefficients for bare metals are 
shown in Table 9-5. A value of 0.5 would be a 


good value to assume for the absorption coef- 


ficient in the absence of more definitive data for 
a specific system under investigation. 

A coating or surface treatment on a 
“metal substrate can alter the effective absorption 
coefficient significantly, and thereby can alter 
the amount of energy absorbed by the metal 
structure. The amount of energy absorbed 
depends on color, thickness, adhesiveness, and 
heat transfer characteristics of the coating. The 
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Table 9-5. 8 Absorption Coefficients for 
are Metals . 


Absorption 
Material Coefficient 
Polished metals 0.25-0.5 
Especially clean aluminum 0.2-0.4 
Clean aluminum 0.5 
Unpolished metals 0.45-0.55 
Average over. aircraft 0.5 


most common type of coating is, of course, 
paint. The importance of paint in altering the 
amount of thermal radiation absorbéd was noted 
in early investigations where it was found that 
the thin aircraft skin under painted insignia had 
melted, whereas the skin of adjacent areas was 


not affected. lt should not be assumed, however, 
that paint or other coating is necessarily a detri- 


ment. In some instances.the smoké that results 
from the breakdown. of the paint can reduce the 
amount of energy absorbed and thereby reduce 
the maximum temperature reached by the sub- 
strate. Figure 9-8 illustrates this factor. The only 
true genera] statement that can be made about 
coatings is that they can alter the amount of 
thermal energy absorbed by the substrate. In 
fact some coatings, even paint systems, can be 
designed to reduce the amount of energy absorb- 
ed. As a general rule, however, most coating 
systems on present tactical military equipment 


‘G Implicit in this statement and throughout this discussion is 
the fact that energy coupling depends on the wavelength of the 


incident radiation. 


Absorptance is a function of time and temperature, al- 
though an average value usually is employed. This function when 
known may be found as an input specification for certain com- 
puter programs. 
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Histories of Aluminum Honeycomb. 
Gray Painted Skin: Core 1/8-in., Cell Size 1/2-in. 
Thick. 21-24.2 cal/em?, 3.8 Mt 


will cause an increase in the amount of energy 
absorbed. The mechanisms of the response of 
coatings and the transfer of energy to the sub- 
Strate are so complex that theoretical descrip- 
tions have not been successful. The response of 
coatings depends upon both the rate at which 
thermal radiation is received and the total inci- 
dent energy. A given amount of radiation re- 
ceived at a slow rate may only cause discolora- 
tion or slight scorching. whereas if it were 
received at a very rapid rate it would cause 
charring or blistering. The latter response is 
more typical in a nuclear environment. The 
change in state of the coating reduces the 
amount of energy available to be transferred to 
the substrate. because these processes:‘are irre- 
versible. However, the change in state usually 
increases the absorption coefficient of the coat- 
ing. and the net result usually is an increase on 
the amount of energy absorbed by the substrate. 
An exception occurs if the absorption coef- 
ficient was initially high, and energy was 
absorbed by the smoke layer generated by the 
decomposing paint. Coatings that respond by 


Facing, 0.016-in. 


physically separating from the surface, Le. 
blister. potentially can reduce the amounts of 
energy transferred-to the substrate. Since the re- 


-sponse of the coating-substrate system to nu- 


clear weapon thermal radiation is so complex, 
effective absorption coefficients: usually are 
determined experimentally. 

The average value of the absorption co- 
efficient. A, depends on the properties of the 
substrate as the flow of thermal energy from the 
coating influences the response of the coating. A 
series of experiments with laser heating provides 
data to illustrate this point. Various metallic 
samples were prepared with a Silicon-carbide- 
synar binder coating (spectrophotometric mea- 
surements gave an average absorption coefficient 
of 0.85 for the coated aluminum sample. for 
equilibrium temperatures of 70 and 600°F fora 
10.6 micron wavelength). The average values of 
the absorption coefficients for a range of irrad- 
iances based on temperature vs time data are 
shown in Table 9-6. The effect of substrate 
properties is significant. 

Values of absorption coefficients are 
shown in Table 9-7, Generally, values range from 
0.3 to 0.6 for light or reflective paints, and from 


Table 9-6. Average Absorption Coefficients 
for Laser Heated Coated Substrates 


Average 

Absorption 

Substrate Coefficient 

2024-T81 clad aluminum alloy 0.7 + 0.1 
6AL-4V annealed titanium 0.35 + 0.1* 
301 full hard stainless steel 0.80 + 0.05 


FB ccarssion of titanium specimens showed con- 
siderable scatter. Erratic behavior probably resulted 
from improper bonding of the coating and cracking or 
flaking of the coating during straining. 


Table 7, Be Representative Values of Absorption Coefficients for 
Metals with Various Coatings or Surface Treatments 


Coating or Absorption 
Substrate Surface Treatment Coefficient 
Metallic skins Aluminum paint 0.5 
Painted metals - White paint 0.3 - 0.5 
(aircraft skins) White paint (clean) 0.25 
- White paint w/oil film 0.30 
Yellow paint. 0.4 --0.55 
Olive paint 0.6 - 0.7 
Gray paint (clean) 0.6 
Black paint 0.65 - 0.95 
Insignia blue : 0.9 
5032-H32AL Black paint OF: = 40 
Aluminums Reflecting white paint | 0.32 
. Camouflage paint 1.00 
6061 Al-Mg alloy Haze gray Navy paint 0.78 + 0.04 
Volcano] , 0.99 + 0.05 
2024-T3 Anodized black 0.67 


0.6. to 1.0 for dark paints or treated surfaces. 
e It is apparent that the selection of an 

absorption coefficient based on the information 
in Table 9-7 for use in analyzing the response of 

_a military system could be in error by as much 
as +50 percent. Currently, the best solution is to 
perform upper and lower bound calculations. 
The upper and lower bound calculations should 
be based on the best estimate of the condition 
of the surface at the time of exposure, i.e., dirty, 
corroded, prior exposures, etc..— 

One additional factor affects the amount 
of energy absorbed. This factor is the orienta- 
tion of the surface to the thermal pulse. Ideally, 
if incident parallel rays are assumed, with no 
scattering, reflection or refraction, the target 
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surface receives all of the radiation if it is 
oriented normal to the rays. If the. surface is 
rotated about an axis normal to the rays, it cuts 
fewer of the rays, and the amount of thermal 
energy received is reduced. The energy received 
by a differential area on a curved surface also 
can be treated in this manner. With such assump- 
tions, the absorbed energy function takes the 
following form: 


Q, = QA, cos 8 
where 


Q = incident energy, usually expressed in 
cal/cm? 


Q, = absorbed energy, usually expressed in 
cal/cm? 


nae 


-l. = effective absorption coefficient 


tod 


& = incident angle. angle between a line from 
the source.to the surface and a line nor- . 
mal to the plane of the surface. 

Under those ideal conditions. the amount of 


absorbed energy follows the cosine law; how- 
ever. there will be conditions where scattering. 
reflection. and refraction, will invalidate the 
parallel ray assumption in addition to altering @. 
This would principally occur at large angles of 
incidence such as shown schematically in Figure 
9-9, A computer program called TRAP has re- 
cently been developed to handle - these 
conditions. 

Figure 9-9 also illustrates the amount of 
energy received by the surface of a cylinder with 
its axis oriented normal to the ray paths. The 
quantity @ then becomes the distanee. in ra- 


15 


Q0/Q 


Figure 9-9, 


w/4 


~ 


dians. along the circumference. from the point 
nearest the source. 

Both the absorption coefficient and 
orientation are important in determining the 
amount of energy absorbed into the target. 
Errors in orientation are of greater significance 
at large angles on incidence; however, orienta- 
tion usually can be inferred, at least for worst 
cases, whereas the absorption coefficient usually 
is not well known unless it has been determined 
experimentally. 


9-30 Energy Losses =z 


The absorption of. thermal radiation 
energy into a material has been discussed in the 
previous paragraph. Before discussing the effect 
of this absorbed energy on material properties or 
change in structural resistance. the form and sig- 
nificance of energy losses from the material 


WITH SCATTERING, ETC. 


T/2 
@ 


Schematic Comparison of Ideal Cosine Law 


Refiection, and Refraction 


with Function Inctuding | 
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must be discussed. The three ways in which ther- 
mal energy can be dissipated are conduction, 
radiation, and convection. For purposes of this 


discussion, conductive dissipation is not really a — 


loss, in that energy is redistributed throughout 
the material system rather than lost from.the 
system. This form of heat dissipation is consid- 
ered in more detail in the following subsection 
that discusses the change in material state, and 
the discussion of temperature profiles. 

a Radiation losses normally are considered 
négligible, because the temperature rises in a ma- 
terial are relatively small, even though the radia- 
tion heat loss is a function of the fourth power 
of absolute temperatures. The heat loss through 
radiation is expressed by: 


Q. = radiative heat loss, cal/cm?-sec, 


o = Stephan Bolzman Coefficient = 1.36 x 


. re . 
10°12 cal/cm? sec °K™*, 
€ = emissivity, nondimensional, 


T = absolute temperature. °K, of radiator 
(may be time dependent), 


T) = absolute temperature, °K, of surround- 
ing air. 

If a 200°C temperature difference is assumed, 

and an emissivity of 0.2 (for oxidized aluminum 

at temperatures of a few hundred °C), and the 


ambient air temperature is 300°K, the radiative. 


heat loss becomes: 


Q, = (1.36 x 107!*)(0.2)(5004 - 300*) 


fl 


1.48 x 1072 cal/cm? «+ sec 


If the temperature rise is 400°C then Q. be- 
comes equal to 6.3 x 10°? cal/cem? sec. Thus, it 
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can be seen that even at high elevated tempera- 
tures, the heat loss through radiative processes is 


ll 
es Convective heat losses occur when there 
is a flow of air over the heated surface. This type 
of heat loss is most important for aircraft in 
flight and is a weak function of position on the 
airfoil, and the speed and altitude of the aircraft. 
The rate of convective heat loss also depends on 
whether the boundary layer flow is laminar or 
turbulent. These processes are generally under- 
stood, and a number of formulas that are paver 
able are given below: 


: V 0.8 
= 0.67(10)4 793 SS 
0.2 
x 
where 
h. = convective heat loss. coefficient for tur- 


bulent flow, 
T = air temperature, °R, | 
V = flight speed, ft/sec, _ 
p = density of air, Ib/ft?, 
x 


= distance from leading edge of airfoil, ft. 


ity = 0.176 pVC,(log,, Rey™S, 

hey = 0.01445 k/x (Re)®. Pr), 

h,, = 0.0282 k/x (Re)? (Bry 

hi = 0.318 k/x (Pr)¥3 (Re)'?, 

hi, = 0.332 pVC, (Pry 73 (Rey ¥?, 
where 


h cy 7 convective heat loss coefficient. for tur- 
bulent flow, cal/cm? - sec + °C 


h’. = convective heat loss coefficient for 
laminar flow, calfem* + sec + °C 
p = density of air, gm/cm? 
V = flight speed, cm/sec 


C. = specific heat of air. cal/gm - °C 


k = thermal conductivity of air, cal/cm + sec 
. °C 

x = distance from leading edge of airfoil, cm 

Pr = Prandtl number. Caulk, dimensionless 

Me = viscosity of air, gm/sec * cm : 
These formulas are based on the assumption that 
_a reasonable approximation may be obtained by 
analysis of airflow over a heated flat plate. 
Experimental. data indicate that the average 
error in using this approximation may be as 
much as 30 percent with a mean error of from 


10 to 15 percent. The value for/,, ranges from 
zero to 0.03 cal/em* + sec * °C for a range of 


reasonable air properties and velocities. General- ~ 


ly, convective heat losses are marginally impor- 
tant for relatively low yield weapons; however, 
for long duration thermal pulses and under con- 
ditions when the temperature difference 
between the aircraft skin and the air boundary 
layer is large. the convective heat loss can be 
significant for high speed aircraft. 

Values of 1, have been calculated (see 


Table 9-8) for purposes of comparing the results © 


of the.equations given above, under the follow- 
ing assumed conditions. 


Aircraft type = fixed wing 

Aircraft speed = 100 mph (4.47 x 10° cm/sec) 
Aircraft altitude = 1,000 ft | 
Air temperature = 70°F, 20°C 


4 ft, 121.92 cm 


x 


p = 7.5(10)? lb/ft? 


Re = Reynolds number, pl’x/p, dimensionless 


Turbulent 


Table 9-8: Values of h. Obtained from 
Various Functions in the Order in Which the 
Equations are Given Above @ 


Flow 
Condition 


2.2710)? 
2.2710)>— 


Laminar 2.15(10)4 


2.82¢10)"4 


1.2010)? gm/em3 


k = 6.0874(10)° cal/em * sec * °C 
C, = 0.2404 éal/gm + °C 
uw = 1.798(10)4 pgm/em + cm 


The use of the convective heat loss coef- 
ficient, #1, is shown below: 


Cee = hey (7, ~f£) 


Q., = convective heat loss, cal/em? = sec 


T_ = surface temperature, °C (may be time 
dependent) 


Ny 
tt 


surrounding air temperature, °C 


For high speed aircraft, it may not be appro- | 
priate to assume that the temperature of the sur- 
rounding air is the temperature of the air 
boundary layer next to the heated surface. This 
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is particularly true if aerodynamic heating 
occurs. For this situation, the recovery tempera- 
ture 7, should replace 7 in the equation given 
above. 

bcd For tactical systems of interest to this 
section, radiative and convective heat losses gen- 
erally can be neglected. However, these are and 
normally should be included in detailed analyses 
and computer calculations for completeness, 
particularly for high speed aircraft and large 
yield weapons. 


CHANGES IN MATERIAL STATE 
AND MATERIALS PROPERTIES 
9-31 Changes in Material State 
The change in material state of primary 
interest is the change in temperature caused by 


absorbed thermal radiation. The temperature _ 


changes are controlled by the energy absorbed in 
the material as a function of time and position 


and by the geometry of the member being ana- 


lyzed. Geometry can have a significant impact 
on heat flow through. the member as well as heat 
absorption and losses. There are a number of 
techniques available for determining tempera- 
ture changes. These techniques vary from simple 
formulas to large computer programs depending 
on the simplifying assumptions that are made 
and the complexity of the system geometry. The 
greatest amount of effort has been devoted to 
the study of temperature changes in flat plates, 
primarily because of the interest in aircraft 
safety/survivability problems. Therefore, the fol- 
lowing discussion is dominated by analysis of 
flat plates, although most of what can be learned 
about the processes affecting plate temperature 


applies to other geometries because the basic 


equations describing the processes are the same. 

Plates can be divided into several types 
according to the temperature profile through the 
plate. A thermally “thin” plate usually is defin- 
ed as.a plate where the temperature gradient 
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between the front and back of the plate is small. 


_ Thermally thin plates simplify the analysis prob- 


lem, because, if a significant thermal gradient 
exists, the thermal stresses become more com- 
plex. Plates usually are assumed to be one 
dimensional, with no heat flow to a supporting 
structure, so heat sink effects can be neglected, 
unless the plate is part of a fuel tank. 

Whether a plate is thermally thin de- 
pends not only on the thermophysical properties 
of the material, but also on the thermal pulse 
characteristics and the time during the pulse that 


. the criterion is applied. Some materials conduct 


heat more quickly than others, and the tempera- 
ture on the back face of these materials will not 
lag behind the temperature on the front as much 
as for materials with lower conductivities. For 
example, the thermal conductivity of aluminum 
is about 0.4 cal/sec + cm * °C, and that of 
titanium is about 0.04 cal/sec + cm - °C. There- 
fore, aluminum conducts. 10 times as much heat 
as titanium in a given amount.of time. Other 
material factors that control the ratio of front to 
back temperature are specific heat, density, and 
thickness. 
The temperature difference also depends 
ulse characteristics. For a pulse that is 
delivered slowly, there may be sufficient time 
for energy to be conducted to the back face, so 
the thermal gradient at the time of interest is 
minimized. Conversely, a pulse that is delivered 


_ Tapidly may cause a steep thermal gradient. For 


a given pulse, material, and thickness, the time 
at which the temperature response is of interest 
determines whether the thermal gradient is 
small. 
= A number of techniques have been 
eveloped to determine whether a plate can be 
considered to be thermally thin. One technique 
that gives a good approximation is illustrated in 
Figure 9-10. The values of n = ¢/t,,,, and 
7 aeibe must be evaluated to use the curves in 
Figure 9-10. From Chapter 3 


= 0.043 14? (p/p)? sec,* 


fT 
Man 


where 7. is the time to the principle thermal 
maximum. I’ is: the weapon yield in kilotons, 
and pp. is the ratio of the air density at the 
burst altitude to the air density at sea level. The 
parameter @ is defined as 


’ 4) 
a= kip C. cm/sec 
mp 


whiere 


k = thermal conductivity of air. cal/em + sec 
co] 
f° G.. 


material density. gmécm>. 


Cy 
and & js the thickness of the plate in centi- 
meters. The point where the values of 7 and 


OPS 6 ‘h? meet determines whether the plate is 


Ut 


specific heat of air. cal’gm * °C 


thin, finite. or thick. A thin plate is one in which 

the difference between back and front tempera- 

tures is less than 10 percent. A thick plate is one 

in which there is no temperature rise on the 
back face. A finite plate falls between the other 

two criteria. 


ae A technique for predicting the tempera- 
ures of thermally thin plates uses the relation- 


ships shown in Figure 9-11. This technique is 
believed to provide good results over a fairly 
wide range of problems of interest. The tech- 
nique assumes that thermophysical properties 
and the rate of convective heat loss are constant. 
that radiation losses are negligible. and that 
there is no heat loss from the back of the plate. . 


B®: €xamples provided in this section will al] be concerned 
wien air bursts that are sufficiently low that the term (pip, 4+ 
is near enough to 1.0 to be neglected (e.g.. at 4,000 feet 


(0'09)°-4? = 0.99). Therefore, the approximation t,,4. 0.043 » 
0.43 will be used throughout the remainder.of this section. 
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Problem 9-1, Calculation of Thermal Thickness of a Metal Plate 


fa Figure 9-10 contains curves that define 
regions where metal plates may be considered 
thermally thin, finite, and thick, respectively. 
The curves are plotted as a function of the 
parameters 7 and ar, ,,/b?, which are defined 
in paragraph 9-31. Their use is demonstrated in 
the following ex : 


ample. 
Exam il 
Given: A 5086 Aluminum alloy plate is an 


important part of a structure under analysis. The 
plate is 0.875 cm thick and the range of air blast 
exposures for the system of which it is a part are 


expected to be from 8 to 15 psi from a 100 kt | 


explosion at a height of burst of 1,000 feet. 

Find. Whether the plate can be considered 
thermally thin. 

Solution: The corresponding height of 

burst for a 1 kt explosion is 


dw'3) (oo)!/3 

From Figures 2-18 and 2-19, Chapter 2, the 
ground distances from a |! kt explosion at a 
height of burst of 215 feet that correspond to 
-15 and 8 psi overpressure are 895 and 1,250 feet 
respectively. From Figure. 2-28, Chapter 2, the 
times of arrival of the blast wave from a 1 kt 
explosion to these distances are 0.37 seconds 
and 0.61 seconds, respectively. The correspond- 
ing times for a 100 kt explosion are . 


_ 1/3 
t= tw 
t = (0.37)(100)"3 = 1.72 sec for 15 psi, 
and 


t = (0.61)(100)/3 = 2.83 sec for 8 psi. 


The time to final maximum for a 100 kt low air 
burst is 
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t_.. & 0.043 yo43 


max ; 
= (0.043)(100)%4? = 0.31 sec. 


Therefore, 


n= tte ax 
1.72 
"= 931 ~ 5.5 for 15 psi, 
and 
83 
n= 031° 9.1 for 8 psi 


The properties of the alloy are 


k = 0.28 cal/sec + °C, 


C, = 0.22 cal/gm + °C, 
p = 2.66 gram/cm?. 
Since 
i ee 
pc, * 
Ol max - Sane 
are 2 
8 pC,b 
(0.28)(0.31) 
(2.66)(0.22)(0.875)? 
= 0.19. 


Answer: From Figure 9-10, at the 15 psi - 


overpressure level where n = 5.5, the plate 
should be considered finite; at the 8 psi over- 
pressure level, where n = 9.1, the plate can be 
considered thermally thin. 


wa Related Material: See paragraph 9-31. 


* 
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Figure 9-10. a Thin Plate, Thick Plate, and Finite Plate Regions for 
Exposure to Nuclear Weapons Thermal! Radiation 
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a The variables used in this technique have 
units in the cgs system and have been defined 
previously; however, they are listed below for 
convenience. 


T 


. temperature of the: plate prior to ex- 


posure, 


AT = temperature rise, °C, 


A = absorption coefficient 
Q = thermal exposure, cal/em? (effective 
9-42 


thermal exposure may be used where Q 


=Q cos 9), 
p,, = Specific gravity of material, gm/cm?, 
C= specific heat, cal/gm + °C, 
b = plate thickness, cm, 
h., = convective heat loss, cal/em? + sec - °C, 
n= tlt ax- 


Use of the technique is illustrated in Problem 
9-2.- 


eRe em me 


’ 


Problem 9-2, Calculation of Temperature Rise in a Thermally Thin Plate 


Figure 9-11 contains a family of curves 
that relate thermophysical properties of a thin 


metal plate to thermal pulse parameters in a ~ 


manner that allows calculation of the tempera- 
ture rise in the plate. The procedures for the 
calculation are made clear in the following 
examples. Symbols for the various parameters 
are listed in paragraph 9-31. 
Example 1 
Given: A 6061-16 Aluminum alloy plate is 
exposed to the thermal pulse from a 100 kt low: 
air burst at a location where the incident radiant 
energy is 76 cal’cm* and arrives at an angle of 
incidence of 30°. Properties of the plate are 


T, = 30°C, 

A= 08. 

po = at gm/cm?, 

Cus 0.216 cal’gm + °C, 
b = 0.15875 cm. 


h. = 0. 


cv 
Find: The temperature of the plate when n 


= 2,4, and 6. ' 
Solution: The effective radiant exposure is 


Q. = Q cos @ = 76 cos 30° = 65.8 cal/em?, 
From Figure 9-11, with 


ht 


cv max 


= 0, 
pcb 


when 


when . 
y= 4; ADC 0.67, 
and when 
is 6, Wie = 0.74. 
The value of AQ/pC,b is 
AQ. _ (0.8)(65.8) 


=. 569°C. 


pC,b ~ (2.7)(0.216)(0.15875) 


The temperature rises are 


AT = (0.47)(569) =.267°C for n = 2. 

AT = (0.67)(569) = 381°C for n = 4, 
and 

AT = (0.74)(369) = 421°C for n = 6. 


Answer: The plate temperatures are 


T= T. 4-AT, 
Therefore, 
Treg = 30 + 267 = 297°C 
T,=4 30 + 38) = 411°C 
Ty. = 30+ 421 = 451°C, 
Example 2 


’ Given: The same conditions as Example 1, 
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| ] 
~ . except thath., =3.0x 10°? cal/em? + °C.* The temperature rises are 
Find: The temperature of the plate when 7 : 
= 2,4, and 6. AT = (0.44)(569) = 250°C for n = 2, 
Solution: as 
. AT = (0.53)(569) = 302°C for n =-4, 
= 043 = 0.43 
Be 00S = (0,043)(100) ad 
= 0.31 sec. . 
oS AT = (0.51)(569) = 290°C for n = 6. 
evlmax © _ (3.0 x 1070031) _ g, 
pcb (2.7)(0.216)(0.15875) ne Answer: The plate temperatures are 
From Figure 9-11, when T= 7, + 7. 
. Therefore, 
es Tye, = 30 + 250 = 280°C, 
when . = 
T,-, = 30 + 302 = 332°C, 
pe ea oe ee nee 
AQ/pC b 
p 
and - 
and when ) 
a = 2 = - ages 
AT. 0.51. The 6 30 + 290 = 320°C. 


Fron Euanpiet eS Related Material: See paragraph 9-31. 


_AQ = 569°C This value corresponds roughly to an airplane speed of 
pC_b , q mi/hr. 


9-44 


TO ee ee eee + 


STE ea 


i } 

' 

\\ 

‘ [ 


So ro) 
a?9¢/0V 
iV 


Figure 9-11. @ Thermal Response of a Thin Plate to Nuclear Weapons Thermal Radiation Pulse & 
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ae Thermally thick plates are those that do 
not have a temperature rise on the back surface 
at the time of interest. Finite plates have a tem- 
perature gradient between the front and back 


s 


Computer codes generally are required 
to calculate the temperature history of ther- 
mally thick and thermally finite plates. Simple 
graphical techniques have been developed to 
determine the temperature rise on the front face 
of a thermally thick plate and on the front face, 
mid place, and back face of thermally finite 
plates. These are not included herein because the 
techniques depend upon thermal source charac- 
teristics that are very different from the source 
‘characteristics described in Chapter 3.* If graph- 
ical techniques that are compatible with the ther- 
mal source data of Chapter 3 become available, 
they will be included in revisions to this manual. 
In the meantime, users who are interested in the 
techniques should consult “Thermal Degrada- 
tion of Structural Resistance to Air Blast’’ (see 
bibliography). 


~ and a temperature rise at both surfaces. 


Requirements to analyze geometries that . 


are more complex than single layer plates (even 
if they are layered thin plates) or to examine 
plates with temperature dependent thermo- 
physical properties or heat losses require the use 
of a computer code. Codes for more complex 
geometries, such as cylinders and tee beams, 
have been developed specifically for the types of 
analysis of concern to this section. As noted 
previously, the TRAP computer program can 
calculate temperatures and stresses for two 
dimensional structures. 

Honeycomb structures or sandwich con- 
struction must be considered as a multi-layered 
system. A five layer system was used to examine 
thin sandwich construction. ‘The five layers con- 
sisted of the core, the two skins, and the two 
bonding layers between the skins and core. Most 
sandwich construction has been found to fail at 
the bonding layer. Consequently, the tempera- 
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ture of the bond and the thermal stresses in the 
face sheets (or skins) can be of equal impor- 
tance. Studies and experiments show that a cal- 
culational. procedure that includes five layers 
was required to duplicate experimental results 
accurately. A code developed for dual layer 
plates has been used to calculate thick single 
layer plates by considering both layers to be the 
same material. 


9.32 Change in Materials Properties a 


The rise in temperature that occurs 
when thermal radiation is absorbed can be 
expected to cause changes in materials prop- 


The important differences in source characteristics are the 
time to final maximum, where ; 


“ 1/2 
imax = 9.032 W!!? sec 


was used in developing the graphical techniques, while a more 


recent value 
fe 0.43), \0.42 
Tce 0.043 W (P/P,) sec 
is given in Chapter 3, and the power at thermal maximum, where 


= 1/2 
F na =4W kt/sec 


was used in developing the graphical techniques while a more 
recent value of 


is given in Chapter 3, The latter determines the maximum irradi- 
ance (thermal energy per unit area per unit time) incident on the 
target. Actually, the procedures described in Problems 9-1 and 
9-2 were developed for use with the older expression for. t.453 
however, since the two expressions agree within a few percent 
for low air bursts and for yields from a few kt to a few Mt, the 
use of the techniques with the value.of f,,,, from Chapter 3 was 


_ illustrated. The values of P,,,,. on the other hand, differ widely 


at all yields, Since the temperature. histories presented in the 
graphical techniques for finite and thick plates depend strongly 
on ittadiance, it is not possible to use the thermal source data 
with the curves derived for use with the older expression for 


Peas 
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erties. Properties of mechanical strength are of 
greatest concem., since these properties influence 
structural resistance. 
Thermal stresses may be. present -in 
either unrestrained or restrained structures or 
elements thereof if thermal gradients, which 
usually are caused by increased temperatures, 
are present. The magnitude of the stresses 
depends upon the amount of expansion the 
structure experiences, which depends on the 
value of the coefficient of thermal expansion. 
a,. for the material. The dimensions of this co- 
efficient are length per length per temperature 
unit. often expressed as in./in.°F or cm/cm/°C. 
The value of a. depends on the temperature. as 
shown in Figures 9-12 through 9-15 for a num- 
ber of materials. Use of an improper value of a, 
can cause serious errors in the estimation of the 
thermal! stress level in some materials. For 
example. the use of the room temperature value 
of a, for AIS] 30] stainless steel when the 
value of 1200°F (~650°C) should have been 
used. would result in a value of thermal stress 
that 1s 20 percent too Jow. For 2024 Aluminum 
Alloy. the error would be about 18 percent, 
whereas for Ti-8Mn Titanium Alloy. the error 
would be 45 percent. Thus. for large tempera- 
ture rises. the assumption of a constant value for 
the coefficient of thermal expansion could lead 
to significant errors. However. for smal] tem- 
perature rises. on the order of 200°F (110°C). 
the error would not be significant. For most 
metallic materials. the error in thermal. stress is 
about 2 percent per !00°F or 3 percent per 
100°C temperature rise above room temperature 
values of a,. Additional data may be obtained 
from ‘“‘Metallic Materials and Elements for Aero- 
space Vehicle Structures” (see bibliography). 
Changes in mechanical strength are the 
prime reason for concern about the effects of 
absorbed thermal radiation. Data from several 
sources are incorporated into Figures 9-16 
through 9-18 to provide some indication of the 


effects of elevated temperature on the tensile 


‘strength properties of 2014-T6 Aluminum 


‘Alloy, Titanium Alloy 6AL-4V, and 301 full. 
hard stainless steel. The various curves in these 
figures provide some indication of the effect of 
elevated temperature and also the effects of the 
speed at which the sample was heated and the 


- tinge it remained at a given temperature. 
on The data that will be discussed comes 


rom several sources. Since the techniques for 
testing that were used in these sources differ 
considerably, it will be of value to discuss the 
various techniques employed. 

The four basic types of loading tests that 
nave been used to determine mechanical prop- 
erties at elevated temperatures are illustrated in 
Figure 9-19 and 9-20. Figures 9-19a and b illus- 
trate an adaptation of the conventional tensile 
test for elevated. temperature. The sample is 
heated rapidly to’ an equilibrium temperature. 
At some finite time after thermal equilibrium. 
called the soak time. a conventional tensile load- 
ing test is conducted with a constant strain rate 
loading mechanism. Two loading rates are shown 


_ in Figure 9-19a. Figure 9-19b illustrates the type 


of results that are obtained from the type of test 
illustrated in Figure 9-}9a. 

Figures 9-19c and d illustrate another 
ype of tensile test where the sample is loaded 
with a constant load, then heated at some rate 
until rupture occurs. Figure 9-19d illustrates the 
type of results that are obtained from the type 
of test illustrated in Figure 9-19c. Figure 9-20 
illustrates two types of short time creep tests 
and their results. 


L Results from. these types of tests are 
shown for three aluminum alloys in Figures 9-21 


through 9-23. Some data shown in. Figure 9-21 
were obtained by the type of the test illustrated 
in Figure 9-19a. The soak time of the sample 
was one-half hour before. a standard constant 
Strain rate of 0.00033 in./in./sec tensile test was 
performed. Other data in Figure 9-21 also were 
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Figure 9-15. BB Coefficient of Linear Expansion as.a Function of Temperature. 
for Aluminum Alloys 
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Figure 9-16, Typical Stress-Strain Curves for Titanium 6AI-4V Alloy at 
Room and Elevated Temperatures eS 
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Figure 9-17, Typical Stress-Strain Curves for 2014-T6 Aluminum Alloy’ at 
Room and Elevated Temperatures eS 
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Figure 9-19. | Types of Rapid Heating Tensile Tests P| 
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Figure 9-20, - Types of Short Time Creep Tests | 
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obtained in the manner illustrated in Figure 
9-19a, but with the equilibrium temperature 
being reached in 10 sec and with a 10 sec soak 
time, and with three different strain rates used 
in the tensile tests. Other data were obtained in 
the manner illustrated in Figure 9-19c with a 
heating rate of 60°F/sec. The final data in Fig- 
ure 9-21 also were obtained in a manner similar 
to that illustrated in Figure 9-19c with a heating 
time and time of rupture edch approximately 10 
se 

The ultimate strength data in Figure 
9-2] tend to show the importance-of tempera- 
ture and strain rate, but the data from the differ- 
ent sources are not consistent and are, therefore, 
difficult to evaluate. For example, the 0.00005 
in./in./sec strain rate curve for 10 sec soak times 
from one source lies close to the 0,00033 
in./in./sec strain rate curve for one-half houi 
soak times from a different source. Interpolation 
between the strain rate data for 10 sec soak 
times would indicate that a 0.00033 in./in./sec 
strain rate curve for a 10 sec soak time would lie 
slightly to the right of the 0.00005 in./in./sec 
curve and, therefore, slightly to the right of the 
one-half hour soak time curve. However, there 


are ample data, as shown. in Figures 9-16. 


through 9-18, to justify the expectation that the 
10.second curve should fall significantly to the 
right of the one-half hour curve; that is, for a 
given temperature and strain rate, the shorter 
the soak time, the smaller the decrease in tensile 


strength. It should also be noted that constant. 


load curves should agree more closely. 

The ultimate strength data in Figure 
9-22 suggest a similar inconsistency although the 
data were taken from two different types of 
loading tests. The one-half hour data is from 
loading tests similar to that shown in Figure 
9-19a, and ‘‘constant load, heating time as indi- 
cated” data is from loading tests similar to that 
shown in Figure 9-19c. The latter data clearly 
shows the influence of heating times or heating 
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rates on rapid heating, constant load tensile test 
results, - . 


. The yield strength data shown in Figure 
9- or 2014-T6 Aluminum Alloy is somewhat 
more consistent. The instantaneous heating data 
were obtained in a manner similar to that illus- 
trated in Figure 9-19a, although times are very 
short. Specimens were heated instantaneously 
by an-electron beam, then a stress wave was sent 
up the specimen. An analysis of the stress wave 
characteristics on either side of the heated 
region provided yield and ultimate strength data 
as well as elastic modulus data. The times indi- 
cated for these data are the delay times between 
time of heating and the arrival of the stress 
wave. Thus, the data are for instantaneous heat- - 
ing, very short soak times, and extremely high 
strain rates (~10°9 sec”! ). The comparable strain 
rate data for 1/2 hour and 10 sec soak times are 
in the relative positions that would be expected. 
The one factor that may explain why the com- 
parable strain rate data of Figures 9-21 and 9-23 
are relatively close together despite the different 
soak times is the heating rate. 


Similar data from a single source are 
shown in Figures 9-24 and 9-25 for 6AL-4V 
Titanium Alloy and full, hard 301 stainless steel. 
The data in these figures also illustrates the 
inconsistencies in data discussed for Figures 9-21 
through 9-23. 


@ One other mechanical strength property 
of interest is the elastic or Young’s modulus of 
the material. One source indicates that for soak 
times of from 1/2 to 10,000 hours, the elastic 
modulus depends only on temperature and not 
on soak time for a number of different ma- 
terials. A separate source indicates that this is 
true for 6061-T6 and 2014-T6 Aluminum Alloys 
with soak times of 0.2, 1, and 10 msec. Al- 
though this conclusion. has not been validated 
for other materials at small soak times, it ap- 
pears reasonable to assume that the elastic 
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Figure 9-25, \_d Uitimate Strength Test Results for Full Hard 
Stainless Steel at Elevated Temperature 
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modulus only depends on temperature for 
metallic materials. 

In evaluating the information available 
on Changes in mechanical strength properties of 
materials, it should be. reemphasized that a 
variety of testing techniques have- been used to 
obtain the data presented in this section. Con- 
Stant strain rate and constant load techniques 
have both been used. A variety of heating tech- 
niques have been used, including ovens, resis- 
tance heating. inductive heating. hot-fluid bath 
heating. radiant heating. lasers, and electron 
beams. A considerable amount of data have been 
ined on a number of different materials. 
The data indicate that there are five 
basic variables that influence the degradation of 
mechanical strength properties at elevated tem- 
perature for particular materials: initial tempera- 


ture. heating rate, final temperature, soak time. — 


and strain rate. Not much is known about the 
effect of initial temperature since most data 
report “room temperature” as the reference 
temperature. Increases in heating rate. final tem- 
perature. or soak -times cause a decrease in ulti- 
mate and yield strengths. An increase in strain 


rate at a particular final temperature apparently . 


increases the strength. Because of the different 
testing techniques employed, the relative impor- 
tance of these variables in detérmining material 
strengths is difficult to determine: however. it 
appears that soak time and final temperature are 
most important, followed closely by strain rate, 
and finally heating rate as long as final tempera- 
tures reached in less than one minute. 

Although there are considerable data 
available on rapid heating effects, it is difficult 
to extrapolate to other than test conditions with 
any confidence. In any event, in estimating 
strength, data available on rapid heating and 
short times appear to be more appropriate than 
those for slow heating. long soak times for the 
general class of problems of interest in this 
manual. 


WB resistance To toa 


The ultimate concern of combined 
thermal/blast effects is the change, if any, in the 
vulnerability/survivability of tactical systems 
compared to vulnerability/survivability to ther- 
mal or blast effects separately. 

Figure 9-26 illustrates typical experi- 
mental stress-strain curves for 2014-T6 Alumi- 
num Alloy. with two approximations to this 
curve: the more or less classical straight line: and 
the so-called Bell’s stress-strain law for dynamic 
deformation: 


o= B ) - a el? 
ie : 


stress, psi, 


€ = strain, in./in., 

T = temperature of interest, °K. 
7, = melting temperature, °K, - 

6 = parabolic coefficient. 


The area under each of these curves is equal to 
the strain energy absorbed by the material. For 
the straignt line approximation the area under 
the room temperature curve equals about 433 
in. + Ib/in.? to a strain (e) of 0.01 in./in. This 
compares well with the area of 432 in. + Ib/in. - 
under the typical curve to the same e. The Bell 
stress-strain law gives an area of 441 in. + Ib/in.?. 
East The change in allowable uniform equiva- 
ent static load that will give the same maximum 
deflection as a uniform dynamic load has been 
determined for degraded 2014-T6 Aluminum 
alloy properties using elastic theory for large 
deflections of simply supported square plates.* 
The allowable load was defined as the maximum 


B. air blast loading of long duration, or a flat top wave, 
ssumed, 
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Figure 9-26. || Stress-Strain Approximation for 2014-16 
Aluminum Alloy 


static load the plate could carry without exceed- 
ing the yield stress. One-half hour soak time data 
were used to determine the strength properties, 
and a 20-in..square plate. 1/4-in. thick was 
assumed. The results are shown in Table 9-9. As 
can be seen, the maximum allowable load is 
almost a linear function of allowable stress. The 
explanation of this can be seen in the following 
equation, which is the elastic stress equation for 
the plate under study. 


o = 7,356.8p + 4.632E13 p23, 


The second term of the equation is for. mem- 
brane stresses and is the only term that contains 
the temperature dependent Young’s modulus. 
As the temperature increases. the value.of E de- 
creases. so the contribution from the second 
term is further decreased. and a fairly linear 
dependence on the pressure, p, remains. 

The allowable stress in the above 
example. was taken from 1/2 hour soak time 
data. This provides a lower bound for the allow- 
able load. Examination of Figure 9-21 through 
9-25 indicates that the rapid heating. short soak 
time data generally fall above the 1/2 hour soak 
time curves. More reliable. estimation of allow- 
able loads will depend on the availability and 
appropriateness of the material data that is re- 
quired for the type of analysis and calculation 


techniques being employed by an investigator. 
The truly dynamic case may be analyzed 


similarly. In both cases, however, the basic fac- 


tor that influences structural degradation is 
allowable stress as long as the elastic. analysis 
approachisused. - , : 
Where deformation or plastic response is 
responsible for failure, allowable stress no longer 
can be used to determine allowable load; either 
total deformation or total energy absorbed must | 
be.psed. 
The energy approach consists-of stating 
at the energy absorbed by a structure as it 
deforms is proportional to the product of the 
overpressure and the impulse values, Ap/, of the 
blast wave that causes the deformation. Thus. 


for a constant value of absorbed energy used as a 


measure of damage, any combination of over- 


“pressure and impulse that yields the same prod- 


uct will cause the same amount.of “damage.” 

The effect of elevated temperature on 
the damage values has been investigated. One 
source indicates that for face-centered cubic ma- 
terials, such as aluminum, the value of the para- 
bolic coefficient of Bell’s stress-strain law varies 
linearly with the temperature. For 2014-T6 
Aluminum, the equation would be 


o = 9.65 x 105? 


Table 9-9. ' Allowable Uniform Static Load for Aluminum 
Piate at Three Temperatures 
Young's Allowable 
Temperature oy Modulus, E Deflection Load 
(CF) (10? psi) (10° psi) (in.) (psi) 
70 57 0.323 7.25 
350 43 9.77 - 0.26) ' $.45 
500 17 0.117 2.10 
9-65 


where the value of 8 was obtained by fitting the 
curve to pass through the room temperature 
yield stress for the material. The Bell’s law equa- 
tion indicates that the energy under the stress- 
strain curve is a linear function of temperature 
by varying values of 6 linearly. This infers that 
the damage value, Ap/, should depend on a 
linear function of temperature, i.e., 


em 
i. 
Apl = (A 
(Ap Mr ( PDT, T, 
T. 
where 
(Ap); = damage value at temperature 7, 
and 
(Appl), = damage value at temperature Fs 
; 


If this relation were true, it would be a simple 
matter to evaluate the vulnerability of a struc- 
ture. Once the value of Ap/ had been deter- 
mined for ambient temperatures, the Ap/ value 
at elevated temperatures could be determined 
easily. 
The Bell’s law equation also indicates 
at yield strength, ultimate strength, and 


Young’s modulus are linear functions of tem- . 


perature for dynamic deformation. The data 
that were used to develop. the equation were 
obtained from impact tests of oven heated 
samples. No. information was provided as to 
what the heating rates or soak times were. The 
only other comparable data, at least in the dy- 
"namic sense, were shown in Figure 9-24, which 
indicate that data from dynamic tests do seem 
to approximate a straight line to some degree. 
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Thus, with further. testing and research, it may 
be that a relatively simple approach to assessing 
structural vulnerability to combined thermal/ 
blast effects will be possible for dynamic 
deformation damage criteria. 

The previous discussion has been con- 
cerned primarily with unstressed, unrestrained 
structural elements. There are three main struc- 
tural conditions that are encountered in the 
analysis of structural response to. combined 


thermal/blast effects. These are: unstressed 


members with no edge or end restraints; un- 
stressed members with restraints; and prestressed 
members with restraints. These conditions. refer 
to the status of a structure before either thermal 
or blast loading. The first condition is the easiest 
to analyze since there are no restraints to dimen- 
sional changes during or after heating so the 


thermal degradation of mechanical strength © 


properties is the main concern. The second con- 
dition is more complex since dimension changes 


are prevented in varying degrees and appreciable _ 


thermal stresses may exist at the time-of blast 
loading. The inclusion of these thermal stresses 
is. not uncommon in structural elastic analysis. 
Thermal stresses caused by thermal gradients in 


the structural member. also are included under 


this condition. The inclusion of stresses induced 
by the thermal gradient in an elastic structural 
analysis generally will require the use of a finite 


element structural analysis computer code.. 


Unfortunately, it cannot be assumed that either 


‘ form of thermal stress will be inconsequential, 


and they should both be addressed specifically 

in the solution of a particular problem. The final 

condition is the most difficult to analyze be- 

cause a member may be stressed prior to thermal 
r blast loading. : 

The final factor to consider in determin- 

ing structural response is the possibility of 


changes in failure modes, particularly. materials ° 
rupture. Visual examination of a 1/4-in. by > 


0.040-in. 2024-T3 Aluminum specimen that fail- 


ee 


key, 


ed indicated little necking down at the failure 
point. This may be due to the small area (about 
17 mm diameter) exposed to thermal radiation. 
Failure occurred on a 45° plane, with no cup- 
ping. Postshot photos of cylindrical/aluminum 
shells. exposed to HE air blast after rapid heat- 
ing. indicated an almost brittle-like fracture 
mode of failure. Failure of similar shells at nor- 
mal temperature under the same type of loading 
was Of a ductile deformation type. The possi- 
bility of a change in failure mode under dynamic 
loading at elevated temperatures has not been 
investigated. but is of concern since residual load 
carrying capabilities of deformed structures are 
sometimes used in determining vulnerability. 
om: summary. if elastic criteria are used to 

etermine survivabilityv/vulnerability, and if the 
necessary materials strength properties are 
known, a fairly, straight forward but not neces- 
sarily simple stress analysis may be performed. If 
dynamic or deformation criteria are necessary. 
the problem becomes very difficult. if not im- 
possible. because of the lack of validated anal- 
ysis techniques. A third criterion, material melt- 
ing. may be used to bound the problem: A tem- 
perature criterion may be chosen that is o 
nearly equal to the melting point of the materia 
with the inference that. at that temperature. the 
material has no strength, and therefore. no resis- 
tance to load. This latter criterion has been used. 
for example. to establish kill criteria for Army 
aircraft. 


SECTION V 


GE X-ksy pawsce errects 
GER tRocuction gy 
es Nuclear weapons aS X-ray sources 
and the environments they produce are de- 
scribed in Chapter 4. This section discusses the 


interaction of the X-ray environment with aero- 
space systems and some typical damage mech- 


anisms that result from the interaction.. The 
damage mechanisms depend on the fluence and 
energy. or temperature, of the X-rays and on the 
protective and structural materials involved in 
the interaction. For purposes of discussion the 
damage mechanisms can be divided conveniently . 
into the effects caused by hot X-rays and those 
caused by cold X-rays. Some familiarity with the 
materia! discussed in Chapter 4 is assumed in the 


following discussion. 
a: X-rays (typically 1 to 3 keV 
ack body temperatures) are absorbed in a thin 


surface layer. At sufficiently high fluence. a 
short pulse of X-rays can heat the surface rapid- 
ly and may. cause it to vaporize and blow off. 
This results in: (1) an impulse imparted to the 
total structure; and (2) generation of a strong 
shock wave that propagates into the structure, 
and which may cause spallation of material at 
free boundaries and internal fracture of ma- 
terials and bonds. These latter effects are pro- 
duced by shock wave propagation through the 
thickness of a surface structure such as the ther- 
mal protection shel! of a reentry vehicle. The 
former effects may produce damage by whole 
icle modes of response to the net impulse. 
The hot X-rays are more penetrat- 
ing. They can cause: (1) thermally. generated 
shock waves in the vehicle structural materials 
and internal components; (2) melting and vapor- 
ization of the substructure: (3) internal deposi- 
tion of energy in electronic components produc- 
ing transient or permanent damage (see Chapter 
6 and Section 7 of this chapter): or (4) produce 
internal EMP signals (see Chapter 7). 

While some nuclear weapons emit 
only co rays, all hot X-ray weapons have a 
cold component. Hence, for exoatmospheric 
events the hot X-ray effects are accompanied by 
cold X-ray effects. On the other hand, for endo- 
atmospheric explosions, the cold X-rays have 
short mean free paths, and the X-ray effects be- 
yond distances of a few tens of meters are. pro- 
duced by hot X-rays alone. 
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Damage to a particular. system depends 
on the. structural details of the system as well as 


the X-ray spectrum. Accurate damage prediction 


even for the simplest geometries, generally re- 


quire elaborate computer programs to calculate 
the energy deposition and the response of the 
system. Ali vulnerability analyses follow similar 
computational steps: 

1. The X-ray energy deposition is computed 
using known processes for the materials and 
structure. This energy is assumed most often to 
be deposited instantaneously. 

2. From the calculated energy deposition and 


the equation of state for the materials in the’ 


structure (if known) for the liquid, solid, and 
vapor phases of the material, a stress wave, 
which propagates through the surface structure, 
is calculated. 

3. Damage to the surface structure that re- 
sults from the stress wave (spallation, internal 
fracturing, delamination and debonding), is 
determined. These calculations generally are 
based on experimentally derived damage re- 
‘sponse curves. A good deal of progress has been 
made in the last few years in understanding and 
predicting fractures in most homogeneous and 
some composite materials. Even a small amount 
of fracturing may cause a significant loss in 
strength of material. One of the most complicat- 
ing factors involved in making reliable damage 


predictions is that of in-depth heating prior to 


stress wave propagation through the material. 

4. The response of the whole structure that 
results from the impulse imparted to it is deter- 
mined. This response also is generally based on 
experimentally determined damage or properties 
degradation measurements. In these computa- 
tions, the effects of the loss of strength of the 
structural material as a result of fracturing 
caused by the immediately preceding stress wave 
or by in-depth heating must be included. The 
structure may be broken by the impulse loading, 
or it may be weakened (buckled) so that it will 
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fail under subsequent stresses, for example, re- 
entry into the atmosphere. 


X-RAY ENERGY DEPOSITION 
CALCULATIONS & 


ee The starting point of all X-ray vulner- 
ability’ analysis is a calculation of the X-ray 
energy deposition. The determination of the 
transmitted X-ray fluence, or shine through, into 
the interior of the aerospace system, is implicit 
in these calculations. Calculation of X-ray 
energy deposition involves the use of energy 
dependent photon absorption cross sections and 


. energy and angle dependent scattering cross sec- 


tions. Corrections for fluorescence, i.e., reradia- 
tion of X-rays, may be significant for high 
energy photons in certain materials. Sophisti- 
cated energy deposition and transport calcula- 
tions, using detailed energy dependent cross 
sections and their angular dependence can be 
made by computers using Monte Carlo tech- 
niques. The results serve as an accuracy check 
and provide “effective” cross section data for 
simpler techniques. Approximation techniques 
greatly reduce the amount of effort and time 
involved in system analysis. Fairly accurate 
approximate values sometimes can be. obtained 
even by hand calculations, as will be shown with 
examples. 


me X-ray Cross Sections | 


The probability of a photon of energy 
hy traversing a distance of absorbing material x 
is e'#*, where p is the linear attenuation coef- 
ficient. This probability also can be written as 
e'(H/P)PX where w/p is the mass attenuation 
coefficient for the material (see paragraph 4-3). 
In this representation, p/p is in cm?/gm and px 
is the thickness in gm/cm?, i.e., the mass of ma- 
terial in the column of 1 square centimeter cross 


ection and x centimeters long. 
as: the monoenergic X-ray fluence inci- 


ent normal (perpendicular) to the material sur- 


face is y,. the direct fluence after traversing a 
thickness px of absorbing material is 


San = Foe MiPIP® cal/em?. 


The total mass attenuation coefficient 
p/p is the sum of several components that repre- 
sent the various mechanisms that can remove a 
primary photon from the direct flux, while it 
traverses a material. These mechanisms are de- 
scribed in paragraph 4-3. Chapter 4, and are: the 
photoelectric process. Compton inelastic scatter- 
ing. Compton elastic scattering (incoherent elas- 
tic scattering) and Rayleigh elastic scattering 
(coherent elastic scattering). The first two pro- 
cesses result in photon energy. being absorbed 
and the production of secondary electrons, the 
photoelectrons and Compton recoil electrons. 
respectively. The kinetic energy of these elec- 
trons is dissipated in the material and heats it. 
Energy deposition usually is calculated to have 
occurred at the depth at which the electrons are 
produced. However. in the case of very thin 
samples, the range of the freed electrons can 
exceed the material thickness, and energy de- 
_ position cannot be considered to be local. The 


energy removed from the primary photon beam 


in the Compton elastic process or the Rayleigh 
elastic process is not locally absorbed. A clear 
distinction must be made between attenuation 
of the primary flux and absorption of energy 
-from the primary flux. For this reason the total 
linear attenuation coefficient usually is written 
as w= yw, + w,. Here, , is an absorption coef- 
ficient. It represents the first two processes men- 
tioned above, which result in energy absorption. 
The second term, p., is an elastic scattering coef- 
ficient. which contributes to the attenuation of 
the primary flux but not to local energy absorp- 
tion. It is sometimes convenient to ignore the 
Rayleigh coherent scattering coefficient in the 
total attenuation coefficient. The coherent ra- 
diation has the same energy’, and nearly the same 


direction, as the primary photon and cannot be 
distinguished from it in the “bad” geometry 


situations that usually occur in nuclear effects 


aj lications.* 
os Mass attenuation coefficients for the ele- 


ments beryllium, aluminum, iron, copper, tung- 
sten, and uranium are given in Tables 9-10 
through 9-15, and Figures 9-27 through 9-32. 
respectively. These are representative of metallic 
materials used in aerospace systems. Mass atten- 
uation coefficients for ablator materials, carbon 
phenolic and tape-wound silicon phenolic are 
shown in Figures 9-33 and 9-34, respectively. In 
these tables and figures, Z is the atomic number. 
H,,/p is the coherent elastic scattering coef- 
ficient, B,./p is the incoherent Compton elastic 
coefficient, u,,/p is the inelastic Compton coef- 
ficient, and uw. /p is the photoelectric coefficient. 
AS designated previously, u,/p and y/p are the 
energy absorption coefficient and the total 
attenuation coefficient.’ 


9-34 X-ray Energy Deposition and 
Shine Through Fiuences 


X-ray energy deposition in a thickness 6 
at a depth x due to direct fluence photons is 
given by 


axa betes) ee 


If u,5 << I, and if y, is in cal/em?, this expres- 
sion can be written as 


: Ky (6 px 
A gin = (=) psd e =) 


cal/em?. 
Pp 


*Bi A “good” geometry is one in which the distinction between 
primary and scattered photons can be made accurately. 


The symbols XK, Ly Ly, etc., in the tables and figures indi- 
cate the binding energies of the various electron shells (see para- ; 
graph 4-3, Chapter 4), 
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Frequently, the absorption is written in terms of 
cal/gm by dividing out the thickness pé, 


x 
cal/gm. 


This expression for the absorption is in terms of 
a dose; however, this assumes that very little of 
the flux is absorbed in the deposition region at 
depth x, i.e., the deposition region considered is 
very thin. Clearly, more energy than is in the 


. incident flux cannot be absorbed. 
The equation for direct fluence (%,,,) - 


given in paragraph 9-32 can be used to represent 
a small energy band of photons in X-ray energy 
spectra such as those tabulated in Table 4-3, 
Chapter 4, for various black body spectra. The 


total energy in the direct X-ray fluence after 
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OE One ee ne ae ee ee 


traversing thickness x is obtained by summing 
over the energy bands. 


-(5) px 
y= 2) ve Pri cal/em?. 

i 

In a like manner, the total direct fluence X-ray 
energy absorption at depth px is obtained by 
summing for each energy band. 


pu 
rm a), O* 
A=2; (=) Yo, € tr) cal/cm?., 


Problems 9-3 and 9-4 illustrate how 
these equations can be used to calculate approxi- 
mate values for energy deposition and. shine 
through. 


OT peep aeipore ager 


Problem 9-3. Calculation of X-ray Energy Deposition at the Surface 


The. information contained in Tables 
9-10 through 9-15 and Figures 9-27 through 
9-34 together with the equation for energy 
absorption given in paragraph 9-34 provide the 
means to obtain the approximate X-ray energy 
absorbed for various spectra and several ma- 
terials. If the energy deposition at the surface on 
which the X-rays are incident is desired, the 
thickness. x, reduces to zero, and the energy 
deposition equation becomes 


A= DO w,/e), 9; cal/em. 


age 9-72. deleted. 


Related Material: See paragraphs 9-33 
and 9-34. See also Tables 9-10 through 9-15 and 
Figures 9-27 through 9-34. See also paragraphs 
4-1 through 4-3 and Tables 4-3 and 4-4, Chapter 
4, 
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Problem 9-4, Calculation of X-ray Energy Deposition at a Depth 
in a Material and the Shine Through Fluence 


\ Oe information contained in Tables 
9-10 through 9-15 and Figures 9-27 through 
9-34 together with the equation for energy 
absorption given in paragraph 9-34 provide the 
means to obtain the approximate X-ray energy 
absorbed for various spectra and several ma- 
terials at a depth x in the materials, as well as 
the fluence at that depth or the fluence emerg- 
ing from the back face of the material (shine 
through). 


roblem 4-1, Chapter 4, describes the method to obtain 
energy distributions for various source temperatures from the 
normalized Pianck distributions given in Table 4-1, Chapter 4, 


974 Fase 9-75 cleleted. 


eane 


Reliability: The reliability of the pro- 
® cedures described above depends on the X-ray 
spectrum and the absorbing material. The ac- 
curacy depends primarily on the relative impor- 
tance of the scattering cross sections, and, for 
some materials, fluorescence. No definite reli- 
ability estimate can be made. 


§ Related Material: See paragraphs 9-33 
and 9-34. See also Tables 9-10 through 9-15 and 
Figures 9-27 through 9-34. See also paragraphs 
4-1 through 4-3 and Tables 4-3 and 4-4, Chapter 
4. 


Table 9-10. | X-ray Cross Sections, Beryllium 7 = 4 (cm? /gm) z= 


Proton 
Eneres 
keV Uon'P H;./P Hip Hy io H,/P Hip 
Gt 7.097. [* 1.[54-3 2.770-7 1.507+4 1.507+47 1.507+4 
OLX 7.094.) 1.33)-3 3.560-7 1.109+4 1.109+4 1,109+4 
O11] 7,094- | 1.331-3 3.560-7 2,323+5 23a3td 2.323+5 
0.18 7.079. | .2.007-3 7257-7 1.105+5 1.105+5 1.105+8 
0:2 7.084-1 2.972-3 1.421-6 5.302+4 §.302+4 $.302+4 
0.3 6.984. ] 5.171-3 3.685-6 1.816+4 1.816+4 1.817+4 
0.4 6.889- | 7.660-3 7.223-6 8.29243 8.29243 8.29243 
0.5 6.770-1 1.039-2 1.218-5 4.456+3 4.45643 4,456+3 
0.6 6.63.1 1,332.2 1.866-5 2.66 2+3 2.662+3 2.66243 
0.8 6.306: | 1.972-2 3650-5 1.165+3 1.165+3 1.16543 
1.0 §.939-] 2.674-2 6.137-5 6.072+2 6.072+2 6.078+2 
13 - 4,986.] 4.648-2 1.591-4 1,820+2 1.820+2 1.825+2 
a0 4.133-] 6.19)-2 2.771-4 7.626+1) 7626+} 7.674+ | 
3.0 3.01 9-1 7946-2 5.309-4 2,199+| 2.199+ | 2.237+1 
4.0 2.37 1-1 9.009-2 7.9) 1}-4 9.004 9.004 9.332 
3.0 | .968-] 9.83}-2 1.066-3 4.481 4.482 4.778 
6.0 1.685-} 1.053-1 1.356-3 2.527 | 2.§28 2.802 
8.0 }.285-] 1,170-] 1.97]-3 1.019 1.023 1.267 
10.0 1.003-] 1,26}-] 2.613-3 5.024-] 5.050-] 7.314-] 
15.0 5,780.2 1.386-] 4.173-3 1.384-] 1.426-1 3.390-1 
20.0 3.677-2 1432-1] 5.589-3 5.539-2 6.098-2 2.410} 
30.0 1. 844-2 ].433-] 8.064-3 1,525-2, 2.331-2 1.851-] 
40.0 ].098-2 1.395-] 1.017-2 6.126-3 1.630-2 1.668-] 
$0.0 Ta51e2 1.348-1 1.199.2 3.027-3 1,502.2 1.570-] 
60.0 §.]32-2 1.299-] 1.357-2 1.705-3 1,528.2 1.503-1 
80.0 2.947-3 1.207-1 1.6}7-2 6.930-4 1.686-2 1.406-} 
100.0 1.898-3 1.126-1] 1.8]9-2 3.464-4 }.854-2 1.331-1 
150.0 8.554.4 9.646-2 2.164-2 9953-5 2.174-2 1.191-] 
200.0 4.81 1-4 8.476-2 2.364-2 4.159-5 2.368-2 1,089. ] 
300.0 2.138-4 6.877-2 2,.$59-2 1.241-5 2.560-2 9,458.2 
400.0 1.203-4 5.830-2 2.630-2 5346-6 2.63 1-2 8.473-2 
500.0 8.019.5 5.080-2 2.645-2 2.811-6 2.645-2 7733-2 
600.0 5.346-5 4.546-2 2.634-2 1.675-6 2.634-2 7.156-2 
800.0 2.807-5 3.708-2 2.571-2 7501-7 2.571-2 6.282.2 
1000.0 1.737-5 3.157-2 2.487-2 4.074-7 2,487-2 $.646-2 


ah 


7.097 -1=7.097x 107, 
1.507 +4 = 1.507 x 104, 
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Table 9-11. 6 X-ray Cross Sections, Aluminum Z = 13 (om?/grm) 


Photon 
Energy i 
keV ,,/0 H;./p hy ./o Uplp ule u/p 
0.1 2.504 5.614-4* 1.348-7 4.462457 4,462+5 4.46245 
OV T8L, 2.503 7.044.4 1,998.7 2.849+5 2.849+5 2.849+5 
0.118 2.503 7,044.4 1.998-7 4.27345 4.273+5 4.27345 
0.15 2.501 9.8724 3.568-7 2.192+5 2.192+5 2.19245 
0.2 2.496 1,473.3 7.102-7 9.926+4 9.92644 9.926+4 
3 2.482 2.59)]-3 1,874.6 3.250+4 3.25044 3.25044 
0.4 2.463 3.867-3 3.706-6 }.472+4 1.47244 1.472+4 
0.3 2.439 5.276-3 6.286-6 ‘7.96243 7,96243 7.96443 
0.6 2.411 6.799-3 9.696-6 4.8}9+3 4.819+3 4.82243 
0.8 2.345 1.015-2 1.917-5 2.183+3 2.18343 2.18543 | 
1.0 2.269 1.384-2 3.251-5 1.181+3 1.181+3 1.183+3 
1.5 2.063 2.43 }-2 8.474.5 3.866+2 3.866+2 3.887+2 
5.560K 2.036 2.533-2 9.169-5 3.473+2 3.473+2- 3.49342 
1.560 2.036 2.533-2 9.4169-5 '4.403+3 4.403+3 4.405+3 
2.0 1.869 3.295-2 1514-4 2.352+3 2.35243 2.35443 
3.0 1,557 4.629-2 3.164-4 8.087+2 8.087+2 8.103+2 
4.0 1.326 5.732.2 5.166-4 3.674+2 3.67442 3.688+2 
5.0 1.143 6.618-2 7.3814 1.960+2 1.960+2 1.972+2 
6.0 9.910-1 7.543-2 1.004-3 1.162+2 1.362+2 1.372+2. 
8.0 7,482. | 9.132-2 1.598-3 5.004+! - 5.00441 5.088+ } 
10.0 5.730-] 1.045-1 2.258-3 2.57141 2.57141 2.6394] 
15.0 3,272.1 1.228-1 3.861-3 7.47) VATS 7.925 
20.0 2.133-1 1.319-1 5.390-3 3.057 3.062 ' 3.407 
30.0 1.148-1 1.388-] 8.161-3 8.505-1 8.587-1 1.112 
40.0 7.319-2 1.394-1 1,052-2 3.394-1.. 3.499-17 5.625-1 
$0.0 5.136-2 1.374-} 1.256-2 1657-1, 1.783-] 3.668-] 
60.0 3.792-2 1.343-1 1.432-2 9199.2 1.063-] 2.785-1 
80.0 2.209-2 1.268-] 1.71922 3.625-2 5.344-2 2.024-1 
100.0 1.466-2 1.194-] 1.9432 1.759-2 3.702-2 1711-4 
150.0 6.855-3 1.035-} 2.324-2 4.728-3 2.797-2 1.383-1 
200.0 * 3.935-2 9,131.2 2,549.2 1,868.3 2.736-2 1.226-] 
300.0 1.775-3 7,425.2 2.775-2 3.092-4 2.826-2 1.043-] 
400.0 1.000-3 6.297-2 2.861-2 2.044-4 2.881-2 9.279-2 
500.0 6.384-4 5.490-2 2.880-2 1.0144 2.890-2 8.444.2 
600.0 4397-4 4.880-2 2.868-2 5.7495 2.874-2 7.797-2 
800.0 2.455-4 4.012-2 2.798-2 2.374-5 2.800-2 6.837-2 
1000.0 1.540-4 3.417-2 2.706-2 1.208-5 6.139-2 


$.614-4=5.614x 104, 


PJ 
Bs: +5 = 4,462 x 105, | 


2.707-2 - 


an 


Table 9-12, } X-ray Cross Sections, Iron Z = 26 (cm*/gm) 2 


Pristor 
Energy 
keV Kg b Hye & ui #, P He up 
Gs 442 2.9] 0-4" 6,9 96-8 1.030+67 1.036+6 1.036+6 
as 4.538 S.217-4 1.879-7 3.4 1345 3.4)3+5 3.413+5 
Or 4.833 7.890344 3,798.7 LSS5245 1.552+5 1. S53+8 
o3 4.8)7 patse3 1.023-6 S.}1St4 $.114+4 5.) 1544 
0.4 4.790 2.440-3 2.003-6 2.32744 232744 2.32744 
0.8 4.769 2950-3 3887-6 1.203+4 1.26344 1.26344 
0.6 4.739 3.835-3 5.522-6 7,065+3 7.005+3 7.67043 
0.7032 ; 4.69 4.807-3 8.224-6 4.87043 4.870+3 4.87443 
O.7US 4.690 4.867-3 §.224-6 1.74 }4+4 1.74144 1.74244 
G7. 4.69) 4.905.3 8,507.6 1.652+4 1.652+4 1.65344 
O72) 7 4.60) 4.995-3 8397-0 2.479+4 2.47944 2AT9e4 
Qs 4.05 4500-3 EPOS-§ 1.838+4 83544 1.83844 
Osteo? tide 6.203 } 265-5 1.So4t4 1.56444 1.So4+4 
Oade | 4.03% 6.2an-3 1.28.8 2.346+4 2.346+4 2.346+4 
ro 4.564 7903.3 1.§07-8 1.449+5 }.449+44 1.450+4 
1s 4.258 P43 1-2 $045.5 4.50543 4,805+3 4,509+3 
2.0 3,980 2.022-2 9.4528 1.96643 1,906+3 1.970+3 
3.0 3.414, 3.115-2 24a 6.) ]2+2 6.11242 6.146+2 
46 2.9)7 4.1 26-2 3 + 2.06>+2 2.H6d+2 2.6097+2 
3.0 AST §.035-2 S713 1.402+2 1.402+2 1.428+2 
6.0 2128 $835.2 7, 8.2934} 8.29342 8.81141 
O40IL AK 1.994 6.147-2 §.830-4 6.8734] 6.8734] - 7,079+] 
TI2K 1.796 6.649-2 1.085-3 5.08141 5.08 1+] 5.267+] 
7A12 1,796 6.649-2 1,055-3 4.19742 4.19742 4.21542 
8.0 1.594 7.257-2 1.28.3 3.09542 3.09542 3.11242 
10.6 1.24] §.321-2 1.832-3 1.71442 LL7i4+2 {.727+2 
15.0 7 6&4} 1.013] 3.263-3 $.622+1 5.622+] §.709+ | 
20.0 §.320-} LJ 18-1 4.688-3 | 2.483+] 7483+] 25484} 
30.0 2.93]-] 1211-1 7.31 h-3 7.632 7.639 6.034 
40.0 1828-1 242-1 9.633-3 3.258 3.268 3.575 
$0.0 1,242.1 1.240-] 1.168-2 1.675 1.68? 1.938 
60.0 8.9742 1.223-] 1.34.2 9.709-} 9.843-1 1.196 
$0.0 §,290-2 1.172-1 1.623-2 4.104-1 4.260-1 5.967-] 
100.0 3.479-2 J. 14-1 1,843-2 2.1 08-1 2.292] 3.754-] 
150.0 1593-2 9.800-2 2220-2 6.353-2 8.573-2 1.997-1 
200.6 9.058-3 8.699-2 2441-2 2.75}-2 $.192-2 1.480-1 
300.0 4.083-3 7.24.2 2,665.2 $699.3 3.535-2 1.106-1 
400.0 2.281-3 6.057-2 2,752.2 3.945-3 3.146-2 9,432.2 
500.0 1.457-3 §.284-2 2.778-2 2.175-3 2.996-2 8.425-2 
600.0 1.008-3 4.698-2 2.770-2 1.354-3 2.905.2 7.704-2 
800.0 5.618-4 3.863-2 2.707-2 6.590-4 2.773-2 6.692-2 
1000.0 3,548.4 3.291-2 2.619-2 3.863-4 2.658-2 §.984-2 


2.910-4=2.910x 107%, 
+P 1.036 + 6 = 1.036 x 10°. 
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Table 9-13, & X-ray Cross Sections, Copper Z = 29 (cm2/gm) eS 


Photon 
Energy 
keV Hoe /P Hi, /? uIP h,le u,/p ulp 
0.1 §.295 1.655-4* 3.975-8 3.00}4+5t 3.001+5 3.001+5 
0.1204, §.294 2.185-4 6.29)-8 1.976+5 1.976+5 1.976+5 
0.120 5.294 2.1854 6,291.8 2.17445. 2.17445 2.174+5 
0.15 §,292 3,090-4 1.132-7 1.293+5 1.293+5 1.293+5 
Q.2 5.288 4.814-4 2315-7 6.64544 6.64544 6.64644 
0.3 §.276 8.9924 6.503-7 2.60246 2,602+4 2.602+4 
0.4 5,259 1.401-3 1,350-6 1.338+4 1.338+4 1.338+4 
0.5 5.238 1.97§-3 2.380-6 7.98443 7.98443 7.99043 
0.6 §.213 2.615-3 3.7744 $.238+3 §.238+3 $.243+3 
0.8 5.15) 4.072-3 71,7684 2.69343 2.69343 2.69843 
0.931L §.302 5.145-3 1,138-5 }.896+3 - 189643 1,90}+3 
0.931 5.102 $.145-3 1.138-5 8.043+3 8.043+3 8.048+3 
O.951L, 5.094 5.314-3 1,201-5 7.$93+3 7.59343 7.598+3 
0.951 ~ 5.094 §.314-3 1,201-5 $.139+4 1,139+4 1,)39+4 
1.0 5.074 5.74]-3 1364-5 9.93243 9,.932+3 9.937+3 
1.097L ; 5.020 6.617-3 Lj 725-3 7.72543 7,725+3 7.73043 
1.097 5.020 6.617-3 1,722.5 1.159+4 1,.t59+4 1.15944 
1.5 4.840 1.071-2 3.782-5 4.93543 4.93543, 4.940+3 
2.0 4.569 1.578-2 7,372.5 2.25443 2.25443 2.25843 
3.0 4.006 2.566-2 1.775-4 7.465+2 7.465+2 7,505+2 
4.0 3.470 3.499-2 3.1914 . 3.40942 3.409+2 3.444+2 
5.0 2.987 4.352-2 4.938-4 1.856+2 1,856+2 1,886+2 
6.0 2.57) $.130-2 © 6.952-4 1.129+2 1, 12942 1,155+2 
8.0 1.934 6,489.2 1.155-3 $.156+1 $.156+1 $.356+1 
8 0482 ,K 1.92) 6.517-2 }.166-3 5.073+} 5.073+1 §.272+1 
8.979K 1.697 7046-2 £,399-3 3.765+1 3.765+) 3.94t+] 
8.979 1.697 7.046-2 1.399.3 2.97842 2.97842 2.995+2 
10.0 1.499 7607-2 1.676-3 2.234+2 2.234+2 ae DOt 2 
15.0 9.085-1 9.494.2 3.065-3 7.373+1 7:373+1 7473+} 
20.0 6.245-] 1.059-] 4,458.3 3285+) 3.285+1 3.35841 
30.0 3.479-] 1466-1 7.071-3 1.02641 1.026+1 1,073+1 
40.0 2.195-1 1.201-1 9.378-3 4.434 4.443 4,783 
50.0 1.502-1 1.203-] §.137-2 2.301 2.312 2.583 
60.0 1.090! 1187-1 1.309-2 1.344 5.357 1.585. 
80.0 6.473-2 1.140-1 1.587-2 5,736-1 5.895-1 7.678 | 
100.0 4.256-2 1,085.1 1,803-2 "2.963-1 3.143-1 4.654-1 
150.0 1.958-2 9.568-2 2.172-2 8.966-2 1,114-] 2.267-t 
200.0 1217-2 8.510-2 2.390-2 3.869-2 6,259.2 1.589-1 
300.0 §.014-3 6.974.2 2.609-2 1.204-2 2.813-2 1,329-] 
400.0 2.8213 $.933-2 2.696-2 $,.346-2 3.231-2 9.466-2 
$00.0 }.804-3 5.177-2 2.721-2 2.882-3 3.009.2 8.367-2 
600.0 1.249.3 4.603-2 2.715-2 1.754-3 2.890-2 7.619-2 
800.0 6.96 7-4 3.785-2 2.654-2 8.161-4 2,736-2 6.590-2 
1000.0 4407-4 3.224-2 2.568-2 4584-4 5,882.2 


1.655-4 = 1.655 x 107, 
tH 3.001 + 5= 3.001 x 10°. 


“2.614-2 


Presi 
Enon 
keV Boy HP Hye Hye uIP Hp 
Gi. }.19}43% 1.6 78-45 4.091-8 3.230+6 3.23046 3.23046 
os 1.1994) 2.987-4 }.070-7 1.089+6 1.089+6 1.089+6 
0.2 1.1907} 42)-4 2123.7 5.03748 ~ §.03745 §.037+5 
0.3 1.1884) 7.7914 5.609-7 1.69945 1.69945 1,6994+5 | 
0.5 1.18453 1.168-3 1.1206 7.857+4 | 7.85744 7.85844 
0.492%, L1&}o3 1553-8 1.84 1-6 4.520+4 4.520+4 4.82244 
0.492 118141 }.953-3 1841-6 4.97244 4.972+4 4.974+4 
0.5 1.1804] 1.590-3 1.915-6 4.78244 4.752+4 4.75344 
0.8984 1.17641 2.02»-3 2.906-6 2.98144 2.98144 2.98244 
0.505 1. e+i 2.0203 2.900-6 3.27944 3,279+4 3.280+4 
Ge j.7s=] 2.0523 2.96 5-6 3.20644 3.20044 3.207+4 
us 1 dese 3.Go0-3 8907-6 1483+ Lagats 1.454+4 
me J yae- 4186-3 1.00.5 &.)S943 8.15243 &.164+3 
if 1104+} 7.3703 2633.5 2 i4y-3 2.74943 2.7603 
Lage Lats} 9,219.3 3,957-8 1.06443 1.66443 1.67445 
} sce 1.70 | 9,2} 9.3 3987-8 4.18973 4,} 5903 4.17043 
petit 1003+] g.509.3 4.202-8 3.79943 3.79a-3 3.8Ga 3 
ieee 1,003- 9.899.3 4.262-5 5.69043 5.69643 £70743 
2G 1.08344 1.039.2 4.92.8 4.768+3 4.76843 4.77943 
eee 1.Qio7? J.198.2 6.459-8 3.38343 3.35243 3.30243 
22s! 1.O}o+] 1.19%-2 6.489.5 4.022+3 4.022+3 4.03343 
2: Ona }3ma.2 &.289.5 2.90743 2.907+3 2.91673 
oy G.d45 }.Aud-2 8,289.5 3.19743 3.19743 3.207+3 
ey Gols 1.8)G2- 9999.8 2.50743 2.50743 2.S1o#3 
a 9.6}3 1310-2 9,990.5 2. 7S7+5 DISTH+H 2.76773 
3, 9.457 Lele-2 1.336-4 3.33545 2.33543 4.34443 
4.0 8.47] 2103-2 2.0)13-4 1.080+3 1.08043 1.08843 
5.0 7 Sau 2o0$1]-2 3092-4 5 93742 5.93742 6.01342 
oc 6.79) 3.J64-2 4.3424 3.64242 3.64242 3.7102 
0 5.500 4.012-2 7,289.4 },684+2 1.684+2 1.7402 
10.0 4.873 4.7802 3.078-3 9.26041 1.260+! 9.7174) 
10.2672 . 4n4 4849-2 }.b}4-3 8.7604) 8.7664) QO .DS3+} 
10.207 © 4.424 4,849.2 144-3 1,.824+2 1,824+2 1.869+2 
11,8441. 3.884 238-2 1.3513 1,300+2 1.300+2 1.33942 
11.844 7 3.884 $238.2 1.35)-3 1.95042 1.950+2 1.989+> 
12.1001 , 3.694 5,392-2 1.453-3 1.7}3+2 1.733+2 1.75142 
12.)00 3.094 $.392.2 1.4533 2.570+2 2.5702 2.605+2 
15,0 2.936 1.162-2 2.0303 }.423+2 1.42342 1.45342 
20.0 2.071 FAY2-2 3.079-3 6.448+} 6.448+ | 6.6634} 
30.0. 1.227 §.2460-2 §.173-3 2.1134) 2.1134] 2.2454) 
40.0 §.196-1 8.822.2 7,145.3 9.575 9.582 1.049+1 
80.0 5.854] 9.076-2 8.918-3 $.182 $.19] 5.867 


Table 9-14, . | X-ray Cross Sections, Tungsten Z = 74 (cm? /gm) BD 


Table 9-14. (Concluded) 


Photon 

Energy 
he\ ; Hig Hie P 
S9.318L 4K 4. ] 9.) S)-2 
60.0 | 4.3731 9.1502 
69.5254 3. 1 9.1002 
69.525 3. ] 9.100-2 
80.0 25 ] 9,037.2 
100.0 1.743-] §.768-2 
150.0 §.216-2 &.023.2 
200.0 4,733-2 7,237.2 
300.0 2.44)-2 6.024.2 
400.0 1,209.2 5.160-2 
500.0 7.7503 4,523.2 
600.0 §.366-3 4.03)-2 
800.0. 3.014-3 3.323-2 
1000.0 1.9163 2.831-2 


LIS} +} = PELOd, 
1084-4 1.684 x 107, 
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3.238 
3.138 
2.092 


1.020+] ° 


7.367 


4.265 
1.471 
6.666- J 
2.144-] 
9.66 1-2 


5.287-2 
3.28 1-2 
1.606-2 
9.6103 


3.248 
3.148 
2.104 
1.021+} 
7.380 


4.280 
3.490 
6.875-} 
2.373-4 
1,203-} 


7674-2 
$.664-2 
3.940-2 


3.224-2 


- Table 9-15. aN X-ray Cross Sections, Uranium Z = 92 (cm? /gm) a 
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Table 9-15. (Concluded) 
Photon 
Energy ; 
keV Hig B HP hy .!P H,/p H,/p HIP 
60.0 §.339.} 8.374-2 9:690-3 6.261 6.273 6.885 
80.0 3.178} §.353.2 1.226-2 2.936 2.948 3.349 
95.4401 .K 2.176! &.189-2 1.416-2 1.700 1.714: 2.014 
100.0 2414-1 8.166-2 1.43}-2 ).632 1.646 1.939 
115.606K },046-] 7.956-2 1.560-2 Lala 1.130 1.374 
115.606 1.646-) 7.956-2 1.5602 5.035 5.046 §,295 
150.0 9.832.2 7573-2 1.799-2 2.166 2.184 2.358 
200.0 _ 5.6272 6.872-2 1.999-2 9.940-] 1.013 1,139 
300.0 2,533-2 §.628.2 2.150-2 3.873-} 4.088-1 4,904.5 
400.0 1,429.2 4,935.2 2.273-2 2.034-} 2.261 2.898: 1 
$00.0 9,349.3 4.3402 2.299-2 1.196-] 1.426 1.982.) 
600.0 6.343-3 3.869-2 2.292-2 7428-2 9.720-2 1.4221 
800.0 3.55743 3:128-2 2.200-2 3.107-2 §.307-2 8.792-2 
1000.0 2,264.2 2.725-2 2.1782 1.381-2 3.559.2 6.51}-2 . 


1.42341 = 14.23, 
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Figure 9-27, B Photon Cross Sections in Berytlium a 
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Figure 9-28. | Photon Cross Sections in Aluminum & 
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Figure 9-29, eS Photon Cross Sections in iron 


MASS ATTENUATION COEFFICIENT (cm*/gm) 


9-87 


Soe yeaa a 


hy (keV) 


10° 
=u/p 
S 3 
oa ee 
ae 
: 2 
b t 
z2 10° 2 
w w 
o © 
t ira 
Oo ° 
S S 
z z 
ie) Co . 
os - 
g 3 ) 
4 10 
ry _§ 
- 
‘= 
< q 
@ : 
< a 
= = 
i0* 
10 
100 1000 


hy (keV) 


Figure 9-30. & Photon Cross Sections in Copper f 
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Figure 9-31. | | Photon Cross Sections in Tungsten B 
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Figure 932. & Photon Cross Sections in Uranium a 
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Figure 9-33. 8 Photon 


Cross Sections in Carbon 
6.08 & 


Phenolic (CP}; Z = 
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Figure 9-34, Photon Cross Sections in Tape Wound 
Silicon Phenolic {TWSP}, Z = 9.01 
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9-35 X-ray Energy Deposition 
Summary S 
The methods described in paragraphs 


9-33 and 9-34 and illustrated in Problems 9-3 
and 9-4 allow the calculation of curves that 
show approximations of energy deposition as a 
function of depth for black body spectra inci- — 
dent on any material, if the cross ions are. 
known for the material 


INITIAL PRESSURIZATION OF 
IALS DUE TO. 
X-RAY DEPOSITION Qe 


An immediate consequence of the 
deposition of X-ray energy is the rapid heating 
of the material. This heating causes an initial 
pressure distribution as a function of depth in 
the structure. The initial pressurization generates 
shock waves that propagate through the thick- 
ness of the shell of the structure. The heating 
can result in a solid material changing phase. 
that is, melting or vaporizing. The melting and 
vaporization cause blowoff, which imparts an 
impulse to the structure and excites whole struc- 
ture modes of response. 


9-36 Phase Changes Induced. 

by X-ray Heating 

In most nuclear weapon X-ray environ- 
ments, the X-ray energy is deposited in a very 
short time, a few nanoseconds to a few hundred 
nanoseconds. The material cannot expand 
appreciably during this time, so the energy 
deposition process can be considered to occur at 

a constant volume or at normal material density, 
P,- Rapid. melting and vaporization are accom- 
panied by enormous pressure increases. Values . 
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Figure 9-37. Energy Deposited by Black Body 
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for enthalpy changes for melting and vaporiza- 
tion for the metals discussed in the previous sub- 
section are given in Table 9-]7. These values are 
for one atmosphere pressure. In most X-ray 


problems of interest the material is initially at. 


very high pressure. so these values can be con- 
sidered to be only approximate. This approach is 
not correct for ablators as a class although it 
might apply to carbon phenolic. in a cold en- 
vironment. Confining the discussion to metals 
will not restrict the transfer of principles:* 

The rising pressure that results from 
heating at constant density is illustrated in Fig- 
ures 9-39 and 9-40 where isoenergy lines of 
aluminum are shown in pressure-density plots. If 
the internal energy is above tne critical energy, 
3,016 cal/gm for aluminum, the material can be 
considered as a vapor. Figure 9-40 shows the 
high pressure, high energy intercepts with the 
normal density abscissa (p, = 2.7 gm/cm?), The 
release adiabats for expansion from density p, 
to low density and pressure also. are shown in 
this figure. Expansion along the adiabat results 
in decreasing internal or potential energy as the 
material develops kinetic energy during “blow- 
off.” For example. a 6.000 cal/gm energy depo- 


sition in aluminum at | eae oe gm/cm? results in 
a pressure. of about 1.5 megabars (Mb). The 
aluminum would expand from that state to low 
pressure and density, with final internal energy 
of about 3,000 cal/gm and about 3,000 cal/gm 
of kinetic energy. The 3,000 cal/gm of internal - 
or potential energy is used to overcome the 
physical and chemical forces that bind the atoms 
together in the solid. This leads to the concept 
of heat of sublimation. The heat of sublimation 
at absolute zero, £,.,, is the energy required to 


’ form the saturated vapor from the solid at a 


temperature of absolute zero, Thus, £., does 
not include any energy of kinetic motion. The 
energy of sublimation generally is a function of 
temperature becoming larger for larger deposi- 
tion energies (temperatures). 


§.. problem of phase changes. in a composite heat shield. 


ablator is more complicated since different deposition profiles, 
material enthalpy, and thermal conductivities are involved in the 
calculations. While some. materials, e.g., tape-wrapped carbon 
phenolic, may behave like metals in a cold environment, the 
techniques described here generally are not applicable to the 
description of the blowoff process in the broad category of com- 
posite materials that use three dimension (3-D) weaves for heat ° 


‘shields or for X-ray shields that use dispersed high Z materials 


for loading. : : 


Table 9-17. | Enthalpy Change for Selected Metals 


(cal/gm) 

Metal- 

Atomic To Through To Through ~ Sublimation 
Weight Melt "Melt Vapor Vapor Energy 
Be 9.013 876.0 1,187.0 2,147.0 0,040.0 8,682.0 
A] 26.98 160.4 255.3 771.1 3,347.0 2,893.0 
Fe 55.85 250.8 315.8 $73.0 2,071.0 1,782.0 
Cu 63.54 110.0 160.0 336.0 1,481.0 1,275.0 
W 183.85 153.9 200.0 304.0 3,353.0 1,110.0 
U- 238.00 , 49.0 64.5 


171.9 596.1 492.1 
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Figure 9-39. Aluminum fsoenergy Lines. 
Parameter is Energy in cal/gm 
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Figure 9-40. e Aluminum Isoenergy Lines and Adiabats = 


9-37 The Griineisen Parameter_and 
the Equation of State 


As a first approximation to the equation: 
of state of material, it is assumed that the pres- 
sure in the material increases linearly with the 
deposited internal energy per unit volume, 


P= Gre, 


where G is the Griineisen ratio for the material, 
€ is the internal energy per unit volume, and 7 = 


pip,, the ratio of the density to a normalized 


density of the material, usually the ambient 
density. In codes used to. calculate shock wave 
propagation, for example the PUFF type codes, 
more elaborate equations of state are used to fit 
the experimentally determined behavior of the 
solid and vapor phases. The equations used in 
the PUFF codes reduce to the equation given 
above when 7 = 1, and generally only one vans 
of G is used to specify a material. 

Although there may be small errors in 
calculating the X-ray energy deposition as a 
function of depth, the significant sources of 
errors in pressure predictions are the accuracy 
and validity of the Griineisen parameter for solid 
and vapor equations of state. Experimental 
values of G obtained by different methods result 
in factors of about two uncertainty even for 
some of the common materials such as alumi- 
num, beryllium, and tungsten. Even though indi- 
cations are that G is not well known as a func- 
tion of deposited energy, some very good corre- 
lations have been obtained for computed and 
measured values of pressure waves in X-ray tests. 
when careful calculations are made employing 


The units for € in the equation given 
above are energy per unit volume, which have 


the same dimensions as pressure. Therefore, the 


energy required in cal/gm for a phasé change can 
be expressed in units of pressure, if the density 
of the material is specified. If the internal 
energy, £, is given per unit mass, the relation to 
€ is 


€ (cal/em3) = E (cal/gm) p, (gm/cm? ). 


The value of € in megabars may then be ob- 
tained by the relation 


€ (Mb) = € (cal/em3) x 4.18 x 107 (2) 
dyne * cm Mb \ 

x | (————] x 1 — 
erg \i0!? dyne/cem?/ 


4.8 x 10% ¢ Se 
em? 


5 cal 
4.18 x 10° p,E (=). 


e (Mb) 


Thus, the previous equation for pressure may be 
written 


: P = 
elastic plastic properties of the materials. (Mb) Fe = ee); 
or 
= Initial pressurization in distended ma- 2 5 
terials such as porous metals and foams present a P (Mb) = 4.18 x 10° Gp E (cal/gm). 


particularly difficult and uncertain condition for 
current analytical techniques. 


The enthalpy changes of the metals shown in 
Table os 17 in cal/gm are given in Table 9-19. 


Table 9-19. = Enthaipy Changes € (Mb) S 


ane To Through To Through Sublimation 
Metal (gm/cm?) Melt Melt sa Vapor Energy, £, 
Be 1.85 0.068 0.0918 0.166 | 0.776 "0.671 
Al 2.70 0.0181 0.0288 0.087 0.378 0.326 
Fe 7.86 0.0824 0.1036 0.188 0.680 0.585 
Cu 8.92 0.0410 0.0596 0.125 0.552 0.475 
W 19.3 0.124 0.161 0.245 1.092 © 0.895 
U 18.7 0.0383 0.0504 0.134 ' 0.466 0.385 
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Change, P. (Mb) a ; 
= 1) 


Tabie 9-20, Pressure 
{0 = py. 
To Through: To Through Sublimatio:: 
Meis! G Melt Melt Vapor Vapor Energev.£. 
be is 0.090 0.133 0.24) 0.12 0.973 
A rome 0.0356 0.0615 0.185 0.804 0,694 
Fe 1.69 0.139 0.178 0.318 135 0.989 
C 2.00 0.082 0.119 0.250 1.10 0.930 
W j43 (0157 0.23 0.350 1,56 1.28 
U 2.03 0.078 0.102 0.273 0.782 


0.946 


The. pressures assouiited with these changes at 
‘ambient density. ie. when p = p,.and P (Mb) = 
Ge (Mb). are shown in Table 9-20. , 

From Table 9-20. aluminum has a sub 
mation pressure of about 0.7 Mb at ambient 


small. and. the total mass of material that is 
vaporized generally is small. 


SHOCK WAVE PROPAG TION 
AND DAMAGE PREDICTIONS 


density. corresponding to sublimation energy of. : 
about 7.900 cal gm (Table 9-17), This point The sequence of events for the gen- 


shown in Figure 9-40, labeled EF, at about 3.000 
cal gm. Table 9-2Q-indicates that the pressures 
assOcluted with vaporization of metals at 
ambient density are with some exceptions about 
} Mb. A survey of more than 30 common metal 
elements indicates that an average of 1 Mb for 
vaporization is a good approximation, especially 
if the Griineisen vale for Uie material is uneer- 
tain. Since a bur corresponds to 14.7 psi a Mb is 
the enormous pressure of about 1.45 x 10° psi. 
Thus. tremendous forces are involved in the 
pressure gradients associated with meta} vapori- 
zation at ambient density. Table 9-17 shows that 
vaporization usually involves several thousand 
calories per gram of energy. High explosive ma- _ 
terials (TNT. etc.) release about 1.000 cal/gm. 
Therefore. on a mass basis there is more energy 
associated with metal vaporization than with 
high explosives. Generally. the thicknesses of 
material evaporated by X-ray absorption is 


eration and propagation of a stress wave through 
the thickness of an aerospace shell and the dam- 
age produced is illustrated in Figure 9-4]. Cold - 
X-rays are deposited primarily in a relatively 
thin sheet of material at the front surface (Fig- 
ure 9-4]a). After the energy is deposited a com- 
pression wave propagates inward from the front 
surface, followed by a rarefaction that.causes 
the vapor and liquid to blow off (Figure 9-41). 
This rarefaction also may cause a spall of solid 
material from the front surface (Figure 9-41c¢). 


. Later the compression wave reflects from the 


back surface and returns as a rarefaction wave. 
This rarefaction wave. or the coincidence of this 
wave with the rearward moving rarefaction may 
cause the rear surface to spall (Figure 9-41d). or 
may cause fracturing or debonding. This process 
occurs within the order of a microsecond and 
generally is complete before the overall struc- 
tural motion occurs. The shock effects are 
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Figure 9-41. a Sequence of Spallation Following 
Radiation Deposition 
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extremely local. depending on the sheet thick- 
ness and not on the overall geometry of the 
structure. so spall damage. including fracturing 
and debonding. usually cannot be scaled accord- 
ing to the laws that govern structural behavior. 
Since the spall damage occurs early, possibly 


weakening the structure. it can have a strong 


effect on the subsequent response to the struc- 
ty the impulse that has been imparted. 


Propagation of stress waves and the dam- 


age predictions are not amenable to hand calcu- 
jations. Most stress wave propagation calcula- 
tions are performed on computers using one 
dimensional (1-D) models to represent a cross 
section of an aérospace vehicle shell. These stress 
wave propagation calculations are used routinely 
in design analysis and are the direct outputs of 
PUFF-like codes or their variations. The stress 
wave propagation predictions are used to: 

1. Obtain a calculated pressure time history 
at the rear surface of a 1-D sample for direct 
comparison to experimental values that are ob- 
tained by laboratory simulation or by under- 
ground nuclear tests. The agreement or lack of 
correlation serve as an important criteria for 
judging the adequacy of analytical tecliniques. 
The conclusions are almost universal that the 
largest sources of error in stress wave propaga- 
tion codes are associated with the mechanical 
aspects of the material behavior. such as failure 
criteria. and the elastic-plastic models. which 
should include strain rate and temperature 
dependence. 

2. Obtain maximum tension as a function of 
depth in the materials to indicate regions of po- 
tential material failure. debonding. and spalling. 
Definitions that are required for the materials 
failure level under these dynamic tensile loads 
are obtained experimentally. 

3. Obtain input data for structural response 
analysis. , 

There also are some 2-D shock propaga- 
ron codes being used for specific calculations 


involving cylinders, nose tips. etc. Peak pressure 
predictions using pure hydrodynamic representa- 
tions of the material in code caiculations are 
generally high by a factor of 2 to 4. Ifa dynamic 
’ elastic-plastic constitutive relationship is includ- 
ed. the predictions are improved significantly. 


9-38 Through-the-Thickness Elastic- 
Plastic Shock Propagation ie 


Considerable progress has been made in 
correlating measured and calculated through-the- 
thickness shock propagation in one dimension 
by using the elastic-plastic materials description 
in PUFF-type codes. The stress strain diagram 
shown in Figure 9-42 illustrates a two-wave 
nature of the elastic-plastic shock wave propaga- 
tion. The total stress. o, is the sum of the hydroa- 
‘static pressure, P, and an elastic stress offset of 
2/3 Y,. where Y, is the yield strength of the 
material (its compressive yield strength if that is 
different from its yield strength in simple ten- 
sion). The effect of the yield stress path OAB is | 
to propagate two stress waves having different 
velocities into the undisturbed material. The 
slope of the elastic portion. OA. is larger than 
the plastic portion, AB. The propagation speed 
is larger if the slope of the stress strain curve is 

larger. Hence, an. elastic precursor shock travel- 
ing with the elastic velocity, c, , runs in front of | 
the slower total stress wave which propagates 
with a bulk sound speed, c, 


where K is the bulk modulus and S* is the shear 
modulus. 


| An elastic rarefaction stress (4/3 Y, or 
Wice the stress offset) propagates from the rear 


of the shock wave at velocity Cy and overtakes 
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Figure 9-42, 


the slower plastic wave (path BC). This rarefac- 
tion wave causes a greater attenuation of the 
main stress pressure wave than is obtained in 
pure hydrodynamic calculations (up and down 
the Hugoniot). This greater attenuation of the 
elastic plastic wave propagation compared to 
purely hydrodynamic calculations results in 
better correlation with measured values of stress. 
The plastic unloading to a stress-free condition. 
path CD in Figure 9-42, also propagates with 
velocity c, and completes the cycle. 

Although the elastic-plastic code calcula- 
10ns improved the correlation of measured and 
computed through-the-thickness © stresses for 
some miaterials, if a material has ‘strong strain 
rate or temperature dependence these features 
also must be incorporated into the calculations. 
Stress wave propagation through materials such 
as foams require an entirely different constitu- 
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Material Description 


tive expression. and anisotropic materials such as 
three-dimensional quartz phenolic, can have 
quite different material properties in different 
directions. Similarly, propagation characteristics 
of inhomogeneous materials or materials that 
contain small radiation absorbers lead to compli- 
cated propagation analysis. Propagation in these 
materials has not been treated in an analytically 
satisfactory manner. Most of the current ana- 


lytical techniques used in the solution of stress — 


wave propagation in the materials are deficient 
in varving degrees in the following wavs. 

1. There is incomplete understanding of the 
material behavior under the types of conditions 
that result from X-ray energy deposition, so the 
material representations are inadequate in most 
cases. 

2. Experimental data for verification and 
modification of the equations of state and con- 


stitutive relations in the analytical techniques 
have been limited and of poor quality. 


9-39 Through-the-Thickness 
Material Response 


Simulation tests. and underground nu- 
clear weapons tests can identify failure and dam- 
age modes that result from shock wave propaga- 
tion and through-the-thickness material response 
on samples that are meaningful in terms of the 
expected response of reentry vehicles or aero- 
space shells. The experimental results are most 
meaningful when appropriate analytical tech- 
niques are used to. correlate the experimental 
results with the stress wave characteristics pro- 
ducing the respons: modes. Hundreds of com- 
ple\ aerospace structure. composites as well as 
samples of metals and plastics have been ex- 
posed during underground nuclear weapons tests 
at various X-ray fluence Jevels with different 
X-ray spectra and deposit } 


-RESPONSE ANALYSIS 


IMPULSE AND STRUCTURAL 


The vapor and/or liquid blowoff follow- 
ing A-ray deposition imparts a reactive impulse 


"to the aerospace structure. The structural re- 


sponse takes place for the most part after the 
shock wave discussed in the last section has dissi- 
pated. This shock wave may. have damaged or 
weakened the ‘structure by spallation. debond- 
ing. or fracture. The damage to the structure 
should be considered in the structural response 
analysis if an accurate prediction is to be made. 
Obviously. this is a complicated problem and is 
not amenable to hand calculation. However. 
some rough estimates of impulse can. be made 
and serve to illuminate the processes that take 
lace. 
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9-40 Impulse Calculations 
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Problem 9-5. Calculation of Vapor Blowoff impulse 


az The information presented in the pre- 
ceding paragraphs provides a means to obtain an “ 
approximation of blowoff impulse for certain 
materials. The following example is presented to 
illustrate the procedures. In practice, these cal- 
culations are performed with computer codes. 


through 9-37, and also Figures 9-35 
through 9-38, 9-43 and 9-44, and Table 9-17. 


on 


Ne nee OI oes nae nem nat cewele ts wee 


Problem 9-6. Calculation of the Location of Phase Change Boundaries 


The preceding paragraphs provide the 
information necessary to calculate the location 
of the various phase change boundary layers. 
The procedure is illustrated in the following 


8 Ups 


ney 9.37, and 9-40. 0. Sa sO Table O17 oa : 
Figure eae: 
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2 erent 


Response Analysis 


Most of the structure response analyses 
of reentry vehicies have been developed for sur- 
vivability. Since low level failure criteria re- 
sponses have been selected, only elastic deforma- 


« tion has been considered. Such survivability . 
"analyses are probably inadequate for the deter- _ 


mination of inelastic deformation failure modes 
and the levels required to define lethality. Essen- 


tially all structural response calculations require - 


data on: (1) structural material properties such 
as stress and strain relations as a function of 
strain rate and temperature, for yield, stress and 


‘strain, fracture points, elastic moduli? and so 
forth; (2) structural response forcing functions, - 


including asymmetric impulse loads, and in- 
deptif heating and thermal loads applied as dis- 
tributed and point loads; and (3) structural 
variation in shell wall thickness and material 
densities in various lavers. 


Bek REENTRY VEHICLE 
HARDENING 


fe 


X-ray hardening refers to techniques 


Ec 


4 


S reduce the susceptibility of reentry vehicles - 


to damage from X-ra 


it is beyond the scope of this manual to + 
provide a complete discussion of reentry vehicle ~ 
hardening techniques. Therefore, a brief general -” 
discussion of the techniques, with a few specific ~ 


examples, is presented in the following para- _ 


graphs. 


9-42 Balanced Hardening with Respect 
to Neutrons and X-rays 
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SECTION VI 


NUCLEAR RADIATION 
SHIELDING 


Ee Any mass of material between a nuclear 
radiation source and personne! or sensitive 
equipment will reduce the dose to the personnel 
or equipment compared to the free field dose at 
the same location relative to the source. The 
dose received by a person behind a building, in a 
building, in a field forti ion, in a tank, or in a 


~ ship will be Jess than that which would be re- 
ceived in an exposed free field position. 

In principle. it would be possible to cal- 
culate the dose behind a shield by using simpli- 
fied neutron and gamma ray transport equations 
similar to those given in Chapter 5: however. 

_ several considerations make hand calculations 
impractical. 

@ The neutrons and gamma rays will arrive at 
the shield from manv directions after hav- 
ing been scattered during their transport 
through the atmosphere. 

@ Calculation of radiation transport from a 
point source through single geometries can 
be done with reasonable accuracy by using 
ui exponential factor and a buildup factor 
similar to those described in’ Section I, 
Chapter 3 for transport through homo- 
geneous air. but the geometries of shielding 
enclosures are seldom simple. and. as men- 
tioned above. the radiation incident on the 
shield is not unidirectional... This has the 
effect of making the problem similar to a 
multi-source problem. 

@ The attenuation of neutrons by inelastic 
scatter or radiative capture results in the 
emission of one or more gamma rays, 
which can contribute to the dose within 
the shielded area even though the neutrons 
are removed from the penetrating radiation 
beam. 

@ The specific location of the burst and the 
receiver can affect the degree of attenua- 
tion, e.g.. different members of a tank crew 
will have different degrees of protection 
from the same explosion, and any par- 
ticular member of the crew will have a dif- 
ferent degree of protection from initial 
radiation from bursts at different locations 
relative to the tank. In other words, each lo- 
cation within each shielded structure must 
be evaluated separately. or at least a suf- 
ficient number of locations must be evalu- 


ated to determine the overall shielding el- 
ficiency of the structure or piece of equip- 
ment.. 


’. These complexities, and others, make it neces- 


sary to use computer codes to determine the 
shielding efficiency of any structure or piece of 
equipment. Even when using the codes. some 
simplification in the description of the shield 
geometry usually is necessary. However. 
satisfactory methods for treating fairly complex 
geometries have been developed. These are de- 
scribed in DASA 1892, “Weapons Radiation 
Shielding Handbook” (see bibliography). Results 
of some sample calculations are also provided. 
Unfortunately, sufficient calculations have not 
been performed to permit the development of a 
generalized simple method for shielding calcula- 
tions that would be suitable for inclusion in this. 
manual. 


Table 9-22 provides some estimates of 
the shielding afforded by various structures and 


' types of military material. These estimates are 


given in terms of a “dose transmission factor.” 
which is defined as the ratio of the dose mea- 
sured behind the shield to the dose that would 
be measured in the absence of the shield. Some 
of the transmission factors shown in Table 9-22 
were obtained by measurements at weapon tests 
or are extrapolations from such measurements. 
Other transmission factors were obtained by 
relatively detailed calculations. while still others 
are mere estimates. Ranges of values are given 
for the dose transmission factors for most struc- 
tures. These ranges result from two causes: un- 
certainty in the estimates themselves, and varia- 
tions in the degree of shielding that may be ob- 
tained at different locations within a structure. 
Separate transmission factors are given for initial 
gamma rays and residual (fallout) gamma rays 
for two reasons: the average energy of the initial 
gamma rays is higher than the average energy of 
the residual gamma rays, and there are different 
source geometries, 
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Table 9-22. @ Estimated Dose Transmission Factors (interior Dose/Exterior Dove 
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y No specific reliability can be attached to 
the dose transmission factors shown in Table 
9-22. They are felt to be reasonable values for 
generalized applications: however, specific 
shielding problems should be addressed by 
methods such as those described in DASA 1892. 


SECTION Vil 
=. TREE — COMPONENT PART 
‘D CIRCUIT RESPONSE 

J This section contains a discussion of the 
effects of nuclear radiation on electronic compo- 
nent parts and circuits. The intent is to provide 
an introduction to the failure modes and re- 
sponses that result from ionization. displace- 
ment. and thermomechanical shock. A more 
detailed description of the effects is contained in 
the TREE Handbook (see bibliography). 

Since semiconductor devices generally 
are Considered to be the most sensitive devices in 
an electronic circuit. emphasis is placed on 
them. However, the discussion does include 
other electronic component parts. 
ena Within the scope of this manual, this sec- 
1on can only aid in establishing levels where 
electronic equipment might experience failures 
as a result of the radiation environment and indi- 
cate the extent of the problems. The numbers 
-presented are only estimates for a generic group- 
ing or subgrouping. It must: be understood that 
it is very difficult to give accurate radiation 
effects data for generic types of electronics. For 
example. transistors can show sufficient amplifi- 
cation loss to be useless after receiving.10!° to 
10!5 n/em? (E > 10 keV, fission). That is 5 
orders of magnitude difference and is relatively 
useless information. Fortunately, transistors can 
be subgrouped to improve their characterization. 
The problem of accurate prediction, however, is 
not completely solved. Even if one particular 
type of transistor is considered, the neutron 
fluence required to reduce the amplification to a 


particular value can vary more than an order of 
magnitude. Additionally, the failure threshold 
for a transistor in one circuit may be greatly 
different from the failure threshold for that 
same transistor in a different circuit. i-e., the 
response of each component part and the circuit 
configuration must be known to make a reason- 
able estimate of circuit response. 

The transient radiation effects are 
classified into three groups: transient. perma- 
nent, and .thermomechanical. The transient 
effects result from nonequilibrium free-charge 
conditions introduced through the ionization 
phenomena (see Chapter 6). The permanent ef- 
fects are attributed to physical property changes 
of the irradiated materials caused by energetic 
particles (neutrons and secondary electrons). 


The thermomechanical effects result primarily 
from the absorption of low energy photons. 

The class of effect is defined by the 
primary effect induced by the radiation. For 
example. the short term currents that result 
from ionization may trigger a digital state 
change or may permanently damage a device in 
some other manner. This phenomenon is treated 
as a transient effect, because the permanent con- 
sequences result from circuit action and/or de- 
vice construction rather than from direct radia- 
tign induced material property change. 

The terms “radiation hardened’’ or 
“radiation resistance” frequently are connotated 
incorrectly. These terms mean that someone or 
some group has examined the component part, 
circuit, or system, made some modification 
and/or used part selection criteria and has deter- 
mined that the electronics will withstand certain 
levels of certain types of radiation. The 
experience of the person making the examina- 
tion as well as the resulting modifications can 
both vary widely. Consequently, when these 
terms are used they should be accompanied by 
the answers to the following questions. 


@ Hardened to what type of radiation? 
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@ Hardened to what levels of radiation? 
@ What procedure was used for hardening? 
@ How has the hardening been verified? 


Based on the answers to these questions, the 
applicability of that “hardened” component 
part, circuit, or system for his intended use can 
be judged. 


SEMICONDUCTOR COMPONENT 
PARTS ‘ 


9-44 Transient Effects i 


The transient effects observed in semi- 
conductor devices exposed to nuclear radiation | 


result from the creation of excess charge carriers 
(ionization) that cause current and voltage 
changes. These changes do not cause permanent 
damage directly to the semiconductor material. 
However, their presence may produce perma- 
nent changes as a result of current overloads to 
some components, loss of information stored in 
memory units, or by the creation of premature 
signals. Furthermore, such transients can cause 
saturation of some circuits for times that are 
long compared to the duration of the ionizing 
radiation pulse, and may thus cause system fail- 
ure during. a significant and possibly critical 


pepiod. 
a Most of the ionization produced in a 


nuclear weapon environment is due to the 
photons and 14-MeV-neutrons, although all nu- 
clear radiation can contribute. Because of this 
predominance of photon excitation of carriers, 


the observed currents and voltages are often- 


referred to as photocurrents and photovoltages, 
even though the prefix may. not always be 
proper. The mechanisms whereby this ionization 
occurs are described in Chapter 6. 

Before discussing the transient effects in 
particular semiconductor devices, it is desirable 
to review briefly the important basic physical 


processes responsible for device behavior. These 
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processes are common to all semiconductor de- 
vices and must be reasonably well understood 
before the description of the device response ca 
llowed intelligently. 
Current is the manifestation of charge 
carrier movement within a material. As. de- 
scribed in paragraph 6-4, the movement may be 
random, random-with a diffusion from regions 
of high concentration to regions of low concen- 
tration, or, if an electric field is present, the 
movement may be a drift in the field super- 
imposed on the normal random scattering. 
Carriers may. be trapped by impurities, which are 
always present in solid state devices. Eventually, 
the trapped carriers may be annihilated by their 
mates (oppositely charged carriers) in a process — 
called recombination. The net result of these 
processes is that the free carriers diffuse and/or 
drift until they are trapped and are then prompt- 
ly recombined. Therefore, any electrical currents 
in the materials are considered to have two com- 
ponents — the drift component and the diffu- 
sion component. The time during which a charge 
carrier is free to move by drift and/or diffusion 
is called the lifetime of the charge carrier. 

The simplest semiconductor device is the 
diode. The diode is made up of pure semicon- 
ductor material doped with impurities so as to 
be considered to have two regions. Using pure 
silicon as the example, its atomic structure is 
considered to have four electrons in the outer 
orbit. In its crystal form, these outer orbiting 
electrons are considered to be shared with 
neighboring silicon atoms (see Figure 9-45). 
When the pure silicon is doped with an impurity, 
the impurity atom replaces one of the silicon 
atoms in the crystal lattice. The impurities are 
selected to have either three or five outer orbit- 
ing electrons. Silicon doped with an impurity 
having five outer orbiting electrons will have 
extra electrons that are more or less free to 
move around, and hence are available for con- 
duction. Silicon doped with impurities with only 
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three orbiting electrons does not have sufficient 
electrons to share with its neighboring silicon 
atoms. This lack of an electron is called a hole. 
For conceptual purposes, the hole can be treated 
like an electron with a positive charge. Conse- 
quently, in semiconductor discussion the term 
charge carriers refers to the electrons and holes 
made available for the conduction -process by 
the impurity doping. 

ee Figure 9-46 illustrates a semiconductor 
junction (diode), with the two impurity regions 
shown separated by a solid line. This represents 
what conceptually would be the junction of the 
two materials. However, in reality the silicon is 
all one crystal. and there is really no physical 
junction. There are. as illustrated. regions. at 
both ends of the crystal which are predominant- 
ly either PP type (excess of free holes) or Ntype 
(excess of free electrons). A transition region, 
called the depletion region. exists in the center. 
Within this region. holes from the P region com- 
bine with equal numbers of electrons from the N 
region. As a result. only a few free charge car- 
riers remain in the region when equilibrium has 
been reached. Depending on the number of 
charges removed from each region, a voltage 
(electric field) will be developed across the 
depletion region. The voltage across the deple- 
tion region under equilibrium is such that any 


holes (positive charge carriers) introduced into 


the region would migrate by the drift process to. 


the P region, and a negative charge carrier (elec- 


- tron) introduced into the depletion region 


would migrate by the drift process to the N 
region. Under normal conditions conventional 


current flow* would allow current to flow for | 


ideal diodes only from the P side to the N side. 
If, however, free carriers were generated in the 
depletion region, conventional current flow 
would dictate that a current, proportional to the 
number of free charge carriers generated, would 
flow from the N region to the P region. This is in 
the reverse direction from the normal conven- 


tional current flow through an ideal diode. This — 


reverse current, if generated in the depletion 
region by ionizing radiation, is called the drift 
component of the photocurrent. in a PN junc- 
tion. The diffusion component of the photocur- 
rent is the result of free carrier generation by the 
ionizing radiation in the P and N regions near 
the depletion region. Those hole-electron pairs 
generated far from the depletion region will be 
trapped and recombined before they can be- 
come effective. The effect of the radiation is to 


5 in this chapter, normal current flow is considered to be 
cOfventional current flow. This is opposite to actual electron 


flow. 


Where junction would be 
if it were a line 


Excess of 
free holes 


Excess of free 
electrons 


pg alee ere this region constitutes the junction 
and is called the depletion region. 


Figure 9-46. & Ulustration of a Semiconductor Junction S&S 
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generate a current in the. reverse direction 
through the diode junction. which tends to for- 
ward bias the diode. An expression that general- 
ly will predict the photocurrent for many cases 
of interest is: ° 


pp max 


I = qK,DA (w+ Ll). 


where 


= the maximum* pri- 
mary photocurrent for 
a PN junction 


pp max 


gq = electronic charge 1.60 
x 10°29 coulombs,_ 


A. = the energy-dependent 
free-charge-carrier- 
generation constant 
(electron-hole pairs’ 
em* * roentgen) 

Kk. (for silicon) = 4.2 x 10!4  electron- 
hole pairs/cm* + R 
K, (for germanium) = 1 x 10! 4 electron-hole 
pairs‘em? + R 

D = the gamma exposure 
rate. R/sec (The gam- 
ma ray. photons have 
energies which would 
give best results with 
the A, constants given 
above.) 


A = junction area in cm- 
(A can typically vary 
from 0.3 x 1074 to 0.2 
x 104 em?) ~ 

we = depletion laver width 
in cm (w varies with 
applied voltage and is 
best expressed as w = 
wy, (Yo - vy? where 
V, is the built-in junc- 


tion potential (~0.35 
volts for germanium 
and ~0.72 volts for 
silicon) w, typically 
varies from 0.5 x 104 
em to 1 x 10° cm, Vy 
is the junction voltage. 
and b typically vanes 
from 0.05 to.0.5. 


L = diffusion length for 
minority carriers on 
the side of the junc- 
tion with the longer 
diffusion length in cm 
(L can vary typically - 
from 0.15 x 10! cm 
to 0.2 x 10% cm). 


The duration. of the photocurrent. de- 
pends on the time required for free holes. and 
free electrons to be trapped and to recombine. 


. Since these times are short. the photocurrent is a 


pulse of current of relatively short duration. 
Typical pulse shapes are shown in Figure 9-47. 


Tunnel diodes are at least an order of 
magnitude more radiation resistant than diodes | 
in general, since they are characterized by small 
geometry, heavily doped PN regions, and narrow 
junctions. The short-circuit photocurrent can be 
estimated from the equations applicable to 
general diodes. Since the lifetime of. minority 
carriers in tunnel diodes is very short, the re- 
sponse of a tunnel. diode to a pulse of radiation 


_ iS expected to follow the radiation pulse. Nor- 


mally, the width of the depletion region for 
tunnel] diodes is so narrow that the photocurrent 
will consist primarily of the diffusion compo- 
nent. The maximum photocurrent under short 


bg For those more theoretically inclined this equation applies 
in I case where T < t,/4 and is the steady-state value, A more 


detaiied equation and development is the TREE Handbook, 
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Relative Primary Photocurrent 
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Figure 9-47, & 


circuit conditions can thus be estimated as 
follows: 


io 5k eGR DAU set) 


pp max 


where 


q, K,. BD and A are aehaed as in the diode 
expression, and L, and L, are the diffusion 
Jengths in the N and P regions. respectively. 
Diffusion lengths are comparable on each side 
of the junction. . 


The radiation induced offset voltage in 
an open circuited tunnel diode is influenced by 
the tunneling current at low dose rates, but at 
higher dose rates a tunnel diode behaves like a 
conventional diode. This behavior occurs. be- 


cause very little induced voltage is required to © 


set up a tunneling current equal to the diffusion 
current of excess carriers. When the tunneling 
current reaches a maximum, the induced voltage 
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Relative Shapes of Diffusion Component 
of Primary Photocurrent 


rises sharply to the value for a conventional 
diode 


A transistor is 4.more complicated de- 
vice than a diode, both structurally and func- 
tionally. The main functional advantage of the 
transistor is the fact the transistor has the ability 
to amplify its input signal. The transistor has 
three impurity regions that alternate, ie., NPN 
or PNP. When an operating transistor* is ex- 
posed to transient ionizing radiation, a current 
pulse is observed in the external circuit. This 
current pulse, which may be orders of magni- 
tude larger than that of a diode with comparable 
dimensions, can reach a peak value at a time 
later than the radiation peak, and can, in some 
cases, continue for milliseconds. 


The following paragraphs are concemed with bipolar tran- 
aistors; not to MOS, junction field-effect transistors, or thin-film 
transistors. ; 


Fo] This characteristic behavior of transis- 
tors is the result of the amplifying properties 
acting on the primary radiation induced photo- 
currents. The electrical action of the transistor 
creates a secondary photocurrent that is typi- 
cally. greater than the primary photocurrent by a 
factor equal to the gain (a measure of ability of 
the transistor to.amplify) of the transistor. Anal- 
ysis of the radiation response of a transistor 
involves the determination of the primary 
photocurrent. followed by a calculation of the 
magnitude and duration of secondary photocur- 
rents under the given circuit conditions. 

In transistors. the primary photocurrent 
is generated in five regions: in the collector- and 
emMitter-junction depletion regions. in the base. 


and in the emitter and collector bodies lying . 
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Depletion regions 


within a diffusion length of the junctions. In 
most cases. the generation of primary photocur- 
rent in the emitter body and within the junc- 


tions can be neglected since the emitter body 


and the junction volumes are a relatively small 


. part of the total generating volume (see Figure 


9-48). 
The secondary photocurrent is produced 
y the accumulation of excess majority carriers 
in the base region as a result of the flow of 
primary photocurrents across the PN junctions 
of the device. This excess charge, which is con- 
fined in the base by the built-in junction fields. 
is the correct polarity to forward bias the emit- 
ter base junction and to cause normal current to. 
flow. The collector current that is produced is 
called secondary photocurrent. This collector 


Collector 


The farger junction area results in a 
larger photocurrent 


NPN Transistor Hiustrating the 


Primary Photocurrent 
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current continues to flow until the excess charge 
stored in the base can either recombine with 
minority carriers or flow out through the external 
b ead. 


The magnitude of the aolicerar current 


pulse will depend on the dose rate if the radia: 
tion pulse is long compared to the lifetimes of 
the transistor base and collector, but it will 
depend on the total dose if the pulse is shorter 
than those lifetimes. This results from the charge 
transferred from the collector by the primary 
photocurrent being stored in the base region. 
Since many transistors have lifetimes that areas 
long aS or longer than typical nuclear weapon 
radiation pulses, the prompt dose is quite often 
the most important factor. 


A simplified linear approximation of the. 


primary photocurrent in a transistor can be esti- 
mated with the following equation: 


D 


Ip ” CONST. UP 0.0), 


D = the gamma ray exposure rate, R/s 


CONST. =.0.83 x 10® for NPN devices and 
0.42 x 108 for PNP devices* 


t. = charge storage time (in ys) for the 
device of interest (1, typically can 
vary from a-few nanoseconds to 
hundreds of microseconds). 


The duration of the photocurrent re- 
sponse depends on the radiation pulse width, the 
radiation storage time (which is different than 
the charge storage time (t,)), and the time for 
holes and electrons to recombine. The radiation 
pulse width is determined by the weapon. The 
time for holes and electrons to recombine is 
iNustrated for a diode in Figure 9-47, which 
applies generally for transistors. The radiation 
storage time, f,,, is defined by Figure 9-49a. ‘The 
radiation storage times for several dose rates are 
shown in Figure 9-49b. These latter curves are 
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Figure 9-49, & Radiation Storage Time | 


@:.- equations only apply for the steady-state estimate of 
the primary photocurrent for silicon pianar and mesa transistors 
with rated maximum continuous collector dissipation below 0.8 
watt at 250C, More detailed information about prediction is 
available in the TREE Handbook (see bibliography). 
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for a particular bias and radiation condition dur- 
ing storage. They are presented as being illus- 
trative. and no attempt should be made to gen- 

eralize from the curves. It is necessary to analyze 
a transistor or diode in the particular circuit con- 
figuration in which it is used to determine the 
threshold for circuit malfunction. As a result of 

“the many different circuit configurations and 
bias conditions that can occur for transistors and 
diodes, failure thresholds. that would be of value 
can not be specified generally for the TREE 
environment. 


A silicon controlled rectifier (SCR) is a 
solid state semiconductor device composed of 
four layers of alternate-impurity semiconductor 
material containing three PN junctions. The SCR 
iS an active switching element that will remain in 
a nonconducting or “off” state until turned on 
or “fired” by a low-level control signal on the 


gate. Jt wil] then remain ‘‘on” without need for. 


a sustaining control signal. The SCR is turned 
“off by reducing its ouiput (anode) current 
below the “dropout” level. Radiation induced 
currents. like those discussed for diodes and 


transistors. are a direct function of the junction ~ 


areas. diffusion lengths. etc.. and thus are diffi- 
cult. to predict since values for the parameters 
usually are not available. Since these currents. 
above a threshold. can induce changes in the 
state of an SCR. some method is required to 
predict the magnitude of the radiation induced 
currents. or, more specifically. the radiation 
threshold above which switching occurs. 


& It has been found that the transient - 
radiation switching thresholds (critical radiation 


exposure rate) for SCR’s are functions of the 
radiation pulse width. The exposure rate re- 
quired to trigger an SCR becomes constant for 
pulse widths longer than a critical value. This 
critical value is a function of the device minority 
carrier lifetime and the device delay “turn on” 
time. For pulse widths less than the critical 
value. the exposure rate required to trigger an 


SCR increases rapidly as the pulse width ap- 
proaches zero. The dependence of the switching 
threshold of a 3A60A SCR on pulse widths and 
gamma ray exposure rate is shown in Figure 


“9-50. The critical pulse width for this device is 


a ximately 2 microseconds. 

Typically the SCR type of device would 
not be expected to fire below 10° rads (Si)/sec. 
For most cases. however. the pulse width is suf- 
ficiently short that the devices ate dose de- 
pendent. Failures occur typically at prompt 
doses between 0.1 and 1 rad (Si). 

Field-effect transistors (FET) are a 
famuy of unipolar devices that have pentode-like 
characteristics. The three major categories with- 
in this family are the junction FET. the metal- 
oXide insulated-gate FET (MOSFET). and the 
thin-film insulated-gate FET (TFT). The geom- 
etry and construction features of typical field- 
effect transistors are shown in Figure 9-51, The 
basic structure of the FET devices involves a 
source. a gate. and a drain in rough functional 
correspondence to the familiar cathode, grid, 


‘and plate of vacuum tube technology. 


The mechanisms by which radiation gen- 
er. photocurrents in an FET are not substan- 
tially different from those for bipolar transistors 
and diodes. The important radiation parameters 
in an FET are the transient gate and drain-to- 
source currents. Possible sources of transient 
currents in FET’s can be grouped into the fol- 
lowing categories: . 

@ Leakage currents across PN junctions that 
behave like PN junction photocurrents dis- 
cussed previously. 

@ Direct modulation of the channel conduc- 
tivity and mobility (usually. applicable at 
high (>10® rads (Si)/s) dose rates). 

@ Leakage currents through the gate oxide 
layer (applicable to the metal oxide and | 
thin-film FET’s).. 

© Secondary emission (see Chapter 6) and 
atmosphere ionization currents. 
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Figure 9-51. & Construction of Typical Junction, MOS, and 
Thin-Film Field-Effect. Transistors 
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a Significant problems generally are not 
caused in field-effect transistors at dose rates be- 
low 10® rads (Si)/sec. 


9-45 Permanent Effects BD 


Permanent effects in semiconductor de- 
vices are those that can be attributed to physical 
property changes that result from the direct 
interaction of the radiation with the material of 
interest. These property changes typically last 
for periods that are long with respect to the re- 
covery times of the components. These property 
changes occur in a very short time period and 
result in a rapid change in the operating charac- 
teristics of the device. A closely related effect is 
called rapid annealing. which is the process by. 
which an initially large change in device param- 
eters recovers very rapidly, approaching the 
smaller change observed several minutes after 
the radiation exposure. 

Most permanent effects in semicon- 
ducior devices subjected to a nuclear radiation 
pulse result from damage to the semiconductor 
material by energetic neutrons (E > 10 keV): 
however, the effects of gamma rays and sec- 
ondary electrons must not be underestimated. In 
certain devices such as MOS field-effect transis- 
tors the effects of ionizing radiation can be the 
Prigcipal causes of permanent failure. 

Permanent effects can be grouped into 
two categories — bulk, and surface effects. Bulk 
effects are changes in the device characteristics 
that can result from damage to the bulk ma¢ 
terial. Surface effects are changes that are gen- 
erally caused by radiation induced ionization 
near the surface of the device. Bulk damage ef- 
fects from neutron radiation usually can be 
predicted within a factor of 2, while surface 
effects are generally unpredictable. 

Bulk effects result from electron, gamma 
ray, and neutron induced lattice displacements 
in the bulk of the material (see Chapter 6). Fast 


neutrons lose energy primarily by elastic colli- 
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sions with the semiconductor atoms and cause 
large disordered clusters to be formed within the 
material. Gamma radiation loses energy by creat- 
ing Compton electrons (see Section I, Chapter 
5), which may cause lattice displacements. Since 
electrons have such a small mass, they primarily 
cause vacancy-interstitial pairs rather than clus- 
ters of defects that are typical of neutron dam- 


age. Lattice damage that results from gamma 


radiation usually is of secondary importance, 
unless a large gamma dose (>10° rads) is 
absorbed by the material. 


Lattice damage: degrades the electrical 


characteristics of semiconductor devices by in- 


_ creasing the number of trapping, srattenns, and 


recombination centers. 
@ The trapping centers remove carniers from 
the conduction process. 


@ The additional scattering centers reduce the 
mean free path of the free carriers. Since 
the mobility is directly proportional to the 
mean. free path, radiation exposure reduces 
the mobility of charge carriers. 

@ The recombination centers decrease the 
minority carrier lifetime according to the 
relationship (see Chapter 6). 


af 
T, 


—+ 
7 Ky, 


7, = minority-carrier lifetime at fluence y 
in seconds, 


_ lifetime), in seconds, 
K = lifetime damage constant, cm?/ 
(neutron-seconds), 
y = fast neutron fluence, neutrons/cm2. 


Permanent effects also can be caused by 
radiation induced changes in the semiconductor 
surface. The changes in the surface conditions 


Tt, = initial minority-carrier lifetime (bulk . 


attributed to radiation that can cause permanent 
effects are surface charging mechanisms and 
changes in the surface recombination velocity. 
The most likely charging mechanisms are the 
collection of ions from a gas in the atmosphere 
surrounding the semiconductor device and the 
ejection of electrons from dielectric materials 
which are either deposited on the semiconductor 
surface or in which the device is encapsulated. 
These ions migrate under the influence of elec- 
tric fields. As a result of the collection of these 


charges. inversion layers can form near. the sur- — 


face. causing large increases in leakage currents. 
In silicon devices these leakage currents result 
from recombination-generation in the enlarged 
depletion region. 

lonizing radiation can cause changes in 


recombination. velocity. which has deleterious 


effects on the effective lifetime according to the 
relation 


Forward voltage 


for constant current 


Leakage current 


or reverse Current 


voltage 


Current, | 


Reverse breakdown 


] 1 ] 
—~ = —+t —, 
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Poa = effective lifetime, 
bulk lifetime, 


= surface lifetime (an inverse function 
of recombination velocity). 
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Surface effects usually are negligible compared 
to bulk lifetime damage for most conventional 
devices in a transient radiation environment. 
Field-effect devices, where radiation induced 
surface changes are the primary damage 


- mechanisms, represent an important exception. 


The general effects of nuclear radiation 

-on semiconductor diodes are summarized below 
and are illustrated in Figure 9-52 and 9-53:. 
@ The forward voltage of the diode at con- 
‘stant current normally will incréase as a 


Preirrodiation 


Postirradiation 


AV 
Slope of curve Ar 


is forward resistance 


Voltage V 


—— 


Figure 9-52. | An Illustrative Diode Characteristic for 
Pre- and Post-irradiation by Neutrons 6 
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Figure 9-53. & An Illustration of Rise Time and Storage Time | 


result of changes in resistivity and mobility 
in the bulk material. 


® The dynamic forward. resistance will in- 


crease as a result of changes in resistivity 
and conductivity modulation. 

@ The surface effects and increases in carrier 
generation in the space-charge region will 
cause the reverse current to increase. 


@ The reverse breakdown voltage normally 
will increase because of an increase in re- 
sistivity. 

@ The switching characteristic of ‘the diode 
generally will be changed. The rise time will 
increase, and the storage time will decrease 
as a result of lifetime damage. 

Diodes generally are an order of magni- 
tude more resistant to radiation than transistors 
of similar type. For this reason, theoretical and 
experimental studies have concentrated on tran- 
sistors rather than diodes. Prediction techniques 
for diodes are complicated by the fact that the 
various bulk and surface damage mechanisms 
interact in an intricate manner. The wide variety 
of diode types, i.e., material, doping level and 
profile, geometry, etc., all tend to make compre- 
hensive prediction schemes difficult and inaccur- 
ate. However, some quantitative trends in the 
changes in diode parameters can be given. The 
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more important changes in diode characteristics. 
from a circuits point of view, are the increase in 
forward voltage at constant current and the in- 
crease of reverse leakage current. Changes in the 


~ dynamic forward resistance, breakdown voltage. 


and switching characteristics usually are of 
secondary importance. 


The forward voltage at constant current 


generally starts to increase at a fluence of 10!3 
to 10'* n/em? (E > 10 keV, fission), though 
some diodes exhibit changes at fluences as low 
as 10!2 n/em? (E > 10 keV, fission) while 
others show no change up to fluences as large as 
10'5 n/em? (E > 10 keV, fission). Usually, flu- 
ences of about an. order of magnitude greater 
than the fluences required to cause the initial 


’ change will double the forward voltage. 


The reverse leakage current usually in- 
creases with exposure, but decreases also have 
been observed. Nominally, changes begin at flu- 
ences from 10'? to 10!4 n/em? (E > 10 keV, 


fission). Gamma ray doses as low as 5 x 10* rads 


(Si) have caused significant leakage currents. 
Germanium devices generally have larger changes 
in leakage current than silicon devices. 


The changes in breakdown voltage are 
typically the largest for diodes with high break- 
down voltage. Reference voltage diodes (zener 


diodes) are relatively resistant to radiation with 
reference voltage changes less than S percent at 
10}= nom? (E> 10 keV. fission). 

€ The storage time is directly proportional 
to lifetime. Hence. the fluence at which the stor- 
age time wil] be reduced to one-half the preir- 
radiation value will be 


where 
7, is the minority. carrier lifetime. 
K is the damage constant (see Chapter 6). 


Other diode types are specifically de- 
Signed for rectifier application where high break- 
down voltage and a low voltage drop are re- 
quired even at high current. These diodes usually 
‘are designed with @ PIN junction. which results 
in a device that is less sensitive to radiation than 
the standard tvpe diodes. 
Selenium rectifiers and hot carrier 
iodes (meta] semiconductor junction} appear to 
be more radiation. resistant than either ger- 
manium or silicon diodes because of their mia- 
terial and structural differences. Reactor tests 
conducted on the HPA-2300 series hot carrier 
diodes confirm the relative radiation hardness of 
these devices. Most units tested remained within 
manufacturer’s specifications at fluences of 3 x 
10'5 nicm? (F > 0.1 MeV. fission) and 9 x 105 
rads (Si). 

Diodes classified as “tunnel diodes” are 
easily recognized by their forward current char- 
acteristic, which shows a region of negative resis- 
tance. The effect of radiation on tunnel diodes is 
observable on the current-voltage (I-V) charac- 
teristic at fluences between 10!* and 10! 
n/em? (EF > 10 keV. fission), and can be sum- 
marized as follows: 

@ The slope of the primary tunneling current 
is not changed: however, the peak current 


is slightly reduced as a result of the redis- 
tribution of the electrons in the defect 
states. - 

'@ The valley current at a given voltage in- 
creases due to the additional tunneling via 
defect states. 


@ For voltages larger than valley. voltage. 
there is an increased current for a given 
voltage due to the excess current. which 
predominates over the diminishing norma] 
diode current. 

Figures 9-54 and 9-55 show representa- 
tive germanium and silicon tunnel diodes under 
neutron irradiation. The figures show that tun- 
nel diodes are still operational at S x JO!*_ 
niem? (E > 10 keV. fission). which indicates the 
relative radiation hardness of these devices. 

The general effects of nuclear radiation 
on bipolar transistors can be summarized as fol- 
lows: : : 
@ The current gain (amplification) of the 

transistor will be degraded as a result of 

lifetime damage to the bulk material. 

Degradation of gain will be greatest immed- 
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Figure 9-54, B Fast Neutron Bombardment Effect 
on Voltage-Current Characteristic 
of a Germanium Tunnel Diode 
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iately following a burst of nuclear radia- 
tion. and the gain will recover rapidly to a 
quasi steady-state value. The annealing can 
go on for weeks. but usuaily the current 
gain recovery is small or negligible after the 
imitial recovery. 

@ The reverse leakage current will increase as 
a result of surface effects and carrier gen- 
eration in the space-charge region. 

@ Changes oceur in the punch-through 
voltage. and the ‘base-to-emitter and 
collector-to-base breakdown voltage as a 
result of changes in resistivity. 

® Increases in base-spreading resistance, col- 
lector body resistance, and saturation volt- 
age result from changes in resistivity and 
conductivity modulation. 

@ The switching characteristics also are 
changed slightly — exemplified by de- 
creased storage time and increased turn-on 
time as a result of changes in lifetime and 


resistivity. 
@& Since transistors usually are the most 
vulnerable devices used in conventional circuits. 
predictions of circuit response under radiation 
conditions wil] be limited by the accuracy with 
which transistor behavior can be predicted. The 
reduction of current gain generally will limit the 
usefulness of the component before the other 
factors listed above become a serious problem. 
Therefore, emphasis is placed on the prediction 
of current gain degradation. However. in some 
applications saturation voltages and/or leakage 
currents across reverse biased junctions may be 
the limiting factors. 

The structure of a device is an icortant 
actor in determining its radiation resistance. A 
general rule is that the thinner base, higher fre- 
quency, and smaller junction area devices 
usually have better radiation resistance. For 
example. the diffused-junction devices usually 
offer a resistance to radiation about one order of 


magnitude better than that offered by the alloy- 


junction devices. 


Experimental data for conventional tran- 
sistors in the form of generalized gain degrada- 
tion curves are shown in Figure 9-56 for some 
common transistor types. The ratio 8 0B, repte- 
sents the ratio of gain at some fast neutron. flu- 
ence to the gain prior to irradiation. Figure 9-56 


‘is representative of preliminary data from 


steady-state reactor experiments. The transistor 


_types have been placed in their respective 


regions on the basis of where a majority of 
samples of a given type fell on the graph. Cau- 
tion should be exercised in using and interpret- 
ing the information in the figure, since irradia- 
tion temperature, irradiation source. measure- 
ment conditions. etc. are not specified. 

A further word of caution should be 
injected concerning the interpretation of gain 
degradation data. A sharp decrease in 6 occurs 
during exposure and then rapidly anneals to a 


final value that is commonly measured. The time 


dependence of gain degradation can best be 
interpreted with the use of the annealing factor, 
F, defined as follows: 


Lae 
Be By Kn) 
a 


BAe) 8, 


where 


6 .(t), A(t) are gain and damage constants as a 
function of time after a fast burst of nuclear 
radiation. 


By(%), K() are the steady-state values of gain. 
and damage constants. 


_ The magnitude and. form of F varies 
with temperature, injection level, doping level, 
and impurity content. Typical values for the 
annealing factor at 1 msec range from approxi- 
mately two to three for NPN transistors and 
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slightiv smaller tor PNP transistors. The anneal- 
ing factor can be as high as 7 for low injection 

Situations and immediately after turning on. a 

device that was off during the neutron pulse. 


Other transistor parameters also. can be 
egraded. Permanent increases occur in leakage 
currents across reverse-biased junctions. In gen- 
eral. the discussion of diode leakage is applicable 
lo transistor leakage. 
eS Changes in. breakdown voltages. punch- 
through voltage. and collector and emitter body 
resistances are negligible. at fluences where the 
gain is still usable. The effects of nuclear radia- 
tion on these parameters can be analyzed from 
comments made. about diodes. : 
The changes in saturation voltage and in 
le switching me as a result of nuclear radia- 
tion are of interest for switching applications. In 
" Many cases the saturation voltage may appear to 
increase at relatively low levels of radiation. but 
actually the transistor is losing base drive and is 
coming out of saturation. The only significant 


increases in saturation voltage are seen at high. 


collector currents. These increases. however. 
occur at high neutron fluences. : 

The switching time of a transistor is re- 
erred to as the turn-on time and the turn-off 
time. which consist of the delay time. 7,. the 
rise time, 7,. the storage time, 7,. and the fali 
time. 7,. First-order theory indicates that. with 
the exception of the rise time. these parameters 
either remain relatively constant or decrease 
with radiation. Normally. the decreases are 
larger than the increases in these parameters: 
thus a net reduction of transistor switching time 
occurs with radiation, which usually is desirable. 
The largest changes occur in the storage time. 
which is proportional to the lifetime. Thus, 


2 dbs, 
" Kt, 


is the fluence at which the storage time is re- 
duced by one half. 


_ The general effects of nuclear radiation 
on field-effect transistors can be summarized as 
follows: <P. 
~ @ Changes occur in the threshold voltage. }’;. 
These changes in threshold voltage affect 
most of the field-effect transistor 
" parameters. — 
@ Increases in leakage current occur. 


@ Changes in channel resistivity and carrier 
“mobility occur. 
Damage in MOS field-effect transistors is 
ue primarily to ionizing radiation. For this 
reason damage is reported in terms of dose (in 
rads} or exposure (in Roentgens) rather than flu- 
ence (in n’cm*). The most sensitive parameter 
to radiation in field-effect’ transistors is the 
threshold voltage. 1’,. In general. degradation in 
Vy proceeds rapidi\: in the range of 10° and 104 
rads (Si), but becomes more gradual above this 
dose. Complete failures. i.e.. complete degrada- 
tion in transconductance. have been observed at 
doses of 10° to 107 rads (Si). 
Considerable interest has been shown in 
ye use of junction field-effect transistors - 
(JFETS) in a radiation environment. since these. 
are unipolar devices and. do not depend on 
minority-carrier lifetime for operation. The 
primary effect of radiation on JFET’s is the © 
removal of carriers in the channel region, Radia- 
tion induced carrier removal is a strong function 
of resistivity. Thus, increased radiation toler- 
ances is expected from JFET’s with a high initial 
carrier concentration. The planar process has: 
allowed heavily doped JFET’s to be manufac- 
tured with the necessary control to make them 
commercially available. These heavily doped 
JFET’s have demonstrated a very significant 
improvernent in radiation hardness. e.g.. approx- 
imately 15 percent degradation in transconduc- 
tance after a neutron fluence of 7x 10!4 n/em? 
(E >J0 keV. fission). 
Thin-film field-effect transistor (TFT) 
devices have been found to be more. radiation © 
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resistant than conventional field-effect transis- Hence, increases in surface and bulk leakage cur- 
tors. Tests indicate that both the cadmium sele- rents induced by radiation may cause the device 
nide and the silicon on sapphire type TFT’s are to conduct continuously. This effect will depend 
operational at 10° rads (Si) or 101° n/em? (E > u the bias level applied in the application. 


10 keV, fission). - Since PNPN devices are used in medium- 
_ The negative resistance characteristic of to high-power applications, they cannot be com- 
the unijunction transistor depends upon the con- pared to high-frequency transistors in radiation 


ductivity modulation of a moderately high resis- _ resistance. Typical failure thresholds for PNPN 
tivity silicon bar. by means of injected minority devices range between 10! 2 and 5x 10!4 nfem? 
carriers from the rectifying emitter contact. This (E > 10 keV, fission). Some narrow base PNPN 
transistor is highly sensitive to radiation induced devices have performed well at 10" n/cem? (E> 
changes in minority-carrier lifetime. and resis- 10 keV, fission). 

tivity. Typical failure thresholds for unijunction 

transistors are of the order of 5 x 10!! to 5 x 9-46 Heating agd Thermomechanical 


10!2 n/cm? (E > 10 keV, fission). The degrada- __ Damage 
tion is manifested by an increase in valley volt- Any. electronic component parts. in | 
age. a decrease in valley current, and increases in which sufficient energy has been deposited by 
the interbase and emitter-base resistances. the electronic system environment will experi- 
The three basic types of silicon PNPN ence a transient rise in its temperature. The per- 
devices are: the silicon-controlled rectifier, SCR: formance characteristics of most component 
the silicon-controlled switch, SCS, and the parts are sensitive to temperature. Therefore, a 
Shockley diode. All of these devices may be con- temporary perturbation in the response of elec- 


sidered to consist of overlapping NPN and PNP tronic component parts can be expected. The 
transistors. the primary difference being the severity of the perturbation is a function of the _ 
external accessibility of the various layers; that  d¢posited energy and the manner in which com- 
is, the Shockley diode provides external access ponents are interconnected and mounted. Semi- 
to only the outer P- and N-layers, the SCR has conductor devices are particularly vulnerable to 
leads to all but the central N-region, and the SCS temperature transients. . 
has leads to all four regions. The “two transis- eh Fortunately, most of the common cir- 
tors” of the PNPN structure operate in a cuit design techniques for compensation for 
positive-feedback configuration, and the temperature rises are directly applicable to the 
current-transfer ratios of the two sections add circumvention of heating effects caused by the 
together for the composite device. TREE environment. Consequently, heating ef- 
As previously discussed for transistors, fects seldom are emphasized, and other effects 
radiation induced defects reduce the current predominate. However, in some particularly 
gain for both “transistors” so that the required sensitive components such as inertial-guidance 
gate current, holding current, and breakover devices, heating effects remain a_ serious 
voltage should: increase with radiation at flu- roblem. = 
ences comparable to the bulk damage fence The thermomechanical-shock effects 
levels in silicon transistors. Theoretical consider- anse Irom the deposition of short pulses of high 
ations of the mechanisms of PNPN device opera- intensity X-ray energy. The component part 
tion also indicate that excessive leakage currents response differs from the effects discussed so far 
will cause premature triggering of the devices. in that the primary manifestation is the loss of 
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mechanical integrity. The processes of spalla- 
tion. blowoff. and delamination combine to pro- 


duce mechanical damage. In most cases. this © 
mechanical damage results in a permanent, cata- . 


strophic electrical failure of the component parts. 


ap All electronic components are poten- 
aw vulnerable to thermomechanical shock. 


but semiconductor devices are among the most 
vulnerable components. The failure. modes for 


transistors exposed to X-rays depend on the ma-. 


terials and geometries employed in their con- 
struction. It is. therefore. worthwhile to con- 
sider device fabrication in some detail. 
Transistors are composed of a combina- 
10M O! materials. and the relationship of these 
materials to each is best illustrated in terms of 
the processes by which these devices are fabri- 
cated. Transistors are produced from single 
crystal semiconducting material. which is pro- 
cessed into regions of desired type and resistivity 
to form the junctions necessary for transistor 
action. The many techniques employed to 
achieve the required junction configurations in- 
clude growing the desired material from suitably 
doped melts and alloving-in the dopant impur- 
ities from appropriately metallized surfaces. 
However. the most prevalent technology being 
used to fabricate silicon transistors is the so- 
called “planar” process in which the required 
dopants are allowed to diffuse through area- 
defining masks formed on the surface of the 
silicon. These masks are made of silicon dioxide, 
which (1) is thermall¥ grown on the surface of 
the silicon. (2) may be etched to form windows 


for diffusion. and (3) is a natural barrier to the. 


diffusion of phosphorous and boron (1he most 
commonly employed dopants for producing N- 
and P-type silicon, respectively). The windows in 
the silicon dioxide are defined by photoetching 
techniques initially developed and employed in 
the fabrication of etched wiring boards and 
‘subsequently refined in resolution to permit 
application in the fabrication of semiconductor 
devices. 


| Figure 9-57 shows the steps employed in 
ne fabrication of a typical transistor by the 
planar process. The oxide layer grown in Step 


(a) is removed in a selected area by photoetching 
_(Step b), and a P-type dopant is allowed to dif- 


_ fuse into the starting N-type. silicon to form 


what will be the base region of the transistor 
(Step c). A new oxide is grown next (Step d) and 
is photoetched to form a window over a smaller 
area. and an N-type dopant in high concentration 
is allowed to diffuse into this portion of the P- 
type region (Step e), converting it to (N* )-type 
silicon. This forms the emitter region of the tran- 


_Sistor. Depth of penetration of the dopant ma- 


terials is controlled by the times and tempera- 
tures at which the diffusion operations are 
conducted. A final oxide growth (Step f) and 
selective photoetching (Step g) expose the emit- 
ter and base contact points. The whole surface is 
subsequently coated with aluminum (Step h) 
and selectively etched to produce the desired 
current paths and Jand areas (tabs) for connec- 
tion (Step i). 

Since transistors. fabricated in the above 
manner usually are prepared in multiple array’s 
containing hundreds of identical devices, the 
next steps are the separation of the individual 
devices by some technique such as diamond 
scribing and the mounting ofeach of the re- 
sultant chips (or dice) into a separate package. 
Typical chips for high-frequency transistors are 
ess than 0.025 x 0.025 in. in area. 

The transistor package can take many 
forms (see Figure 9-58a), but in each form there 
is usually a section of the package that contains 
the lead wires to which the emitter, base. and 
collector portions of the device must be con- 
nected. This section is called the header. In a. 
typical package this is the bottom section of the 
can, which contains the three lead wires. A cap 
is attached to complete the package. The header 
usually is made of Kovar, a nickel-iron alloy fre- 
quently used for glass-to-metal seals because of — 
their comparable thermal-expansion coefficients. 
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The three Kovar leadout wires are electrically 
insulated from each other by glass feed-throughs 
in the Kovar header. Since Kovar oxidizes easily, 
the external leads, and often the internal leads, 
on the transistor usually are plated to facilitate 
su uent soldering or welding. 

The processed semiconductor chip is 
attached to the header by a fabrication step 
known as die bonding. Die bonding involves 
soldering, brazing. or glazing the die to the 
header. This attachment serves as a mechanical 
contact. a thermal path, and in some cases (such 
as the example being employed) an electrical 
contact. When the connection is to a metal por- 
tion of the header, the attachment usually is 
accomplished by brazing using a eutectic brazing 


alloy (generally containing gold). This is superior ' 


to a soft solder because the high melting point 
of the eutectic permits use of higher tempera- 
tures in sealing the top on the header and in 
preconditioning the devices if desired. Some dis- 
advantages in the use of eutectic alloys for die 
bonding are that they generally are expensive, 
quite brittle. difficult to form into unusual 
shapes, and they cannot be vacuum evaporated 
to form thin films. A preform. of the eutectic 
alloy is placed on the header on a platform di- 
rectly connected to the collector leadout wire. 
In many applications, a simple gold plating on 
the Kovar platform suffices as the preform, since 
on heating in contact with the silicon. chip, a 
gold-silicon eutectic will form. Heat is applied 
(390°-400°C), the transistor chip is placed on 
the preform, and the eutectic is chilled with a jet 
of nitrogen. This operation attaches the chip to 
header and forms the chip-to-header bond (Step 
j, Figure 9-57). 

Wire bonding makes electrical contact to 
e emitter and base regions. In some devices 
this connection is made by attaching the bridg- 
ing wires directly to suitably metallized etched- 
open areas on the silicon chip. In many devices, 
however, these areas are much too small to per- 
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mit forming quality connections, so the connec- 
tions are made to extended tabs formed in the 
aluminum metallization on top of the silicon’ 
oxide. This latter configuration is illustrated in 
Figure 9-58b. The bridging wires, 1/2 to 3 mils 


_in diameter, may be composed of any of a num- 


ber of materials, but generally they are either 
gold or aluminum. Gold is used because it forms 
ohmic contact when alloyed with silicon, is 
available as extremely fine wire with reasonable 
strength, and may be bonded by thermocom- 
pression, Aluminum wire is used because it pro- 
vides a one metal system when aluminum metal- 
lization is employed for extended tabs and 
avoids the formation of gold/aluminum alloys 
(“purple” and ‘‘white” plagues), it is available as 
fine wire, and it makes ohmic contact to both 


P-type and (N* )-type silicon. 
= In thermocompression bonding, the two 


metals (such as the wire and the bonding pad) 
are made to seize without a third intermediate 
phase (such as solder) and without melting. This 
is accomplished with high pressures and tem- 
peratures. The high temperature keeps the 
metals in the annealing range as they flow into 
atomic intimacy under the bond. One form of 
thermocompression bonding (called ball bonding 
or nail-head bonding) employs gold wire that is 
heated at. the tip with a hydrogen flame to form 
a ball, which is subsequently driven against the 
heated chip under pressure to form the bond 
and concurrently is flattened into a nailhead 
configuration. In stitch bonding (another form 


of thermocompression bonding), the wire is bent 


under pressure from the tool head and. forced 
into the heated pad in such a way that the wire 
is wedged flat at the point of contact. Either 
gold or aluminum may be employed in stitch 
bonding. In ultrasonic bonding, ultrasonic 
energy is employed to force the two metals into 
atomic intimacy. The resultant bond resembles a 
stitch bond. These. three bonding techniques are 
illustrated in Figures 9-59a, b, and c, respectively. 


Several general. modes of failure 
(see Figure 9-60) are apparent: (a) the wires may 
break: (b) the emitter or base bonds may 
separate from the chip: (c) the silicon chip might 
fracture: (d) the chip may separate from the 


header: (e) the aluminum metallization. consti- - 


tuting the extended tabs might become discon- 
tinuous; and (f) the wire bonds at the post may 
eparate _ 
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Fe OTHER ELECTRONIC 


COMPONENT PARTS 
9-47 Electron Tubes 


The electron tubes discussed in this sub- 
section are divided into three categories:: vacuum 
tubes, gas-filled tubes. and phototubes. These 
categories include nearly all types of conven- 
tional electron tubes with the exception of 
microwave tubes. Because of the known dif- 
ferences in interaction of nuclear radiation with 
metals, gases. and photoemissive materials, as 
well as with glasses and other ceramics, the ob- 
servable effects of these interactions might be 


’ observed in 


expected to differ for each category of tube. 
Permanent damage to vacuum tubes is 
associated with thermal and fast neutron expo- 
sure. No significant permanent damage has been 
conventional vacuum tubes or 
Nuvistors exposed to test radiation pulses or to a 
thermal neutron fluence of less than 10!? 
n/em? (fission). Moderate damage may occur in 
standard size glass or metal tubes at thermal! nev- 
tron fluences from 10} 4 to 10!© n/em? (fis- 
sion), and severe damage may occur at fluences 
greater than 10!° n/cm? (fission). “Severe dam- 
age” is often the failure of the glass envelope. 
usually at a glass-to-metal seal: “‘moderate dam- 
age” is generally a permanent change in tube 
operating characteristics. Miniature, sub- 
miniature; and ceramic-type tubes are less sub- | 
ject to both moderate and permanent damage 
than standard size tubes, largely because of the 
smaller area-and mass of the tube parts. How- 
ever, thermal-neutron fluences sufficient to 
cause permanent damage to high-vacuum tubes 
are greater than neutron fluerices expected to be 
encountered in the nuclear weapon environ- 
ment. e.g., >10! 6 n/em? (fission). 

The principal transient effect that results 
from exposure of vacuum tubes. including the 
Nuvistor, to nuclear radiation is produced by 
Compton scattering of electrons from structural 
members by gamma rays. Leakage currents 
caused by ionization of air between external 
electrodes, or reduction of resistance of insulat- 
ing. materials such as glass, ceramics, and mica, 
caused by electron excitation within the ma- 
terial, may have minor effects on tube perfor- 
mance. There is no appreciable radiation in- 
duced liberation of gas from tube parts, or ioni- 


zation of residual gases within the tube. 
as Most of the Compton electrons pro- 


uced in the structural parts of the tube and 
ejected into the evacuated region are too ener- 
getic to be influenced significantly by the elec- 
tronic fields in the tube. However, the impact of 
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high energy electrons on the interior surfaces of 


the tube assembly produces low energy secon- ‘ 


dary electrons that can be influenced by the 
_ existing electric fields, and thus can alter the 
normal operating characteristics of the tube. The 
ejected electrons can be collected by one or 
more electrodes, depending upon the electrode 
potential and position. The magnitudes of the 


resulting transient voltages that appear at the’ 


respective electrodes depend on the magnitudes 
of the transient current and the circuit resis- 
tance. The grid circuit is particularly affected by 
this phenomenon, since it usually. suffers a net 
loss of electrons and therefore may assume a 
positive charge. The resulting increase in plate 
current is determined largely by the grid resis- 
tance and the gain of the tube. 

Gas - filled tubes (thyratrons), under 
extensive neutron bombardment, most likely 
will fail by breaking of the glass envelope or 
glass-to-metal seals. Such damage occurs at ther- 
mal neutron fluences exceeding 10!® nfcm? 
(fission). 

The principal transient effect in a thyra- 
ron subjected to a nuclear radiation pulse is 
spurious firing caused by ionization of the filling 
gas. The filling gas, in these cases xenon, be- 
comes partially ionized, primarily by gamma 
rays. Additional ions are created by ion-neutral 
molecule collisions in the electric field between 
the plate and the grid. A positive ion sheath can 
form around the negatively-biased control grid, 
which will neutralize the grid charge and will 
permit electrons to be accelerated from the 
cathode space charge toward the plate. As ion 
density increases, a sustaining discharge ensues, 
which can be shut off only by eee the 
plate voltage. 


. wnt Phototubes are designed to exhibit peak 
nsitivities to electromagnetic radiation in the 
visible and near-infrared regions by choosing 


photosensitive cathodes with low work func- 
tions. The presence of photosensitive material in 
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a vacuum tube introduces effects other than 
ane found in receiving type eens tubes in a 
igtion field. 

Severe permanent damage to baste 
tubes, as with other vacuum tubes, is attribut- 
able to large thermal neutron fluences. Severe 
damage starts in phototubes at about 10! 
n/cm? (fission). However, moderate permanent 
damage may occur. at thermal neutron fluences 
two orders of magnitude lower. Moderate dam- 
age in some cases is an increase in dark current, 
or a decrease in anode luminous sensitivity, but 
in most cases it is a darkening of the glass 
envelope. This darkening effectively reduces 
tube sensitivity. The glass discoloration has been 
observed at fast neutron fluences of 5.5 x 10!? 


(E > 10 keV, fission), and gamma ray doses of 


6.3 x.10® rads (C). Thus, both gamma and nev- 
tron components of mixed radiation. contribute 
to permanent damage. 


The principal transient. effect in photo- 
multiplier exposed to pulsed X-ray radiation 
(and presumably to gamma radiation) is an in- 
crease in anode current. The increase can be as 
much as the space-charge-limiting value for a 
given tube. Furthermore, the duration of the 
current increase is much greater than: that of the 
radiation pulse. At first, the current mcreases as 
a result of currents initiated by luminescence of 
those areas of the glass envelope that are opti- 
cally coupled to the cathode of the photomulti- 
plier, This current increase also has been demon- 
Strated. in steady-state gamma fields. Glass 
luminesces at wavelengths and intensities deter- 
mined by the glass composition and by the dose 
rate. Since a great deal of the radiant energy is in 
the visible portion of the light spectrum, where 
common photomultipliers are sensitive, it would 
be expected that the degree of photomultiplier 
response to X-ray and gamma radiation would 
depend upon both the type of glass used for the 
envelope and the spectral sensitivity of the 
cathode material. 


A second mechanism is required to ex- 
plain the relatively slow decay of anode current 
after the radiation field is removed, since glass 
luminescence decays rapidly. It is believed that 
decay of anode current may be retarded by 
electric-field changes in the tube when electrical 
insulators become charged as a result of large 
initial current pulse 


9-48 Capacitors oe . 


Nuclear radiation affects most of the 
electronic properties of capacitors to some ex- 
tent. Changes in the capacitance value, dissipa- 
tion factor, and leakage resistance have been 
observed during’ steady-state. reactor experi- 
ments. These effects generally are not consid- 
ered severe for fast neutron fluences less than 
10/* n/cm? (E > 10 keV, fission), and for most 
capacitors this limit is about 10!? n/em? (£ > 
10 keV, fission). , 

During a high-intensity pulse of nuclear 
radiation, the most pronounced effect in a 
capacitor is a transient change in the conductiv- 
ity of the dielectric material with a correspond- 


ing increase in the leakage currents through the - 


capacitor. The most recent concept of ionization 
effects in insulating and dielectric materials indi- 
cates that the ionizing particles create ionizing 
tracks in the irradiated materia]. This means 


_ that, microscopically, the material is not uni- 


formly ionized (an exception to this occurs at 
very high dose rates, =10!? rads (Si)/sec, where 
there should be sufficient overlap of the ionized 
tracks for the material to be considered uniform- 
ly ionized). 


The excess conductivity induced ina’ 
material irradiated with a short pulse of ionizing 
radiation is generally classified in two compo- 
nents: the prompt component, and the delayed 
component. The prompt component is primarily 
the result of excess carrier concentration from 
direct ionization by the radiation and the con- 
current recombination and trapping of these car- 


riers. The delayed component is that component 


of conductivity that remains after the termina- 
tion of the ionizing pulse. This does not mean to 
imply that it does not make a small contribution 
during the radiation pulse. The delayed compo- 
nent is the result of thermal generation of excess 
carriers from shallow traps, in which they are 
caught during the prompt pulse. and. their con- 
current loss to recombination and retrapping. 
The rate at which these carriers are thermally 
regenerated depends on the energy level of the 


_trap site, the concentration of filled traps, and 


the temperature: As there is usually more than 
one energy level trap in a material, more than 
one regeneration rate usually is observed in the 
delayed component. 


The excess conductivity is proportional 
to the number of carriers available to drift under 
the influence of the applied electric field. 
However, the microscopic nonuniformity of the . 
carrier concentration must be considered for the 
pulsed irradiation case. For irradiation with 
ionizing particles with a low specific ionization 
(the ratio of the number of ion pairs produced 
per unit path length to the number produced per 
unit path length by a minimum ionizing par- 
ticle), the excess prompt conductivity will be 
the same as if these carriers were generated uni- 
formly throughout the material. However, if the 
specific ionization is increased (by bombard- 
ment with more heavily ionizing particles), a 
point will be approached where the separation 
of ionization sites is less than the distance travel- 
ed by the electron before it has become thermal- 
ized and is able to drift under the influence of 
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the applied electric field. As. this point is 
approached, the probability that the electron 
will be captured in the field of a neighboring ion 


increases, and the contribution to excess con- . 


ductivity will be reduced. Thus, a plot of 
prompt conductivity as a function of specific 
ionization would show the prompt conductivity 
constant at low specific ionization and decreas- 
ing slowly after some threshold values of specific 
ionization is reached. 
The rate at which carriers are lost con- 
current with their generation by the ionizing 
radiation is proportional to the concentration of 
recombination centers and unfilled trapping 
centers. While an insignificant number of the 
total traps in the material might be filled at low 
doses, the concentration of filled traps within a 
track depends only on the specific ionization. 
Thus, the trapping rate is affected by the 
specific ionization. The result of this effect is to 
cause an increase in prompt conductivity with 
specific ionization, which would serve in part to 
compensate for the decreasing effect mentioned 
above. However, this effect on the carrier loss 
rate should be slight, since most of the carriers 
are lost to recombination rather than trapping. 
When the radiation is delivered in a time 
ort compared with the regeneration time of 
carriers from the traps, and the dose delivered in 
the pulse is large enough that significant num- 
bers of tracks near the end of the pulse overlap 
tracks generated earlier, the concentration of 
' filled traps in a track late in the pulse is different 
from that in a track created earlier in the pulse. 
When this occurs, the observed prompt conduc- 
tivity becomes a function of the total dose de- 
livered in the pulse, as well as of the specific 
ionization of the irradiating particle. It should 
be noted that a sufficient fraction of traps in a 
track must be filled to significantly affect the 
response. Hence, it is quite possible that trap 
densities in many insulators are high enough so 
that this condition is not realized in most pulse 
experiments. 
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The delayed conductivity component 
epends on the rates of carrier regeneration and 
retrapping from trap sites, which depend on the 
concentration of filled traps and the energy 
levels of the traps. The concentration of filled 
traps is a function of the specific ionization of 
the irradiating particle and, in the case of over- 
lapping tracks, the total dose. As the initial con- 
centration of filled traps within a track is usually 
a significant fraction of the total concentration 
of traps, the retrapping probability changes dur- 
ing the time the traps are emptying, thus altering 
the characteristic time for emptying the remain- 
der of the filled traps. As a result, the decay of 
the delayed component does not usually follow 
a simple law. Only in certain cases, where the 
trapping probability is negligibly perturbed by 


_ the radiation, will simple exponential decays b 


observed. ' 4 
= Neutrons produce ionization by a num- 
ber OF collision processes that give rise to ioniz- 
ing secondary particles. These processes include: 
@ Elastic scattering when the recoil atom 


receives sufficient energy to produce ioni- 
zation. 


® Inelastic scattering. producing a recoil atom 
that may or may not ionize but that emits 
a gamma photon that can produce a 
secondary ionization. 

®@ Capture, resulting in the emission of a 
photon and/or an ionizing secondary par- 
ticle (primarily thermal neutrons). 

@ Reactions resulting in an ionizing particle, 
e.g., (",pP), or (m,@) reactions (high-energy - 
neutrons). , 


There are, therefore, many possible different 
specific ionizations associated with ionized 
tracks in neutron bombarded materials. 

} In hydrogenous materials, the principal 
ionization is caused by recoil protons, which 
have a high specific ionization. For this reason, 
neutron induced conductivity in hydrogenous 
dielectrics has been found to be approximately 


one-fifth to one-half that of gamma ray induced 
conductivity tor equal ionization energy deposi- 
tion rates. For nonhydrogenous dielectrics. the 
most important contribution to neutron induced 
ionization. is by the interactions of very high- 


energy neutrons (E > 2 MeV). 


Be A “polarization effect” that is attributed 
space charge buildup within the dielectric. 


material due to nonuniform trapping has been 
observed with some capacitors, particularly with 
Mylar. mica. poly carbonate, tantalum oxide and 
Vitamin Q devices. This effect is manifested in 
several ways. One is an apparent decrease in the 
induced conductance with sequential radiation 
pulsing. Charge transfer across the dielectric dur- 
ing a radiation pulse builds up a space charge 
field opposing the applied electric field. Hf the 
applied electric field is then removed, subse- 
quent radiation pulses result in a current in the 
extemal circuit opposite in direction to that 
observed with the field applied. This is caused 
by the discharge of the space charge field. 
Similarly. if the electric field is reversed rather 
than removed after the space charge has been 


built up. the space charge field enhances the. . 


applied field. and a larger current results than 
would be observed normally. 

Saturation of the polarization effect, 
where no further decrease in the charge transfer 
is observed with subsequent radiation pulses, 
occurs after one or more pulses, depending on 
the capacitor and on the dose delivered in each 
pulse. Decreases of 50 to 70 percent for mica, 
10 to 20 percent for tantalum oxide, and 30 
percent for Mylar have been observed due to this 
ace charge buildup during radiation pulsing. 


9-49 Resistors ee 


ae Radiation effects in resistors are gen- 

erally small compared with effects in semicon- 
ductors and capacitors and are usually neglected. 
However, in circuits requiring high precision 
resistors transient effects may be significant at 
dose rates of as low as 107 rads (C)/sec and at 
neutron fluences of 10!4 n/cm? (E > 10 keV, 
fission ). 

The transient: effects are generally at- 
tributed to gamma rays that interact with ma- 
terials to produce electrons, primarily by the 
Compton process: however, energetic neutrons 
can. also produce significant ionization. 


Transient effects include (1) a change in the ef- _ 


fective resistance due to radiation induced leak- 
age in the insulating material and the surround- 
ing medium, (2) induced current that is the 
result of the difference between the emission 
and absorption of secondary electrons by. the 
resistor materials. and (3) change in the con- 
ductivity in the bulk material of the resistor. 
There is no substantial evidence, however, that 
the third effect is a first-order transient effect. 
The permanent effects are generally 
caused by the displacement of atoms by neu- 
trons, causing a. change in the resistivity of the 
material. a7 3 


9-50 Batteries and Cables 


Batteries are affected much less by radia- 
tion than other component parts. The effects of 
radiation on nickel-cadmium batteries appear to 
be insignificant at dose rates up to 10’ rads 
(air)/sec. No radiation damage was apparent in a 
number of batteries and standard cells that were 
subjected to 10!? n/em? (£ > 10 keV, fission). 


_ Transient radiation effects on an ammonia fuze 


indicated that pulsed gamma ray irradiation of 
108 rads (air)/sec had no effect on the operation 
of the battery. 
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It has been recognized for some time 
that intense pulses of radiation produce signifi- 
cant perturbations in electrical cables and wir- 
ing, including coaxial and triaxial signal cables. 
Even with no voltage applied to a cable, a signal 
is seen when the cable is exposed to a pulsed 
radiation environment. The current associated 
with this signal is defined as a replacement cur- 
rent, since it is most hkely a current in an 
external circuit that is necessary to replace elec- 
trons or other charged particles that are- knocked 
out of their usual positions by the radiation. The 
replacement current definition also applies to 
the effect of charged carriers associated with the 
incident radiation embedded in a test sample. 

The magnitude of the radiation induced 
signal varies with the voltage applied to the 
cable. This voltage-dependent portion of the 
signal, i.e., the total signal exclusive of the re- 
placement current, is called conduction current, 
thus it is ascribed to the conductivity induced in 
the insulating dielectric by the radiation. How- 
ever, it may also contain major contributions 
from polarization or depolarization processes in 
the dielectric. These can usuaily be identified by 
their gradual disappearance (saturation) after 
repetitive exposures and by their reappearance 
in additional “‘shots” in which the applied volt- 


lonizing radiation of any type produces 


ity of the material. Hence, insulators are ex- 
pected to have a transient enhanced conductiv- 
ity in an ionizing radiation environment. 


Conduction in the insulator is frequently charac- 


terized by two components: for very short radia- 

tion pulses, a prompt component whose magni- 

tude is a function of only the instantaneous 

exposure rate, and frequently at the end of the 

short radiation exposure, a delayed component 
_ having approximately exponential decay. 


Although the replacement, conduction, 
and polarization currents are fairly well under- 


9-15 


age is changed greatly, e.g., removed or reversed. 


Tee electrons that contribute to the conductiv- 


stood in terms of the interactions between the 
ionizing radiation and the metal-dielectric target 
system, it is not yet possible to predict quantita- 
tively the. response for a given cable in a speci- 
fied environment. In a mixed neutron-gamma 
environment, the induced replacement. current 
usually contains positive and negative compo- 
nents, and may therefore assume eithey polarity. 
The conduction current sometimes exhibits a 
rather complicated time dependence consisting 
of prompt- and delayed-conductivity contribu- 
tions. The polarization current appears to be 
greatly affected by the properties of the metal- 


dielectric interface. 
Permanent damage effects in cables and 


wiring are manifested as changes in the physical 
and electrical properties of the insulating ma- 
terials. When such damage becomes appreciable, 
e.g., when the insulation resistance is reduced 
severely, electrical characteristics may be af- 
fected. The extent of the damage to insulating 
materials is an increasing function of neutron 
fluence, exposure or dose, humidity, and irradia- 
tion temperature. Certain types of wire insula- 
tion are quite susceptible. to permanent damage. 
For. example, silicon rubber becomes severely 
cracked and powdered after approximately 2 x 
10'5 n/cm? (E > 10 keV, fission). The approxi- 
mate damage thresholds for three common types 
of cable insulation are: polyethylene, 1 x 107 
rads (C); Teflon TFE, 1 x 10% rads (C); and 
Teflon FEP, 2 x 10° rads (C). On the other 
hand, some irradiated polyolefins are capable of 
withstanding up to 5 x 10? rads (C). A consider- 
able degree of annealing has been observed with 
respect to insulation resistance, which implies 
the possibility of adequate electrical service- 
ability after moderate physical damage. 

It is not expected that radiation effects 
on wiring with thin insulation will exhibit the 
strange behavior observed.in coaxial cables. In- 
particular, the very limited measurements that 
have been performed indicate that the replace- 


ment current is primarily a function of the 
gamma environment. To a good approximation, 
it can be assumed that the replacement current 
for a wire ur most other objects placed in the 
radiation environment will amount to the emis- 
sion of a number of electrons between 1 and 5 x 
10°° times the number of gamma photons 
traversing the object. 

| The conduction current is a very sensi- 
tive function of the amount of insulation around 
the wire and its immediate environment. For a 
bare wire in air with a ground plane nearby, con- 
duction is due predominantly to the ionization 
produced in the air. Placing insulation around 
the wire reduces this conduction. but at the 
price of increasing the area of the wire and 
hence the effective replacement current. | 


9-51 Quartz Crystals oe 


The radiation response of a quartz 
crvstal oscillator is primarily-a function of radia- 
tion dose. The type of material from which the 
oscillator is fabricated, e.g., natural quartz, 
Z-growth synthetic. Z-growth swept synthetic, 
etc., and to a lesser extent the type of cut, ¢.g.. 
AT. BT. etc.. and the frequency and mode of 
operation determine the sensitivity of the oscil- 
lator to the radiation. The primary effect of the 
radiation is a shift in the frequency of the oscil- 
lator. Both transient and steady-state shifts have 
been observed. , 

The steady-state frequency offsets are a 
re of changes in the elastic stiffness constants 
of the crystal. For example. perturbations in the 
crystal bonds due to charge trapping at defects 
or to formation of new defect complexes will 
result in steady-state frequency offsets. Of ma- 
terials tested to date, Z-growth swept-synthetic 
quartz has been the most radiation tolerant to 
steady-state frequency offsets. Swept natural 
quartz is slightly more sensitive and unswept 
natural quartz and unswept synthetic quartz, 
respectively, are even more sensitive. The 


response of a particular crystal also varies with 
the manner in which the crystal is mounted, the 
material used for the electrodes, and to. a lesser 
extent differs for each quartz bar grown, even 


‘among bars grown under similar conditions. 


Both swept synthetic and natural quartz crystals 
can recover 80 percent to 90 percent of their 
original frequency change. after annealing at 
500°C for times on the order of 100 to 160 
hours. 

bs Transient shifts in the frequency of an 
oscillator and reduction or cessation of the out- 
put result from energy deposition and. any subse- 
quent temperature rise in the crystal. Tempera- 
ture gradients in the crystal due to faster 
removal of heat near the support wires and due 
to the electrodes absorbing more energy than 
the crystal, can give rise to frequency shifts 
induced by the resulting strain. 

Irradiation of general purpose crystal 
units has shown that they do not suffer signifi- 
cant permanent effects, within the limits of their 
stability, at a-neutron fluence of 10!3 n/om? (E 
> 10 keV, fission) and a gamma dose of 4.4 x 
10? rads (air). Transient phase and amplitude 
changes resulting from this environment are not 
of sufficient magnitude to cause concern about 
their operation under such conditions. 

Moderate precision crystal units display 
negligible frequency and amplitude changes 


when subjected to a reactor pulse. Weapon tests ~ 
as well as steady-state gamma source tests indi- 


cate that a gamma dose greater than 104 rads 
(air) is required to induce significant permanent 
frequency changes in these devices. Frequency 
changes up ‘to almost one part in 10° were ob- 
served after exposure to 7.5 x 104 rads (air) and 
10!2 n/em? at a weapon test. The possibility 
that causes other than radiation contributed to 
the changes observed at this test’ cannot be 
excluded. 
High precision natural quartz-crystal 
units either stop oscillating or exhibit appreci- 
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able decreases in the amplitude of the output 
signal during nuclear pulses of approximately 
10!? n/em? (E > 3 MeV, fission), and 3 x 103 
rads (HQ). If oscillation stops, its cessation per- 
sists for minutes. The cessation of oscillation is 
apparently independent of the voltage, current, 
or power at which the units are being driven. 
Resumption of oscillation occurs at a reduced 
drive current and lower frequency. The drive 
current is tens of microamperes below the speci- 
fied rated drive when oscillation resumes. Fre- 
_ quency changes as high as | part in 107 have 
been observed when this type of crystal was 
exposed to 7.9 x 10!! n/cm? (E > 3 MeV, fis- 
sion), and 2.9 x 10? rads (HO). 


9-52 Solder Joints 


9-53 Infrared Detectors es 


The infrared detectors that exhibit the 
greatest sensitivity to infrared radiation are also 
the most sensitive to nuclear radiation. For 
* radiation pulses that are short compared to the 
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relaxation time, the response of a photoconduc- 
tive type of infrared detector cell is an excess 
conductance that is proportional to the radia- 


_ tion exposure. For long pulses the excess con- 


ductance is proportional to the. radiation 
intensity and the carrier recombination. time. 
Neutron bombardment causes permanent 
degradation of output-signal level and signal-to- 
noise ratio. 

Irradiation of a lead sulfide device to 1.3 
x 10°? n/cm? (E > 0.48 eV, fission) at 134°F 
revealed a 67 percent reduction of output signal 
level and greater than 40 percent reduction in 
the signal-to-noise ratio. Damage was essentially 
catastrophic after 6.9 x 10!° n/cm? (E > 0.48 


what more tolerant to neutron irradiation than 
Jead sulfide detectors. After a neutron fluence of 
1.2 x 10!4 n/cm? (E > 0.48 eV, fission) at 
135°F, the output signal level of a lead selenide 
cell was reduced by 36 percent, and the signal- 
to-noise ratio was down more than 46 percent. 
The output level was down by 96 percent after 
1.8 x 10!§ n/cm? (E > 0.48 eV, fission). Lead 
selenide cells that are designed to operate at low 
temperatures are. more sensitive to radiation 
than those that are not designed for low 
temperatures. ; 

Indium antimode photovoltaic cells that 
operate at liquid nitrogen temperature showed 
significant voltage signals at doses less than 0.88 
rads (air). The cells exhibited radiation induced 
voltages roughly proportional to the logarithm 
of the dose when the radiation was delivered in 
short pulses. Recovery to within 2 percent of 
the maximum voltage occurred within 175 psec 
after the highest intensity radiation pulses. A 
complete loss of output from these cells has 
been observed after a neutron fluence of 2 x 
1016 n/cm? (E> 0.48 eV, fission). _ 

Thermistor-bolometer infrared detectors 
are the most neutron tolerant of the devices that 


have been tested. After 9.6 x 10!3 n/om? (E> 


e ssion). 
Lead selenide devices appear to be some- . 


0.48 eV. fission). the output signal level was 


down by 28 percent. and the signal-to-noise © 


ratio was down about 11 percent. After 1.6 x 
10!© njcm* (E > 0.48 eV, fission), the output 
level was down 57 percent. and the signal-to- 
noise ratio was down 66 percent. The detector 


may be usable for some applications under these 


nditions. 

Thermomechanical shock effects in 
infrared detectors will be similar to those dis- 
cussed for semiconductor devices. and wil] occur 
at about the same levels. 
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9-54 Radiation Response of Discrete- 
Component-Part Circuits ea 


GD dcterinination of the response of a cir- 


cuit is complex because of the large variations in 
circuit configuration and in the component parts 
that can be used within a circuit configuration. 
- Therefore. determining circuit response becomes 
a problem of detailed circuit analysis and/or 
testing. The radiation effects material necessary 
for this kind of analysis and testing are beyond 
the scope of this manual. Guidance may be ob- 
tained from the TREE Handbvuok and the TREE 
Preferred Procedures (see bibliography ). General 
circuit effects and typical analysis techniques are 
discussed in this manual]. 

The transient effects that can cause a 
system to malfunction can result in circuit 
responses that. like component-part responses, 
can be both dose and dose rate dependent. If the 
radiation pulses are short with respect to 
component-part recovery times the circuit time 
constants, the circuits will integrate the effects 
and will be sensitive to the total dose rather than 
the dose rate. However, when the pulse widths 
are wide, the circuits are dose rate sensitive. This 
can be illustrated by reference to Figure 9-61. 
Assume that for the circuit response plotted, the 
malfunction threshold is 1.5 volts. Therefore, a 


0.15 ws pulse at 1 x 10° rad (Si)/sec will cause a 
malfunction: It is obvious from the curve that 
the malfunction dose rate is much lower for 


-wider pulses. 


Discrete digital circuits and circuits that 
contain silicon controlled rectifiers (SCR’s) have 
displayed failure at doses as low as 0.1 to 1.0 rad 


' (Si) for short pulse widths. It is difficult to 


design a circuit specially that will not malfunc- 
tion above a prompt dose of 100 rads (Si). 
Another effect that can result from the 
ionizing radiation is the initiation of a cata- 
strophic action or catastrophic failure. An 
example would be the firing of a pyrotechnic 
device or the premature initiation of a firing 
signal. A second type of catastrophic action 
occurs if a circuit destrovs itself as a result. of: 
the effects caused by ionization (burnout). Fig- 
ure 9-62 shows an output stage of a power 
supply inverter. In normal operation. Q1 and Q2 
are turned on alternately. The output at the 
transformer secondary is a square wave. Ionizing 
radiation may cause Q] and Q2 to tum on 
simultaneously. After the radiation pulse, the 


‘transistors will recover to normal operation, and 


one of the transistors will attempt to turn off. 
At this time a large voltage will be induced 
across the transistor that is attempting to turn 
off. If this voltage exceeds the breakdown volt- 
age, the transistor may be damaged. 
General statements applicable to both 
permanent and transient effects include: 
® Circuits that use low frequency, thick base 
semiconductors usually are more suscepti- 
ble to radiation effects than those circuits 
that contain high frequency, thin base 
devices. . 
@ Germanium devices generally will show 
larger photocurrents and leakage currents 
than comparable silicon devices 
@ High impedance circuitry generally will be 
more susceptible to radiation effects than 
low impedance circuitry 
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Figure 9-61. 


Circuit Response asa Function of 


Radiation-Pulse Width 


@ Magnetic memory devices will not be sensi- 
tive to the neutron and gamma environments 
typically specified as the environment in 
which an electronic system must survive. 


The primary permanent effects on cir- 
cuits will be the degradation of semiconductor 
devices. The solid-state power supplies and 
regulators with their low frequency transistors 
will fail to perform their required function when 
exposed to fluence between 10!! and 10!3 
n/om? (E > 10 keV, fission) depending on the 
circuit configuration, component parts, and 


design margins in the system. Circuits that use — 


MOSFET can fail at gamma doses between 107 
and 10° rads (Si). 
The failure threshold for thermo- 
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mechanical shock will be established by the 
t old of the component parts. 

It is desirable to have a method of anal- 
ysis that can be used to predict the response of 
components, circuits, and systems. The analysis 
methods that are used consist of established 
techniques to calculate circuit and system 


responses by replacing radiation effects by their - 


corresponding electrical effects. Thus, the prob- 
lem that is solved eventually is wholly electrical 
in character. A major advantage of analysis as a 
simulation tool is that it is not necessary for a 
circuit or system to exist in a physical state 
before it can be analyzed. In addition, the 
analyst has control over the “environment,” and 
it is theoretically possible to simulate the total 


Figure 9-62, G Circuit in Which Burnout Could Occur SB 


environment during a single (but complicated) 
analysis. Perhaps the major disadvantage of anal- 
ysis (for TREE) is the relatively low confidence 
in the results. This low confidence level usually 
arises from the simplifying assumptions that are 
often made to expedite the analysis and which, 
themselves, are subject to verification. typically 
by environmental testing. , 
Certain requirements or inputs are 
needed for any circuit or system analysis. First, 
an accurate mathematical description of the 
electrical characteristics must -be obtained. Such 
a description usually is checked by comparing 
the computed electrical response to the mea- 
sured electrical response of a circuit or system. 
Second, the radiation effects on electronic ma- 
terials and devices must be represented, or 
modeled, by electrical effects. For example, dis- 


placement effects may be modeled by. making 


transistor current gain a function of time. This 
step generally requires environmental. testing 
and/or sound analytical procedures to obtain the 
required radiation effects data. Third, an anal- 
ysis method must be employed to make -an 
accurate calculation of the steady-state and 
transient responses of the electronic network of 


interest. This may be done by hand or with the 
aid of a computer.* 

Hand analysis techniques are useful for a 
quick qualitative and, to a limited degree, 
quantitative appraisal of the sensitivity of linear 
circuits and logic circuits to radiation environ- 
ments. Manual techniques are valuable in the 


prediction of permanent effects of transient 


radiation, particularly when the’ relevant elec- 
trical parameters assume constant, degraded 
values after irradiation. The hand analysis tech- 
niques are often quite adequate to establish the 
initial estimate of radiation induced voltage and 
current transients and of the steady-state perfor- 
mance degradation. This type of analysis is suit- 
able for a rough estimate of the peak-amplitude 
radiation response. In practice, only small, 


' simple circuits can be handled. 
<s In doing analysis, there is an inevitable 


choice between time (man-hours) and accuracy. . 
If great confidence in the results is not required 
(e.g., when the analysis results are to be used to 


§: more detailed discussion of hand and computer analysis 
tecfiniques is contained in the “TREE Handbook” (see bibliog- 


raphy), 
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plan further environmental tests), a large num- 
ber of simplifving assumptions can be made, and 
the_analysis can be carried out quickly by hand. 
If the results of an analysis are to be 
used directly in a survivability assessment, high 
accuracy is desirable. This implies that few 
simplifying assumptions may be made. If a cir- 
cuit contains more than one or two active de- 
vices (€.g., transistors), the circuit model will be 
complex, and the speed and accuracy of a digital 
computer should be used. Both hand analysis 
and computer aided analysis require an 
equivalent-circuit model that represents elec- 
trical and radiation induced phenomena. A 
notable difference is that the computer can 
‘handle an equivalent circuit model in its entirety 
and can generate the desired response function 
without making engineering approximations. 
as Both analog and digital computers have 
een used for response predictions, and each has 
advantages and disadvantages; however, recently 
developed mathematical techniques and pro- 
gramming capabilities make the digital computer 
‘preferable for. most problems. Several digital 
programs are available for circuit analysis. These 
are described in some detail in the “TREE Hand- 
book” (see bibliography). These codes do not do 
the analysis. They do perform the tedious, error- 
prone calculations. The individual who uses 
them must provide the accurate description of 
the equivalent circuit model and must interpret 
the results. 


Two factors contro] the accuracy of the’ 


s of any TREE analysis method. The first 
factor is the accuracy and completeness of the 
description of the response of individual compo- 
nents to a particular radiation environment. This 
information is basic to the success of any anal- 
ysis technique and frequently has been the 
major stumbling block for analysis attempts. 

The second factor that affects the accu- 
racy of the analysis results is the assumptions 
that are made to simplify the analysis problem. 
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Certain aspects of the radiation response of 
individual! components generally are neglected 
by assuming that they will contribute only a 
small or negligible portion to the radiation 
response of the circuit. These assumptions are 
made on the basis of general experience with the 
radiation response of circuits and with knowl- 
edge of the limits over which they might be 
valid. Although these assumptions are generally 
correct,. there are specific instances and specific 
circuit configurations for which they may not 
hold, and care must be taken in making such 


assumptions. 
Tm practice most of the analysis ap- 
proaches result in a fairly reasonable correlation 


with experimental results. However, it is com- 


mon knowledge that any experimental result can. 


be explained by a theory, but the theory will 
not always predict the correct result in a new 
situation. Therefore, caution should be exercised 
when accepting a component representation or 
an analysis technique that predicts the results 
for a pulsed reactor environment reasonably 
accurately, if these results are to be applied to a 
nuclear weapon environment. It is also possible 
that techniques applicable to switching circuits 
or nonlinear circuits will not apply to all linear 
circuit analysis. 


9-55 Radiatio sponse of Integrated 
Circuits : 


Integrated. circuits include many circuit - 


types differing in construction materials and 
methods. The four construction types are the 
monolithic semiconductor, thin film, multichip, 
and hybrid integrated circuits. The scope of this 
section will be limited to monolithic and thin 
film circuits since the radiation response of 
both multichip and hybrid circuits can be in- 
ferred from the discussion of monolithic and 
film circuits or discrete devices. The discussion 
includes junction isolated, dielectrically isolated, 
and air isolated integrated circuits. 


“ 


In junction isolated circuits, the com- 
ponents are defined within a single crystalline 
substrate by regions of alternate doping that are 
electrically isolated by reverse biased PN- 
junction boundaries. The doped regions are 


’ formed by the geometrically controlled diffu- 


sion of appropriate impurities into the substrate. 
One or more uniformly doped, epitaxial layers 
may be grown upon the substrate prior to diffu- 
sion (planar epitaxial). The dielectric-isolated 
circuit is distinguished by the use of a dielectric 
(silicon dioxide or ceramic) instead of a PN- 
junction isolation between. critical components. 
A single component or a small number of com- 
ponents are formed within individual single 
crystalline islands (called tubs) that are im- 
bedded in a polycrystalline substrate. The active 
elements in both the junction and dielectric iso- 
lated circuits can be bipolar transistors, junction 
FET's. or insulated gate (MOS) FET’s. The air 
isolated circuit usually employs aggregates of 
-unipolar (i.e.. field effect) transistors of metal- 
oxide-silicon construction.* Since this type of 
transistor may be used as a bias-dependent resis- 
tor. complete designs usually are constructed 
without the. use of other circuit elements. Air 
isolated MOS integrated circuits are. fabricated 
by growing silicon on sapphire (SOS). Portions 
of the silicon are etched away. leaving isolated 
islands of silicon upon which transistors are 
fabricated. 

Thin-film integrated circuits employ 
‘geometrically controlled surface films of con- 
ductive and dielectric materials upon a glass or 


_ ceramic substrate to define passive circuit ele- 


ments.and interconnections. The active-elements 
may be formed as an integral part of the process 
(thin film, insulated gate, field effect transistors) 
or welded to the circuit (conventional, discrete 
transistors). Circuits of the latter type are re- 
ferred to as hybrid thin film circuits. 

The categorization of circuit types given 
above is somewhat arbitrary. It is based on the 


present. developments in: the integrated circuit 
industry, rather than on strict lines of variance 
between the types. The categorization is used 
for convenience of discussion with respect to the 


‘effects of transient radiation, where distinctly 


different effects may occur. For example, the - 
PN junction used for isolation is a source of 
photocurrents during an ionizing radiation pulse. 


and may result in large substrate (and hence. - 


power-supply) currents in junction isolated... 
monolithic circuits. This effect. is absent in 
either the air or dielectric-isolated monolithic 
circuit or the thin film, hybrid circuit. 

Transient radiation may cause both 
transient, permanent and thermomechanical 
shock effects. On a qualitative basis. the primary: 
electrical. effects introduced in integrated cir- 
cuits by transient radiation are similar to the 
effects described for conventional solid-state cir- 
cuitry. The magnitude, duration, and ‘electrical 
consequences of these effects, however, do not - 
follow directly from conventional circuit experi- 
enge. | | 

The effects of radiation on an integrated 
circuit are more closely related to its geometrical 
and physical characteristics than to its electrical 
function or circuit configuration, The proximity 
of circuit elements within the device, and in 
some cases its integral structure, make several 
modes of secondary interaction possible. This is 
especially important in the case of junction iso- 
lated integrated circuits. 

The transient effects observed in inte- 


grated circuits result from the generation of ex- 


cess charge carriers that cause photocurrents and 
voltage. changes. As previously described for 
transistors, the motion of excess carriers is 
governed by the response of carriers to electric 
fields and concentration gradients. The charge 
carriers will cause currents to occur until they 


» Wa bipolar integrated circuits have also been con- 
ted. 
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are swept out by external fields and electron- 
hole recombination. The peak photocurrents can 
be a function of either dose rate or dose depend- 


ing on the duration of the radiation pulse. As a — 


result, if the radiation pulse is long compared to 
the circuit radiation response time, the microcir- 
cuit response is dose rate dependent. However, 
for radiation pulses that are short compared to 
the circuit radiation response time, microcircuit 
response will be dose dependent. Thus, the 
width of the radiation pulse can be of consider- 
able significance, since the peak photocurrent 
generated can be a function of the duration of 


the pulse as well as its amplitude. 
me In a transient-radiation environment, the 


semiconductor. integrated circuit reacts to sev- 
eral mechanisms that have been discussed in con- 
nection with conventional circuitry. One point 


of departure that is made necessary by the | 


monolithic nature of the circuit, is the signifi- 


~ cant and often predominant interelement effects | 


that occur. in addition to the intraelement ef- 


fects. It is important to consider current paths: 


between as well as within component parts of 
the microcircuit. 

Quite generally, the transient effects in 
any integrated electronic: device are a conse- 
quence of a sequence of events that may be 
described as follows: 

® The radiation interacts with the circuit 
material and surrounding encapsulant to 
introduce charge carriers and to establish a 
nonequilibrium charge distribution, 

@ Acting under nonequilibrium electric fields 
and. concentration gradients, mobile. car- 
riers flow in the direction that restores 
equilibrium and thereby produce primary 
electrical currents. These electrical currents 
may be semiconductor-junction photocur- 
rents, replacement currents, dielectric- 
leakage currents, gas-ionization currents, 
etc. 

®@ The nonequilibrium charge distribution and 
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the primary electrical currents may interact 
with the electrical circuit to produce 
secondary effects, e.g., secondary photo- 
currents. Under certain circumstances, the 
secondary effects may be sufficiently re- 
generative to be self-sustaining, and a new 
stable circuit state will result. In addition, 
localized electrical stresses may introduce 
permanent damage. 
@ The cumulative effects of the radiation 
induced currents and circuit action are volt- 
age, current, and impedance changes of 


variable duration at the terminals of the 


integrated circuit. 
In junction isolated circuits, the -pre- 
dominant effect is PN-junction photocurrents 


resulting from ionization in the semiconductor . 


material. Important secondary effects include 
secondary photocurrents produced by transistor 


action in any three adjacent doped regions, po- 


tentially large substrate currents, and “‘latchup.” 

The predominant effect in dielectric iso- 
ated circuits is also PN-junction photocurrents. 
The major difference from other monolithic 
structures is the absence of. the extra PN- 
junction between. the components and the sub- 
strate and its associated photocurrent. Also, 
c t paths are more restrictive. 

In air isolated integrated circuits, the 
important primary transient effects are PN- 
junction photocurrents, replacement currents 
resulting from charge scattered from device lead 
wires and the case, and ionization currents 
through the surrounding encapsulant. The pre- 
dominant secondary effect is a secondary photo- 
current resulting from the radiation induced gate 
current and transistor action. In addition, photo- 
current is generated.in the Zener diode employ- 
ed for protection in the gate lead of MOS 
circuits. 


In MOS integrated circuits, the impor- 
tant primary transient effects are drain-substrate 
and source-substrate PN-junction photocurrents, 


teplacement currents resulting from charge scat- 
tered from device lead wires and the case, and 
jonization currents through the surrounding 
encapsulant. The predominant secondary effect 
is a secondary drain current resulting from the 
radiation induced gate current. In addition, 
photocurrent is generated in the Zener diode 


employed for protection in the gate load. : 
= For the most part, thin film circuits may 


e treated as conventional circuits with extreme- 
ly small geometries. Ionization. currents within 
and between elements have specific importance, 
especially in the high-impedance circuitry asso- 
ciated with thin film circuits that employ field 
effect transistors. Nevertheless. in most cases of 
‘interest, the transistor is the predominant ele- 
ment that determines the transient radiation 
response of the circuit. ; 

Semiconductor integrated circuits of the 
planar-diffused (or planar-epitaxial) type ex- 
‘perience transient effects that may be attributed 
to the interaction of the circuit elements 
through the active substrate. Two predominant 
interelement effects that must be considered are 
the presence of large substrate currents and the 
occurrence of latchup. 
In practice. the high packing density of 
elements on a substrate chip results in the 
presence of isolation diodes over most of the 
area of the chip. Thus. a chip 40 x 40-mils may 
have 1,500 mils” of effective isolation-diode 
area. The total substrate photocurrent may be 
- 100 times that of a typical diode in the circuit. 
Since the substrate is connected to the power 
supply system, the substrate currents will be re- 
flected in large currents appearing in the power 
supply leads. Radiation induced power supply 
currents of the order of 1 ampere, with dura- 
tions of a few microseconds have been observed 
at prompt doses of approximately 10 rads (Si). 


The potential hazards to the power supply sys- 


tem, which must supply many such circuits, are 
evident. 


wa Present evidence indicates that the very 
large power supply currents occur primarily in 
those circuits where transistor action through 
_the substrate is possible. In circuits of this type. 
radiation thresholds above which the current 
increases suddenly have been observed. The 
thresholds have been attributed to the turning 
on of an equivalent four layer device. Other pos- 
sibilities include second breakdown and sustain- 
ing voltage breakdown. 
Transients induced in integrated circuits 
by pulsed ionizing radiation last from less than a 
microsecond in high-speed digital circuits to 
several tens of microseconds in slower circuits. 
Occasionally, pulsed-radiation effects with. con- 
siderably longer recovery. times can be explained 
by circuit time constants. In the extreme case, 
the abnormal state persists until the dc power is 
interrupted. When this occurs, normal circuit 
operation is inhibited and latchup has occurred. 
In some cases, which are referred to as incipient 
latchup, the condition lasts only for periods that 
are long with respect to normal recovery times 
of the circuits. Latchup can be induced in three 
ways: by exposure to ionizing radiation; by par- 
ticular sequences of applying voltage to circuits 
employing more than one power supply; and 
other electrical stimulations such as high voltage 
pulses. Only radiation induced latchup is consid- 
ered here. Radiation induced latchup has been 
observed in only a small percentage of the device 
types that have been irradiated with pulsed 
ionizing radiation, and of these device types 
usually only a small percentage of the samples 
exhibit latchup. In a few cases. a majority of the 
samples of a certain part from a manufacturer of 
integrated circuit part types have. exhibited 


latchup. 

= Integrated circuit latchup is always 
caused by one or more normally reverse biased 
PN junctions becoming conductive, either by the 


initiation of a breakdown mechanism or by be- 
coming forward biased. In either case, a sustain- 
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ing mechanism must act to maintain the break- 
down or forward bias condition, once it has 
been initiated. Three latchup mechanisms have 
been postulated and observed in junction iso- 
lated integrated circuits. These are: PNPN 
action, second breakdown, and transistor sus- 
taining voltage breakdown. No other mechan- 
isms for sustaining latchup are known, although 
severgl others have been postulated. 

ee Latchup is normally associated with 
junction isolation since it usually involves some 
type of interaction with the silicon substrate 
through an isolating PN junction. Dielectric iso- 
lation is effective in isolating elements from one 
another and from the silicon substrate and, thus, 
is an important step in reducing the latchup vul- 
nerability of integrated circuits. However, anal- 
yses of the structural characteristics of certain 
dielectric-isolated circuit types have indicated 
that the possibility of latchup cannot be ruled 
out in dielectric-isolated circuits. It is possible 
that some of the same mechanisms that are re- 
sponsible for latchup in junction-isolated circuits 
also can exist in a dielectric-isolated circuit. 
These mechanisms include second breakdown, 
sustaining voltage breakdown of transistors and 
PNPN action if 4 layer structures are included 
within a dielectrically-isolated region. 

No latchup mechanisms have been found 
that are peculiar to dielectric isolation. While 
photocurrents can be generated in dielectric iso- 
lation, there are no sustaining mechanisms for 
these currents unless the isolation is defective or 
is subjected to destructively high voltages. ‘ 
- Hybrid thin film circuits may be ex- 

ed to be as tolerant of radiation as their 
conventional circuit counterparts. The radiation 
response is determined primarily by the active 
elements in the circuits. 

Undoubtedly, dielectric isolated circuits 
are much. less vulnerable to latchup than are 
their junction isolated counterparts. However, 


dielectric isolated circuits probably are more 
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latchup-prone than discrete component circuits 
because: dielectric isolated process limitations 
occasionally permit four layer structures; dif- 
fused resistors are present in some dielectric iso- 
lated circuits; and protective surface coatings 
occasionally are used in special purpose potting 
compounds or encapsulants, which might com- 


promise the isolation. Component isolation in a - 


dielectric isolated circuit, while much superior 
to that in a junction-isolated circuit, is still 
somewhat less complete than that in a discrete 
component circuit because of photocurrents 
thr the dielectric. 
A large number of integrated circuits 
have been irradiated, but the testing has been 
concentrated on specific microcircuit types, and 
a broad base of experimental data on the re- 
sponse of microcircuits to radiation is not avail- 
able. This lack of data. is especially true for 
linear circuits. Representative radiation failure 
levels for some common digital junction-isolated 
types are shown in Table 9-23. The levels are 
listed either in terms of gamma dose or of dose 


rate, depending upon whether the circuit is nor-. - 


mally dose or dose rate dependent. Design and 
production changes in integrated circuits are 
common in the industry. The broad ranges given 
reflect highly variable experimental results and 
indicate the necessity of considering each circuit 
as a separate problem. 

The proximity and intercoupling of ele- 
ments do not assume importance in the produc- 
tion of permanent effects by nuclear radiation. 


Integrated circuits may be treated as conven- 


tional circuits of small dimensions. The primary 
factor that determines the tolerance of the cir- 
cuits to radiation induced permanent effects is 
the degradation of the active elements with 


accumulated radiation exposure. 
= Changes in the electrical parameters of 
1odes and transistors that result from radiation 


have been discussed. Experiments have shown 
that the circuits will experience failure when the 


3 


Table 9-23. Representative Radiation-Failure 
Levels” for Digital Junction-lsolated : 
Semiconductor Integrated Circuits 


Radiation- 
Typet Failure Level 
Large area. slow § -20 rads (Si) 
speed 
Large area, moderate 10°-10" rads (Si)/sec 
speed 20 -60 rads (Si) 


Small area, moderate 


107-10? rads (Si)/sec 
or high speed , 


Here, the ‘‘failure”’ level corresponds to exceeding 


the circuit’s noise margin. ‘ 
‘* Type designations defined as follows: Large area ~ 
cHip area > 5,000 mil’, Small area - chip area < 5,000 
mil?, Slow speed - propagation delay > 100 nsec, Moder- 
ate speed - propagation delay between 25 and 100 nsec, 
High speed ~ propagation delay < 25 nsec. 


gain of the transistors has dropped to. the point 
that they will no longer support proper circuit 
action. The radiation resistance of the circuits is 
determined by the stability of the gain of the 
transistor elements with respect to radiation 
exposure and the tolerance of the circuit design 
with respect to gain degradation. Although no 
class of integrated circuits has been shown to be 
inherently superior to another, those ..circuits 
employing faster transistors usually can with- 
stand a greater neutron fluence. Epitaxial tran- 
sistors usually, but not exclusively, represent the 
faster transistor types. Radiation failure levels 
have been shown to vary from 10!? to 10!* 
n/cm? (E > 10 keV, fission), with the faster 
circuit types at the high end. The normally con- 
servative design and digital function of most 
integrated circuits accounts for the circuit 


longevity beyond what would be considered the 
inimum useful transistor gain point. 

The radiation induced circuit response 
of microcircuits is manifested by changes in 
both the dc and switching characteristics. The 
effect on the integrated circuit parameters of 
changes in components after irradiation will; of 
course, depend on the specific circuit configura- 
tion involved. Frequently, the radiation sensitiv- 
ity of the circuit is determined by the tolerance 
of the circuit design with respect to gain 
degradation. 

~The most radiation sensitive circuit 
parameter of digital gates and flip-flop circuits is 
the output low voltage. Circuit failures result 
when normally ON transistors leave saturation. 

_ The amount of current that an output transistor 
can sink is directly proportional to the current 
gain.. The changes in the output transistor cur- 
rent gain are reflected directly in the current- 
drive capability (fan-out) of both digital gates 
and flip-flops. a 

Changes resulting from radiation are 
observed for other digital-circuit. parameters. 
The saturation voltage of transistors, Veat: 
increases with neutron fluence, even though suf- 
ficient base drive is supplied to maintain the 
transistor in saturation, as a result of an increase 
in the saturation resistance. These changes in 
saturation resistance usually are negligible at 
threshold fluences applicable for maximum fan- 
out. The input threshold voltage of gate circuits 
normally will increase with radiation exposure as 
a result of changes in the base-emitter voltage of 
the output transistor and increased diode for- 
ward voltage. Changes in this parameter, how- 
ever, are not considered significant. Increases in 
leakage currents also have been observed with 
radiation exposure, but the changes in this 
parameter usually will not affect circuit 
‘performance, 


The switching characteristics of typical 
Igital circuits also are affected by radiation. 
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The rise and fall times of the transistor elements 
of the microcircuit are increased and its storage 
time is decreased after neutron exposure. These 
effects combine. and usually a smal] net increase 
of switching time is observed. 


card A good estimate of the radiation toler- 
ance (at different fan-out conditions) of digital 


circuits can be made by measuring the output 
current-voltage characteristics. The gain degrada- 
tion can be calculated, and the degraded charac- 
teristics can be plotted with the measured char- 
acteristic. An example prediction with experi- 
mental results is shown in Figure 9-63 for the 
RD 308 flip-flop.* 


Failure in MOS logic circuits. results 


rom changes.in the threshold voltage of the - 


transistors caused by ionizing radiation. Since a 
large. negative supply voltage permits greater 
degradation in threshold voltage before circuit 
failure occurs, the radiation failure threshold of 
MOS circuits depends on the maximum supply 
voltage rating. Experiments indicate that MOS 
digital] microcircuits fail at radiation levels from 
105 to 6 x 10? rads (Si) at manufacturer’s rated 
supply voltages. Such an exposure can be associ- 
ated with a neutron fluence of 1 to 2 x 10!4 
nicm? (E > 10 keV, fission) for the mixed 
neutron-gamma flux of a typical fast-burst re- 
actor. The characteristics of a particular MOS 
integrated circuit must be established with rea- 
sonable confidence before meaningful ‘predic- 
tions can be made. 

As with digital areuits. the pmmary 
cause of linear-microcircuit failure is transistor 
gain degradation. The degradation of perfor- 
mance of a linear circuit is characterized by 
radiation induced changes in the transfer charac- 
teristics. In linear circuits. the functional 


dependence of the overall circuit performance . 


on individual transistor elements can be deter- 
mined only by a detailed circuit analysis. This 
analysis usually is frustrated by circuit complex- 
ity and the inability to measure individual 
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mMicrocircuit. elements as a result of the lack of 
accessible terminals. Therefore, prediction of 
the performance of linear circuits under irradia- 
tion is difficult. The large variety of linear. cir- 
cuits precludes. a discussion of each type: how- 
ever, some general comments can be made con- 
cerning the performance of some devices. — 
Pd radiation responses of both differ- 
ential and operational amplifiers have been 
studied. A typical transfer characteristic of a dif- 
ferential amplifier is shown in Figure 9-64. The 
gain of the circuit began to decrease at a fluence 
of 3 x 10!3 n/cm? (F > 10 keV, fission) and 
was degraded to roughly 50 percent of its initial 
value after an order-of-magnitude-larger fluence. 
These amplifiers were found to maintain their 
balance. during irradiation significantly better 


~ th mplifiers made from discrete devices. 
wae The largest changes in operational ampli- 


lers induced by radiation were observed in the 
open loop voltage gain and input bias current. 
The reduction in the open loop gain is a direct ~ 
consequence of degradation of transistor gain. 
The use of lateral and substrate PNP ‘transistors 
results in a relatively low radiation tolerance of 
these amplifiers compared to logic circuits. 
These PNP transistors are widebase units that are: 
degraded at lower fluences than vertical NPN 
transistors. Changes in voltage gain begin to be 
observed (5 percent changes) at fluence levels 
near 107 nfem? (E > 0.1 MeV, fission) for 
709-ty;pe operational amplifiers. The neutron 
fluence where the voltage gain has decreased by 
50 percent is about 8 x 10!? n/em? (E > 10 
keV. fission) for units with lateral and substrate 
transistors. Amplifiers that have eliminated 
lateral and substrate transistors show improved 
performance in the presence of radiation. The 
degradation in gain that is induced by radiation 
in operational amplifiers depends on the elec- 


: Ko details are piven in the “TREE Handbook” (see 
biography). F 
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Legend 


Load-terminal I-V characteristic, experimental 
— ——oad-terminal I-V characteristic, predicted 
————— -Load line at R, =360.2 


"On' segment 
eS 7 he S 3 
$e gi =1.4 x 10° n&em* (E>10 keV, fission) 
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Figure 9-63. eS Load-Terminal |1-V Characteristic (“On Segment) 
“of RD 308 in Neutron Environment 
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tricul design. For example. it is important to 
determine whether the open-loop gain is deter- 
mined by resistor ratios or by transistor gains. 
The changes in the input bias current can be 
correlated directly with changes in the common 
base current gain of the input transistors. It 
should be noted that the input transistors of 
operational amplifiers operate at very low cur- 
rents (high input impedance), thus degradation 
of the base-transport factor is accompanied by 
degradation of the emitter efficiency as well. 
Factor-of-two increases in bias current have been 
observed for 709 amplifiers after 3 x 10!3 
niem* (£ > 10 keV, fission). Offset voltage and 
offset current were found to increase after ir- 
radiation. These changes result from emitter- 
base-voltage and current-gain mismatches after 
irradiation. The changes in both current and 
voltage offsets were small at fluences where the 
gain and the bias current were degraded appreci- 


ably, which indicates the uniformity of active 
elements on the same chip. 

Even though a large number of inte- 
grated microcircuits have been tested, data on 
the effects of neutron irradiation on microcir- 
cuits are still sparse. This is especially true for 
the linear types of integrated circuits. Radiation 


‘ experiments indicate that the failure threshold 


of digital microcircuits is fairly independent of 
the construction technique. For buffered cir- 
cuits the failure threshold, at unity fan-out, is 
near 10!5 n/em? (E > 10 keV, fission), while it 
is somewhat lower for the nonbuffered circuit 
types: however, the failure level at rated fan-out 
(~10) usually occurs over ai order of magnitude 
lower in fluence than the failure level at unity 
fan-out. Representative radiation failure levels 
for some common digital microcircuits are 
shown in Table 9-24. The failure level is speci- 
fied when the output voltage of the test circuit 


Table 9-24. 8 Failure Thresholds for Typical Digital Microcircuits & 


Failure Level, 


Designation Function Construction _ (nicm*) 

MC 20] DTL Gate Junction Isolation 1S x10" 12x07 

DT wl 932 DTL Gate ° Junction Isolation 3.0 x 105 1.5 x 1914 

RD 209 DIL Gate Oxide Isolation 3.0 x 10!° 45 x 1014 

MC 507 TTL Gate - Junction Isolation | 0.8 x 1035 0.8 x 10/4 

SN 54932 TTL Gate Junction Isolation nsx 105 4.2.x 1014 

DT ul 945 DTL Flip-Flop Junction Isolation 1.5 x 1015 1.2 x 104 

SE 124 DTL Flip-Flop Junction Isolation 0.85 x 10/5 0.8 x 104 7 
RD 208 DTL Flip-Flop Oxide Isolation 0.85 x 10!5 0.8 x rls ca 


+ 


* Failure level at fan-out of 1; neutron fluence specified as (E > 10 keV, fission). | 
Failure level at fan-out of 10: neutron fluence specified as (E > 10 keV, fission). 
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exceeds the noise margin of the following cir- 
cuit. Typical radiation failure levels for some 
linear circuits are shown in Table 9-25. Table 
9-26 contains irradiation test results for some 
representative MOS integrated circuits. The fail- 


‘ure level is specified as the point when the cir- 


cuit would not operate or when the threshold 
voltage exceeded the supply voltage. 
The test results presented are only in- 
tended to provide a broad range of failure levels 
for order-of-magnitude reference purposes and 


should be treated with appropriate caution. It 

should be borne in mind that design and produc- 

tion changes in integrated circuits are common 

in industry. For this reason, each circuit should 
sidered a separate problem. 

The thermomechanical shock éffects for 
all integrated circuit types are the same as those 
effects on. discrete semiconductor parts. The 
only difference to consider is the increased num- 
ber of bonds used in each device package, which 
increases the change of bond failure. . 


Tabie 9-25. a Failure Thresholds for Typical Linear Microcircuits F 


Designation Function 

uA 709 Operational Amplifier 
RA 909 Operational! Amplifier 
Ph 709 Operational Amplifier 
MC 1709 Operational Amplifier 
MC 1525 Differential Amplifier 
NM 1024 Differential Amplifier 
NM 1006 Differential Amplifier 
RA 138 Amplifier 


Bp: level — gain degradation $0 percent, 
t 


Neutron fluence specified as (E > 10 keV, fission). 
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Failure Level* 


~ Construction (n/cm?}t 

Junction {Isolation 0.8 x 10!4 
Oxide Isolation 0.8 x 1014 
Oxide Isolation 3.0 x 104 
Oxide Isolation . 30x 10! 
Junction Isolation 3.0 x 104 
Oxide Isolation | 3.0 x 1014 
Junction Isolation 4.0 x 10/4 
Oxide Isolation 1.5 x 10!4 


Table 9-26. @ Failure Thresholds for Typical MOS Digital Microcircuits & 


Failure Level 


Gamma, Neutron.* 
Designation Function : (rads (Si)) (n/em* ) 
SC 117] NAND gate 2x 10°+ 
MEM 529 Binary element 14x 1052 Phe - 
SC 1171 Binary element | ee oe _ 
MEM. SOF. Shift register 1 x 10°§ = 
“MEM $90 Chopper ; Not measured** a8. 70!" 
SC 1149). Flip-Flop » Not measured** ‘sx 10/4 
MC 3155  AND/OR gate 2 x 10° (Cobalt-60)++ 
3300 25-bit static shift ->5 x 10? (FXR)it 
register >8 x 10* (TRIGA)Ii 
3003. ~ - - JOO-bit shift register >2 x 104 (FXR)it 
os >5 x 104 (TRIGA)i+ 
1406 100-bit-shift register >10° (FXR)it 
—<2.5 x 10* (TRIGA) 
1101 — 256-random access memory 4x 104 (FXR) 
2x 104 (TRIGA) 3 x 10!} 
+ Neutron fluence specified as (E > 10 keV, fission). 
t Supply voltage - 20 volts. 
i Supply voltage - 15 volts. 


Clock voltage - 10 volts. 

Supply voltage - 10 volts. 

Type of facility in which test was performed. 
z No failures at these levels. 


a 2 
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SECTION VII 


Ls | ELECTROMAGNETIC PULSE (EMP) 
DAMAGE MECHANISMS 


As described in Chapter 7, the nuclear 
electromagnetic pulse (EMP) is part of a com- 
plex environment produced by a nuclear envi- 
ronment. The EMP contains only a very small 
part of the total energy produced by a nuclear 
explosion; however, under the proper circum- 
‘stances, EMP is capable of causing severe dis- 
ruption and sometimes damage to electrical and 
electronic systems at distances where all other 
effects are absent. 

As with the EMP generation described in 
apter 7, the complexity of the calculation of 
EMP damage mechanisms requires that heavy 
reliance be placed on computer code calcula- 
tions for specific problems, and even these calcu- 
lations must be supplemented by testing in most 
cases. Consequently, the information presented 
herein is largely. qualitative and will only serve as 
an introduction to the subject. More complete 
treatments of EMP damage mechanisms may be 
found in the “DNA EMP (Electromagnetic 
Pulse) Handbook” (see bibliography ). 

Figure 7-18, Chapter 7, provides a 
matnix that provides some indication of whether 
EMP constitutes a threat in a given situation 
relative to the hardness of a system to blast over- 
pressure. This section provides a brief descrip- 
tion of EMP energy coupling, component dam- 
age, EMP hardening, and testing. 


@B evercy covuriine @B 


‘9-56 Basic Coupling Modes 


There are three basic modes of coupling 
the energy contained in an electromagnetic wave 
into the conductors that make up an electric or 
electronic system: electric induction, magnetic 
induction, and resistive coupling. 

2 Electric induction arises as the charges in 
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a conductor move under the influence of the 
tangential component of an impinging electric 
field. The overall result is that of a voltage 
source distribution along the conductor. One 
such point-voltage source is shown in Figure 
9-65: for a simple conducting wire, where the 
current lis produced as a result of the tangential 


component E; ,,,, of the incident electric field 


it 


_ COPPER WIRE 


CHARGE SEPARATION 


Figure 9-65. ea Electric Induction in a 


Copper Wire be 


| Magnetic induction occurs in conductors 


shaped fo form a closed loop when the compo- 
nent of the impinging magnetic field perpen- 
dicular to the plane of the loop varies in time, 
causing charges to flow in the loop. This effect is 
illustrated in Figure 9-66 for a simple wire loop. 
Here the magnetic field is shown coming out of 
the plane of the loop. The loop need not be 
circular, and magnetic induction may occur with 
any set of conducting components assembled so 
as to form a loop. 
Resistive. coupling comes about indi- 
rectly as a conductor that is immersed in a con- 
ducting medium, such as ionized air or the 
ground, is influenced by the currents induced in 
the medium by the other coupling modes. In 
effect the conductor shares part of the current 


‘as an alternate conducting path. This effect is 


illustrated in Figure 9-67 for the simple case of a 


. 


LOOP 
ANTENNA 


1 


Figure 9-66. a Magnetic Induction in a 
Simple Loop 


conductor immersed in the ground. The tangen- 
tial component of the incident electric field £, 
induces a current density J in the ground. A 


distributed voltage drop appears along the wire ~ 


as a result of the current flow in the ground, and 
this incremental voltage causes current flow | in 
the wire. Current also may be induced in the 
wire directly by the tangential component of the 
refracted electric ‘field, shown as £.. The re- 
flected EMP, E.. A is also shown in Figure 
9-67. The potential importance of these ‘re- 
flected fields is discussed below. 


9-57 Resonant Configurations 


The coupling of energy to a conductor is 
particularly efficient when the maximum dimen- 
sion of the conductor configuration is about the 
same size as the wavelength of the radiation. In 
this event the voltages that are induced along the 
conductor at various points are all approximate- 
ly in phase, so the total voltage induced on the 
conductor is a maximum. The conductor is said 
to be resonant, or to behave as an antenna, for 


Av- WIRE 
+ : | y) 
AZ 


Figure 9-67, @ Resistive Coupling as a Result 
of Currents in the Ground 


frequencies corresponding to near this wave- 
length. Since EMP has a broad spectrum of fre- 
quencies (see Chapter 7), only a portion of this 
spectrum will couple most efficiently into a 
specific. conductor configuration. Thus, a par- 
ticular system of interest must be examined with 
regard to its overall configuration as well as its 
component configuration. Eash aspect will have 
characteristic dimensions that determine what 
part of the pulse (strength and frequencies) con- 
stitutes the principal threat. 

Gross system features that are not nor- 
ma considered antennas, such as structural 
features, beams, girders, buried cable, overhead 
conduit or ducting, wings. fuselage. missile skins, 
and any wall apertures, must be. considered to be 
potential collectors and conductors of energy 
into the system. In particular, radiation that 
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enters through an aperture is analogous to radia- 
tion that originates from a plate of the same size 
and shape as the aperture. Thus, it is resonant, 
and the aperture is resonant, and it admits a 
maximum of energy from the pulse for those 
frequencies near its resonance. 

When the EMP strikes the ground, a por- 
tion of the pulse will be transmitted through the 
. interface, inducing currents in the ground or any 
system components buried there, and a portion 
normally will. be reflected as shown in Figure 
9-68. Thus, a system that is above the ground 
will receive the reflected pulse as well as the 
direct pulse. These may cancel one another par- 
tially, but in the worst case they may’ reinforce 
and may constitute a greater threat level. 

It can be seen that most practical sys- 
ems in their operational environment present 
exceedingly complex coupling problems for an 
arbitrary explosion. The solution for any com- 
bination of svstem and environment is probably 
unique and will be very sensitive to even minor 


Ej Er 


Figure 9-68. Reflected and Refracted Waves 
at the Air-Ground Interface 
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changes in the parameters. Two approximate 
approaches have been tried: computer studies as 
mentioned in Chapter 7, and threat simulation, 
which will be discussed in succeeding 
paragraphs. 


9-58 Types of Damage 


BR component vamace Gy 
a 


Degradation of system performance may 
occur as a result of functional damage or opera- 
tional upset. A system will suffer damage if it is 
damaged permanently as a result of a large elec- 
trical transient. For example, a catastrophic fail- 
ure of a. device or component will render its 
eperation unsatisfactory in any circuit. A para- 
metric failure of a device occurs when degrada- 
tion of some parameter has proceeded to a point 
where the circuit will continue to operate but at 
reduced efficiency. These latter failures are 
classed as functional damage. On the other hand, 
a. system suffering operational upset is only im- 
paired temporarily. 

= Electronic components that are sensitive 

to functional damage or burnout are listed 
below in the order of decréasing sensitivity to 
damage effects: 


" microwave semiconductor diodes, 
field-effect transistors, 
radio-frequency transistors, 

~ audio transistors, 
silicon-controlled rectifiers, 
power rectifier semiconductor diodes, 
vacuum tubes. | 


Thus, on the basis of components alone, yacuum 
tubes are less susceptible to EMP damage effects 


than transistors. . 
Electronic or electrical systems that are 


subject to malfunction include: 
Most susceptible: 
e@ Low power, high speed digital computer 


(upset) either transistorized or vacuum 
tube 


Systems employing transistors or semicon- 
ductor rectifiers (eitner silicon or selen- 
ium), such as 


computers 
computer power supplies 


transistorized components terminating long 
cable runs, especially between sites 


alarm systems 
intercom systems 
life-support system controls 


some telephone equipment which is par- 
tially transistorized 


transistorized receivers 
transistorized transmitters — . 
| transistorized 60 to 400 cps converters 
transistorized process control systems 
power system controls: communication 
links 

Less ‘susceptible: 


Vacuum tube equipment (does not include 
high speed digital equipment and equip- 
ment with semiconductor or selenium récti- 
fiers) 


‘transmitters intercoms 


receivers teletype-telephone 


alarm systems power supplies 


Equipment employing low. current 
switches. relays. meters 


panel indicators, 
status boards 


alarms 


life-support svstems process controls 


power system control ~ 
panels 


@ Hazardous equipment.containing 


detonators explosive mixtures 
squibs rocket fuels 
pyrotechnical 

_ devices 

® Other 


Long power cable runs employing dielectric 
insulation, equipment associated with high 
energy storage capacitors or inductors: 


Least susceptible: 
@ High-voltage 60 cps equipment 


transformers. 
motors 


rotary converters 


heavy duty relays, 


lamps, filament 
circuit breakers 


air-insulated power 
cable runs 


heaters 


The Jess susceptible equipment or com-. 
ponents would be made more susceptible if they 
are connected to long exposed cable runs, such 
as intersite wiring or overhead exposed power or 
telephone cables. The equipment can be made 
less vulnerable if it is protected. 


9-59 Damage Levels Pod 


The nature of a circuit has a strong bear- 
ing on the transients that cause damage; how- 
ever, in general pulse lengths of microsecond and 
submicrosecond duration are required to cause 
problems. Table 9-27 shows a list of common 
active devices and the approximate energy re- 
quired to cause functional damage. The wide 
range of energies should be noted. 

The minimum energy required to 
damage meters or to ignite fuel vapors is about 
the same as that required to damage semicon- 
ductors as shown in Table 9-28. Good composi- 


_tion resistors can withstand pulse powers more 


than 10,000 times their power rating for micro- 
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Table 9-27, eB Minimum Observed Joule Energy to Cause Burnout = 


Minimum 
Joule 
Type Energy Material Other Data 

2N36 4.0 x 10? Ge PNP Audio Transistor 

2N327A 1.6 x 10? Si PNP Audio Transistor 

2N1041 2.0 x 10? Ge PNP Audio Transistor 

2N1308 5.0 x 30° Ge NPN Switching Transistors 

2N706 6.0 x 10° Si -NPN Switching Transistors 

2N594 6.0 x 10° Ge NPN Switching Transistors 

2N398 8.0 x 104 Ge PNP Switching Transistors 

2N240 1.0 x 10°? Ge PNP Switching Transistors 

MC715 8.0 x 10° Si Data Input Gate Integrated Circuit 
- 2N4220 1.0 x 105 Si RF General Purpose. FET 

2N4224 3.0 x 10° Si _ VHF Amp and Mixer FET 
1N3659 8.0 x 107 Si Automotive Rectifier Diode 

1N277 2.0 x 10° Ge High Speed Switching Diode 
1N3720 50:x 107 Tunnel Diode . 

1N238 1.0 x 107 Si Microwave Diode 

2N3528 3.0 x 10° Si Silicon Controlled Rectifier 
67D-5010 1.0 x 10% G.E. Varistar (30-joule Rating) 
6AF4 1.0 x 10° UHF Oscillator Vacuum Tube | 
66N8 2.0 x 10° General Purpose Triode Vacuum Tube 
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Table 9-28. Minimum Joule Energy to Cause Permanent 
" Degradation Indicated 


Minimum 
Joule 

Designation Energy Malfunction Other Data 

Relay 2x 107 Welded Contact Potter-Brumfield (539) 
low-current relay 

Relay 1x jo! Welded Contact ' Sigma (JIF) one-ampere 
relay 

Microammeter 3.x 10° Slammed Meter Simpson Microammeter 
(Model 1212C) 

Expivsive Pe aoe Ignition EBW 8 amp for 10 usec 

Bolt detonator, MK1 

Squib 2x 107 Ignition Electric Squib. N& 
3.5 watts for 5 usec detonator 

Fuel Vapors 3x 10° Ignition Propane-air mixture 


1.75 mm ignition gap 


second pulses. Capacitors are also fairly hard 
components. The approximate energies required 
for degradation of several common components 
are shown in Table 9-28. 

The minimum energy necessary for oper- 
ational upset is a factor of 10 to 100 less than 
that which is required to damage the most sensi- 
tive semiconductor component. Table 9-29 
shows the levels required to cause operational 
upset to some common components to illustrate 
this factor. 


quired to damage several classes of electrical 
equipment is provided in Figure 9-69. 

The large range of damage levels empha- 
sizes the fact that it is important to consider 
EMP damage criteria early during the design 
stage of any piece of equipment that might be 


A gross comparison of the energy re- 


susceptible. It is also important to realize that 
energy collected in one part of a system may be 
transmitted to other parts of the system as a 
result of the currents that are induced. Thus, it 
is not necessary that the EMP couple directly to 
a sensitive component: energy coupled to vari- 
ous parts of a system may ultimately reach a 
particular component in sufficient quantity to 
cause malfunction. With the current state of the 
art in EMP vulnerability evaluation, the design 
and hardening of complicated systems requires 
the joint efforts of systems engineers and profes- 
sional EMP effects personnel. 


WB cme narpenine 


9-60 System Analysis 
r | A general approach to the examination 
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Minimum Joule Energy to Cause Circuit 


of a system with regard to its EMP vulnerability 
could include the following steps. First informa- 
tion concerning the system components and de- 
vices is collected. The information is categorized 
methodically into physical zones based on the 
susceptibility and worst case exposure for these 
items. Using objective criteria, problem areas are 
identified, analyzed, and tested. Suitable 
changes are made as necessary to correct de- 
ficiencies, and the modified system is examined 
and tested. The approach may be followed on 
proposed systems or those already in place, al- 
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Table 9-29. 
Upset or Interference 
Minimum 
Joule 
Designation Energy Malfunction Other Data 
Logic Card 3 x. 10% . Circuit Typical logic transistor 
Upset inverter gate 
Logic Card 1x 10° Circuit Typical flip-flop transistor 
Upset assembly 
Integrated 4x 107°... Circuit Sylvania J-K flip-flop -monolithic 
Circuit Upset ‘integrated circuit (SF50) 
Memory Core 2x 10° Core Erasure Burroughs fast computer core 
Via Wiring memory (FC2001) 
' Memory Core sx 10° Core Erasure Burroughs medium speed computer 
Via Wiring core memory (FC8001) 
Memory Core 3x 10° : Core Erasure RCA medium speed, core memory 
Via Wiring (269M 1) 
Memory Core 2x 108 Minimum observable energy in a 
typical high-gain subsystem 
Amplifier 4x 107! Interference Minimum observable energy in’ a 


typical high-gain amplifier 


though experience indicates that the cost of 
retrofitting EMP protection is usually over- 
whelming. 


9-61 Recommended Practices oe] 


Within the scope of this manual it is 
only possible to mention a few of the practices 
that may be employed in hardening a system to 
EMP. The following discussion is intended to 
convey some impression of the extra effort 
involved in hardening a system to the EMP 
rather than to provide a comprehensive treat- 
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Energy Required to Damage Various 


Classes of Equipment 


ment of what is a highly technical and special- 
ized field. a 
Some general methods for reduction of 
the EMP environment include geometric arrange- 
ment of the equipment, shielding, geographic 


relocation, and proper grounding. 
<a Circuit layout recommendations include 
1e use Of common ground points. twisted cable 
pairs, system and intrasystem wiring in “tree” 
format (radial spikes) avoiding loop layouts and 
circuit routes coupling to other circuits, use of 
conduit or cope trays, and shielded-.isolated 
transformers. Avoiding ground return in cable 
shields is also recommended. Many specific prac- 
tices carry over from communications and 
power engineering while many do not. Each 


must be examined carefully. 
Ca Good shielding practices include the use 


of independent zone shields, several thin shields 


to replace a thick one, continuous shield joints, 
and keeping sensitive equipment away from 
‘shield corners. Avoiding shield apertures, and 
avoiding the use of the shield as a ground or 
return conductor is also recommended. The 
shielding effectiveness of many enclosures fre- 
quently is defeated by energy carned by cables 
or pipes (including water pipes, sewage lines. 


etc.) into the enclosure. 
Cabling recommendations include’ the 


use of deeply buried intersystem cables (more 
than 3 feet), shield layer continuity at splices, 
and good junction box contact. Ordinary braid 
shielding should be avoided. Cable design. repre- 
sents an extension of shielding and circuit prac- 
tices from the viewpoint of EMP protection. It is - 
an area where compromises frequently are made 
in the interests of economy, and thus is an area 
where professional EMP effects personnel can be 


9-177 


of considerable assistance. 

Good grounding practices must be em- 
ed. In general, a “ground” is thought of as a 
part of a circuit that has a relatively low im- 
pedance to the local earth surface. A particular 
ground arrangement that satisfies such a defini- 
tion may not be optimum, and may be worse 
than no ground from the EMP viewpoint. A 
ground can be identified as: the chassis of an 
electronic circuit, the “low” side of an antenna 
system, a common bus, or a metal rod driven 
into the earth. The last depends critically on 
local soil conditions, and it may result in resis- 
tive induced currents in the ground circuit. A 
good starting point is to provide a single point 


ground for a circuit cluster, usually at the lowest ~ 


impedance element — the biggest piece of the 
system that is electrically immersed in the earth, 
e.g., the water supply system. It is beyond the 
scope of this manual even to list all of the 
grounding recommendations. Once again, this is 
an area where professional EMP effects per- 
el can be of considerable assistance. 

Finally, various protective devices repre- 
_ sent a means to counter other protective short- 
comings indirectly. Filters, absorbers, limiters, 
decouplers, switching devices, arc arresters, 
fuses, etc., are part of this class of components. 
When other design practices cannot be used or 
are not adequate, such devices must be added. 
Typically they are found in an “EMP room” at 
the cable entrance to underground installations, 


in aircraft antenna feeds, in telephone lines, at. 


power entry panels to shielded rooms, etc. On 4 
smaller scale, diodes, nonlinear resistors, SCR 
clamps, and other such items are built into cir- 
cuit boards or cabinet entry panels. Few of these 
devices by themselves are sufficient as a com- 
plete solution to a specific problem area, be- 
cause each has some limitation in speed of 
response, voltage rating, power dissipation 
capacity or reset time. Thus, most protective 
devices are hybrids. A few prepackaged hybrids 
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are available for the protection of audio and 
power lines from lightning strokes and power 
surges and, if modified, may be used for EMP 
protection. No such packages are readily avail- 
able for high frequency lines, multiple wire 
cables, antennas, etc., and usually must be 
custom designed for each application. 


GP testinc 


9-62 {importance of Testing a 


Even with present day sophistication in 
analytic techniques, it is clear from the complex- 
ities described above that sole reliance cannot be 
placed in analysis and prediction alone. Testing 
has a.number of important functions. . 

Testing is essential to verify prior anal- 
ysis of devices, components, and complete sys- 
tems early in the design stage. Testing also is the 
only known method that can be used to identify 
surprises. Surprises can be unexpected coupling 
or interaction modes or weaknesses that were 
overlooked during the design. Nonlinear effects 
in interaction are a form of surprise that only 
can be found by testing. After the test, many of 
the original approximations made in analysis can 
be refined and improved for future analysis, and 
the data can improve the analytic capability for 


m omplex problems. - 
“ Testing quickly locates weak or sus- 


cepuble points in components or systems early 


enough for economic improvement. After the 
improvements, testing quickly verifies that the 
improvements bring the performance up to 
standards. . 


Testing provides assurance and con- 
idence that the complete system is actually 
hardened to EMP to the specified threat level. 
Actual certification can only be obtained by 
providing the actual nuclear threat environment. 
Further, periodic testing insures that system 
hardness is not degraded as a result of environ- - 
mental or human factors. 


9-63 Simulation Facilities . 


As a result of the limited test ban treaty. 
heavy reliance must be placed on simulation to 
test the EMP hardness of systems. A brief de- 
scription of generic simulation techniques is 
given below. A more thorough description of 
these techniques is contained in the “DNA EMP 
(Electromagnetic Pulse) Handbook, Volume 2, 
: Ivsis and Testing” (see bibliography). 

The classes of EMP tests are: 

@ Low level current mapping. 
@ High level current mapping. 
® High level fields. 
Low level current mapping is a good test for the 
‘beginning of any program. With the system 
power turned off and a low-level field. the mag- 
nitudes and signatures on internal. cables are 


determined. This provides an insight .on the- 


work that must follow. After this test and the 
indicated. improvements are made. a high-level 
current can be injected directly into the system 
with the system power on to explore. for non- 
linearities, and to uncover initial indications of 
system effects. If subsystems malfunction, it 
may be desirable to conduct extensive subsys- 
tem tests in the laboratory. Finally, a high level 


field test is essential. 
The type. of excitation must be defined 


in any type of test. The two principle choices 
are: 
@ Waveform simulations (time-domdin data), 
® Continuous wave -signals (frequency- 
domain data). 
If the intent is to match a system analysis in the 
frequency domain to measured system response, 
continuous wave (CW) signals may be more suit- 
able. If the desire is to compare the test results 
to known electronic thresholds, it is frequently 
necessary to test in the time domain. For a com- 
plete analysis, it is advisable to consider both 
ty pes of tests. 
ze Large-scale simulators are required for 


the final test. of large systems. The two principal 
kinds of large simulators are: 
@ Metallic structures that guide an EM wave 
past a test object, 
@ Antennas that radiate an EM field toa test 
object. 
Guided wave simulators use pulse generators 


' that simulate EMP waveforms and operate in the . 


time domain. Radiating antennas use either 
pulse generators (time domain) or continuous 
wave (CW) signal generators (frequency do- 
main). Pulse generators themselves can be either 
high level single shot or low level repetitive. 

(U) The essential elements of a guided wave. 
or transmission line simulator include: 


@ Pulser 

@ Transition section 
@ Working volume - 
e Termination. 


‘An electromagnetic wave of suitable amplitude 


and waveshape is generated by the pulser. This 
wave is guided by a tapered section of transmis- 
sion line (the transition section) from the small 
cross sectional dimension of the pulser otuput to 
the working volume. The working volume, 
where the test object is located, should be large 
enough to provide a certain degree of field uni- 
formity over the test object. A test object di- 
mension one-fifth that of the working volume 
satisfies this condition. The termination region 
prevents the reflection of the guided wave back 
into the test volume, and consists of a transition 
section that guides the incident wave to a geo- 
metrically small resistive load whose impedance 
is equal to the characteristic impedance of the 
transmission line structure. 

The basic types of radiating simulators 
are: 

@ Long wire 
® Biconic dipole or conical monopole. 

The long wire is usually a long dipole oriented 
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parallel to the earth’s surface. It is supported 
above the ground by nonconducting poles with 
high-voltage insulators. The two arms of the di- 
pole are symmetric about the center and con- 
structed from sections of light weight cylindrical 
conductor, such as irrigation pipe. Pipe sections 
decrease in diameter with increasing distance 
from the center, and resistors are placed be- 
tween the pipe sections to shape the current 
wave and reduce resonances. The two arms of 
the dipole are oppositely charged, and. when the 
voltage across the spark gap at the dipole center 
reaches the breakdown voltage, the gap begins 
conducting and a current wave front propagates 
awav_ from the gap. 


pulsers, such as Marx generators,-or CW trans- 
mitters instead of relying on the discharge of 
static surface charges. The antennas are fabri- 
cated out of lightweight conducting surfaces or 


wire grids. 
- Differences between guided wave and 


radiating simulators are listed in Table 9-30. 
Electromagnetic scale modeling is an im- 
portant alternative to full scale testing under the 
following conditions: 

@ Test facilities or available equipment are at 
a premium, 

® The system to be tested is very large, 

@ The system dedication cost for full scale 
testing is high. 

In addition to the advantages of modeling under 
these conditions, benefits can be derived as for 
lows: 

@ Perhaps sensors can be placed better during 
full scale testing as a result of model 
experiments, 

® Design modifications or cable reroutes can 

' be made prior to full scale testing, 

e EM angles-of-arrival can be determined for 
worst- and best-case conditions, 

@ Effects of changes in the conductivity of 
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Conical and biconical antennas use — 


the surrounding media can be explored to 
an extent, 

@ Estimates can be made of some of the 
responses of a complex system prior to full 
scale testing, 

@ Design confirmation of costly systems can 
be made prior to system fabrication and 
costs can be reduced, 

@ Quantitative data can be obtained to vali- 
date analysis. 

= It should be pointed out that because of 
the difficulty in introducing minute openings or 
poor bonds into models, and since these often 
control interior fields, the usefulness of model- 
ing ordinarily is limited to the measurement of 
external fields, voltages, and currents. Once the 
exterior fields, voltages, and currents are known 
for a complex structure, perhaps having cable 
runs, analysis can often provide internal field 


uantities.of interest. 
= In actually setting up a scale model test, 
he fol 


lowing should be kept in mind: . 

@ Broadband pulse response determination 
involves much more than a steady-state, 
single-frequency response test does. 

@ Special EM illumination sources that are 
coherent, have uniform time delay, and use 
antennas with constant effective height are 
required. 

® Special modeling techniques are required to 
study exposed conductors that pass over or 
within a lossy dielectric, such as earth. 

A pulse-type waveform theoretically can ‘be re- 
placed by a continuous wave (CW) source with a 
sensing system which references the sensed CW 
signal to a reference phase from -the source. 
Complex Fourier transfer functions can be de- 
veloped by processing the recorded data on a 
computer. However, long sweep times are re- 
quired to ensure that all narrow band responses 
are explored adequately, and the actual physical 
implementation of such an approach in the 


Table 9-30. | Comparison of Guided Wave and Radiating Simulators Qj 


Guided Wave 


Radiating 


Energy use 


Test volume 


Polarization 


Angle of incidence 
Earth reflection effects 


Geometric. attenuation 
of EM wave amplitude 


Efficient — mostly 
directed toward 
test object 


Limited by size of 
simulator — difficult 
to construct large 
enough for. sizable 
test objects 


Fixed (or bipolar, 


e.g., ARES) 
Fixedt 
No 


No — relatively 
uniform within 
test volume 


Energy radiated 
symmetrically about 
axis — only fraction 
directed toward 

test object 


Limited by desired 
field intensity 


Variable* 


_ Variable* 


Yes 


Yes (1/R) 


Planar wave capability Yes Only at distance 
Interference with nearby Limited Yes 

electronics 

* ‘hese are, however, limited by the height of the antenna unless it is airborne. 


+ Polarization is fixed relative to earth coordinates; however, a range of polarizations and angles of incidence can be 
provided in some facilities by changing the position and orientation of the object that is being tested. For example, a 
missile can be rotated on several axes in ARES to change these to items relative to the missile. 
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microwave band poses additional difficulties. On 
the other hand, the use of scaled real time wave- 
forms allows quick development of actual re- 
sponses, from which complex Fourier transfer 
functions also can be developed with the aid of 
computers. , 


Several variations of each type of simula- _ 


tion technique described above are currently 
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operational. Each has some advantages and dis- 
advantages when compared to others. As men- 
tioned previously, it is beyond the scope of this 
manual to describe. individual facilities. The 
interested user should consult the “DNA EMP 
(Electromagnetic Pulse) Handbook, Volume 2, 
Analysis and Testing,” (see bibliography), and 
the references listed therein. 
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